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HAUPTAUFSÄTZE 


Theorie der stationären Leewellenströmung in freier 
Atmosphäre*). 
Von @. Lyra in Göttingen. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


Ziel der vorliegenden Arbeit ist die Herleitung der auf der Leeseite von 
quer überströmten Höhenzügen beobachteten stationären Wellenströmung 
einer stabil geschichteten Atmosphäre aus den (irundgleichungen der Me- 
chanik. Das Problem läßt sich auf die erste Randwertaufgabe der Wellen- 
gleichung für die Halbebene zurückführen. Die eindeutige Bestimmung der 
durch die Natur realisierten Lösung dieser an sich unbestimmten Rand- 
wertaufgabe bildet den Gegenstand der malhematischen Untersuchung des 
zweiten Teils, wo drei verschiedene Methoden zur Lösung des Eindeutigkeits- 
problems dargelegt werden. Die vollständige Lösung der Aufgabe wird im 
dritten Teil auf Profile spezieller Bodenformen angewandt und mit den 
bisher vorliegenden Beobachtungen verglichen. 


Einleitung. 


Während der letzten Jahre hat L. Prandtl verschiedene Probleme der Mechanik der 
freien Atmosphäre mit der sonst in der angewandten Mechanik üblichen Strenge der Methode 
von neuem in Angriff genommen und einige Ergebnisse in der Neuauflage seines „Abriß der 
Strömungslehre“* mitgeteilt. Aus der Prandtlschen Schule ist auch die vorliegende Unter- 
suchung hervorgegangen. Sie hat jene stationäre, wellenförmige Strömung der freien Atmo- 
sphäre zum Gegenstand, die sich auf der Leeseite von quer überströmten Höhenzügen, Plateau- 
hängen und anderen langgestreckten Bodenerhebungen ausbildet (vgl. Bild 1). Indem „Wellen- 
berge“ dieser leeseitigen Wellen über die Kondensationshöhe hinausragen, wird in Lee des 
Gebirgshindernisses diese stationäre Luftströmung als stehende „Moazagotl-Wolke“ in Form 
einer oder mehrerer zum Gebirgskamm paralleler Wolkenbänke sichtbar, welche fortwährend 
am Luvrand ihrer Unterseite entstehen und am Leerand ihrer Unterseite vergehen. Grunauer 
Segelflieger haben 1933 die in Lee des Riesengebirges über dem Hirschberger Tal auftretende 
Wellenströmung der Atmosphäre als neue Segelflugmöglichkeit entdeckt und 1938 im „Moazagotl- 
Aufwind“ dieser Strömung mit dem Segelflugzeug Höhen bis zu 8600 m über Meer erreicht’). 


* Erweiterte Wiedergabe eines Vortrages, gehalten auf der Tagung der Gesellschaft für angew. Math. u. Mech. 
in München am 13. Juni 1942. 

1) Reiches Beobachtungsmaterial auf Grund segelfliegerischer Erfahrung und ausführlicher Literaturnachweis 
findet sich bei J. Küttner: Moazagotl und Fühnwelle. Beitrag z. Phys. d. fr. Atmosphäre, Bd. 2 " S. 79 bis 114 
u. S. 297. 
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Bild 1. 


Das Phänomen der Leewellen tritt bekanntlich auch schon an der freien Oberfläche eines 
Flüssigkeitsstromes unter dem Einfluß stationärer Störungen auf. Die Wellen, die in Lee 
eines auf die Oberfläche einer homogenen, inkompressiblen Strömung ausgeübten stationären 
Druckes oder — was auf dasselbe hinauskommt — hinter einer entlang der Oberfläche einer 
ruhenden Flüssigkeit gleichförmig fortschreitenden Druckstörung auftreten, sind zuerst von 
Lord Rayleigh?) 1883 für den Fall großer Tiefe theoretisch erklärt worden. Dieselbe Er- 
scheinung bildet sich aus, wenn eine inkompressible und homogene schwere Flüssigkeit von 
endlicher Tiefe eine langgestreckte Bodenerhebung quer überströmt. Die Theorie dieser 
stationären Wellenströmung in Lee einer schmalen Unebenheit des Bodens hat Lord Kelvin?) 
1886 in einer Reihe von Arbeiten entwickelt und damit u. a. eine 50 Jahre ältere praktische 
Beobachtung von Scott Russel‘) erklärt. Noch auffallender ist eine ähnliche, in den nor- 
wegischen Fjorden bekannte Beobachtung, wonach Schiffe in der über dem Salzwasser lagern- 
den seichten Süßwasserschicht bei langsamer Fahrt einen abnormen Widerstand erfahren, der 
von den einseitigen internen Grenzflächenwellen herrührt. Diese Erscheinung des „Totwassers“ 
hat V. W. Ekman?°) eingehend untersucht und H. Lamb‘) durch seine theoretische Be- 
handlung der Wellen 'in der Strömung zweier übereinander geschichteter, inkompressibler 
Flüssigkeiten verschiedener Dichte erklärt. Eine Untersuchung der Leewellenströmung in 
einem kompressiblen und kontinuierlich geschichteten Medium wurde zuerst 1940 vom Ver- 
fasser in den Beiträgen z. Phys. d. fr. Atmosphäre”) mitgeteilt. Doch ist daselbst nur eine 
erste Näherung für den Fall isothermer Atmosphäre enthalten und noch keine Lösung des 
Eindeutigkeitsproblems gelungen. 


Schon in dem einfachsten Fall der reibungsfreien, inkompressiblen und homogenen 
schweren Flüssigkeit mit freier Oberfläche besteht für die mathematische Behandlung der 
Leewellenströmung die Schwierigkeit, daß für Anströmgeschwindigkeiten, welche kleiner als 
die Grundwellengeschwindigkeit sind, das Problem unbestimmt wird. Denn man kann jeder 
partikulären Lösung des Problems wegen der unbegrenzten Ausdehnung des Stromes einen 
Zug freier Wellen überlagern, deren Geschwindigkeit derjenigen des Stromes entgegengesetzt 
gleich ist. Rayleigh hat a.a.O. diese Schwierigkeit, aus der unendlichen Mannigfaltigkeit 
möglicher Lösungen die eine durch die Natur eindeutig bestimmte Lösung auszusondern, da- 


2) Lord Rayleigh: The Form of Standing Waves on the Surface of Running Water. Proc. London Mathem. 
Soc. Bd. 15 (1883), S. 69 bis 78, oder in Papers Bd. 2, S. 258. Vgl. auch die Darstellung bei H. Lamb: Lehrbuch der 
Hydrodynamik, 2. Aufl., Leipzig 1931, S. 446 ff. 

s) W. Thomson: On Stationary Waves in Flowing Water. Part III. Phil. Mag. (5) Bd. 22 (1886), S. 517 bis 530, 
oder in Papers, Bd. 4, S. 282. Vgl. auch die Darstellung bei Lamb: Lehrb. d. Hydr., 2. Aufl., Lpz. 1931, S. 459 ff. 

% Scott Russel: Experimental Researches into the Laws of certain Hydrodynamical Phenomena... Trans. 
Roy. Soe. Edinburgh. Bd. 14 (1840), S. 79. 

5) V. W. Ekman: On Dead-Water, Scientifie Results of the Norwegian North Polar Expedition 1893 bis 1896, 
15. Teil, Christiania 1904. 

6) H. Lamb: On Waves due to a Travelling Disturbance, with an Application to Waves in Superposed Fluids. 
Phil. Mag. (6), Bd. 31 (1916), S. 386 bis 398. 


?) G. Lyra: Über den Einfluß von Bodenerhebungen auf die Strömung einer stabil geschichteten Atmosphäre. 
Beitr. z. Phys. d. fr."Atmosphäre, Bd. 26 (1940), S. 197 bis 206. 
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durch behoben, daß er zunächst durch Einführung gewisser Reibungskräfte eine allgemeinere, 
aber eindeutig bestimmte Aufgabe gestellt hat, um dann in deren eindeutig vorhandenen Lösung 
mit dem Reibungskoeffizienten den Grenzübergang zu Null zu vollziehen. Im Gegensatz zu 
diesem Verfahren, das in einer „Limitierung“* gewisser divergenter Integrale besteht, hat 
Kelvin die Eindeutigkeit einfach durch die physikalisch begründete Zusatzforderung her- 
gestellt, daß in der wirklichen Lösung stromaufwärts in hinreichender Entfernung von der 
Störungsstelle ale Wellen gelöscht sind. Auch der wesentlich allgemeinere Fall eines reibungs- 
freien, kompressiblen und kontinuierlich geschichteten Mediums führt auf das Problem der 
Eindeutigkeit, das naturgemäß weniger elementare Hilfsmittel zu seiner Lösung erfordert. 
Das geht schon daraus hervor, daß die theoretische Behandlung der Wellen in zwei überein- 
ander geschichteten inkompressiblen Flüssigkeiten verschiedener Dichte bei dem Grenzüber- 
gang zu verschwindendem Dichteunterschied versagt. Auch kann die Kompressibilität der 
aus transversalen Schwerewellen und longitudinalen Elastizitätswellen zusammengesetzten 
atmosphärischen Leewellen nicht mehr vernachlässigt werden, wenn man über eine mehr 
qualitative Betrachtung hinaus zu allgemeinen quantitativen Aussagen gelangen will. 

Es geht also um eine befriedigende Erklärung der beobachteten Strömungserscheinung 
der freien Atmophäre, d. h. um eine strenge Herleitung der stationären Leewellenströmung 
in einem kompressiblen und kontinuierlich geschichteten Medium mit beliebig linearem 
Temperaturfeld aus den Grundgleichungen der Mechanik und Thermodynamik. Für den Fall 
der isothermen Atmosphäre wurde die strenge Lösung 1941 von H. Stümke (z. Z. im Felde) 
in einem noch unveröffentlichten Bericht aus dem Kaiser- Wilhelm-Institut für Strömungs- 
forschung gegeben, wo die Eindeutigkeit des Problems durch das Rayleighsche Limitierungs- 
verfahren erreicht wird. Abschnitt 7 enthält im wesentlichen die Darstellung der Stümke- 
schen Arbeit in Übertragung auf den Fall polytroper Schichtung. Für diesen allgemeinen 
Fall einer Atmosphäre mit beliebig linearem Temperaturfeld hat der Verfasser die Theorie 
der atmosphärischen Leewellen auf einem Wege entwickelt, welcher in gewissem Sinne als 
eine „Methode der unendlich vielen Variablen“ bezeichnet werden kann. Diese Theorie ist 
in Abschnitt 6 dargestellt. Im 8. Abschnitt ist noch eine überaus kurze Herleitung der atmo- 
sphärischen Leewellen auf ganz anderem Wege gegeben. Schließlich wird im III. Teil die 
Theorie noch soweit ausgeführt, daß sie eine unmittelbare Anwendung auf praktische Auf- 
gaben der Meteorologie "gestattet. Mit der Beantwortung der Frage nach der Eindeutigkeit 
der Leewellenlösung und ihrer Erweiterung auf beliebig polytrope Schichtung ist die Theorie 
der atmosphärischen Leewellenströmung zu einem vorläufigen Abschluß gekommen °). 


I. Zurückführung des Leewellenproblems auf eine Randwertaufgabe 
der Wellengleichung. 


1. Die Zustandsänderungen in der Atmosphäre. Einem Prandtlschen Grundsatz ent- 
sprechend, soll die Aufgabe zunächst auf eine nur noch das Wesentliche des zu berechnenden 
Vorgangs enthaltende Form vereinfacht und das hierbei verbleibende Problem der reinen 
Mechanik streng formuliert und gelöst werden. Bei der atmosphärischen Leewellenströmung 
handelt es sich um ein vorwiegend zweidimensionales Problem der Hydrodynamik kom- 
pressibler und stabil geschichteter Medien, bei welchem Reibung, Erddrehung, Wärmekonvektion, 
Strahlung und Feuchtigkeit keine wesentliche Rolle spielen. Da ferner der Vorgang in erster 
Näherung als Störung einer konstanten horizontalen Gründströmung mit linearem Temperatur 
feld aufgefaßt werden kann, so bleibt der Zusammenhang mit der Erfahrung gewahrt, wenn 


8) In der meteorologischen Literatur liegen zahlreiche zum Teil sich widerstreitende Versuche zu einer Erklärung 
der atmosphärischen Leewellen vor. Doch besteht keiner dieser Versuche in einer einheitlichen Zurückführung des 
Gesamtbildes der Leewellen auf die Erhaltungssätze der Physik. Erwähnßsei z. %. die von J. Küttnera.a.O. gegebene 
Erklärung, welehe von der Annahme ausgeht, daß sich die beobachtete Leewellenströmung nur bei Erreichung einer 
bestimmten kritischen Geschwindigkeit ausbildet. Da aber die kritische Wellengeschwindigkeit für kompressible und 
kontinuierlich geschichtete Medien auch schon bei geringer Stabilität weit oberhalb der praktisch in Betracht 
kommenden Windgeschwindigkeiten gelegen ist, nimmt Küttner eine interne Grenzfläche an, deren kritische 
Wellengeschwindigkeit bei geeigneter Annahme über Dicke der „Inversionsschicht“ u. a. die gewünschte Größenordnung 
hat. Gegen die behauptete Notwendigkeit solcher Grenzflächen, die keineswegs aus dem aerologischen Material 
hervorgeht, spricht nicht nur die Tatsache, daß nach vorliegender Theorie die Leewellenströmung auch bei durchaus 
stetiger Temperaturschichtung aus den hydrodynamisch-thermodynamischen Grundgleichungen gefolgert werden kann. 
Nach Küttners Erklärung müßte schon eine geringe Überschreitung der kritischen Windgeschwindigkeit ein 
völliges Erlöschen der Leewellen zur Folge haben, was noch nie beobachtet wurde. Vielmehr treten Leewellen bei 
jeder Windstärke auf. Im übrigen beschränken sich Küttners Beobachtungen nur auf die „Moazagotl- Wetterlage‘‘, 
d. h. einseitig auf diejenigen Fälle von Leewellenströmung, welche vun der Erscheinung der Leewellenwolke 
begleitet sind. Hierdurch entstand das Mißverständnis, wonach die Leewellenströmung eine typische Föhnerscheinung 
darstellen soll. Es sei noch auf eine dem Verf. erst nach Abschluß der vorliegenden Arbeit zur Kenntnis gekommene 
Untersuchung von P. Queney (Comptes rendus, Bd. 213 (1941), S. 588 bis 591) hingewiesen. Doch kann die a. a. 0. 
Seite 590 ohne Beweis mitgeteilte Formel für die Aufwindverteilung nicht richtig sein, da die darin vorkommenden 

1 1 
jedenfalls — au E2 


. 1 1 1 i 
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der Berechnung die atmosphärischen Störungsgleichungen für stationäre, reibungsfreie Strö- 
mung in polytroper Atmosphäre zugrunde gelegt werden. 

Zur Beschreibung der Zustandsänderungen und Bewegungsvorgänge werde in einem 
Punkt der Erdobeıfläche als Koordinatenanfang ein rechtwinkliges Koordinatensystem so ge- 
wählt, daß die z- Achse vertikal nach oben und die «-Achse in Richtung der Grundströmung 
weist,-so daß sich das Bodenhindernis parallel zur y- Achse erstreckt. Dichte, Druck, absolute 
Temperatur und Geschwindigkeitskomponenten der zweidimensionalen gestörten Strömung 
seien beziehungsweise gleich o+o, p+p, T+T und (ä+u,m), worin 0, p, T die ent- 
sprechenden bekannten, nur von 2 abhängigen Größen des Grundzustandes und (@, 0) die 
konstanten Geschwindigkeitskomponenten der Grundströmung bezeichnen. Diese ungestörte 
Luftströmung stellt wegen des räumlich und zeitlich konstanten Feldes horizontaler Ge- 
schwindigkeiten eine Atmosphäre in statischem Gleichgewicht dar. Wir setzen voraus, daß 
sie eine polytrope Atmosphäre sei, d. h. ein lineares Temperaturfeld 


N RE 


an 01 
habe, wobei T, die absolute Temperatur am Boden 2=0 und y=-7, den konstanten 


Temperaturgradienten bezeichnet. Durch die lineare Temperaturabnahme mit der Höhe ist 
— von bodennahen Schichten abgesehen — die Troposphäre gekennzeichnet, deren Höhe in 
mittleren Breiten etwa 10 km und deren vertikaler Temperaturgradient im Mittel 0,6° C|100 m 
beträgt. Um neben der „geometrischen Zustandskurve“ (1) auch die räumliche Verteilung 
von Druck und Dichte zu kennen, integriert man die statische Grundgleichung 


_2P 


ds ye ” . . . . . . . . . . . . . . . (2) 


vermittels (1) und der Gasgleichung #_RT. Man erhält dann unter Einführung der „Höhe 


u=) 


der homogenen Atmosphäre“ H=RT,/ga=p,/9 0, und der Dimensionslosen n= g/(g— Ry) für 
die Abhängigkeit der drei Zustandsgrößen von der Höhe 





1 
HER — ( n—i a) 
0e=% n H 
on Te A 
B=3.[1- n A) 
„»_al, _r 1 

= ‚(1- n 3) 


worin die zweite Beziehung in der nach 2 aufgelösten Form die, „polytrope Höhenformel“ 

darstellt. Danach bestehen für diese polytrope oder „nach der Polytropen der Klasse n auf- 

gebaute“ Atmosphäre zwischen je zwei der drei Zustandsgrößen die geometrischen Zustands- 
n—1 1 


gleichungen 
E-). 4-8”. E47” 
Po \00) ' & \% 


Dem vertikalen Temperaturgradienten y "1 z entspricht der Exponent n=g/(g— Rp). 


als 
Fleı 


Der Fall n=x=c,/c,=1,4 liefert die adiabatische Atmosphäre und die Grenzübergänge 
n>1 und n>x in bekannter Weise die isotherme bzw. die homogene (gleichförmige) Atmo- 
sphäre. 

Wegen der schnellen Druckänderungen, welche die Luftmassen in der gestörten Strö- 
mung erfahren, setzen wir für die individuellen Zustandsänderungen eines Luftteilchens vor- 
aus, daß ihm von außen Wärme weder zugeführt noch entzogen werde. Diese adiabatische 
Zustandsänderung ist nach dem ersten Hauptsatz der Thermodynamik durch die Diffe- 
rentialgleichung 


je 
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d=— 
T 
+ 1B+n 
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gegeben, wo Fr) die „individuelle“ oder „substantielle* Differentiation nach der Zeit be- 
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P=R(T+ T) liefert 


als „individuelle Zustandsgleichungen“ zwischen je zwei Zustandsgrößen die drei Formen der 
Adiabate 


zeichnet. Die Integration dieser Gleichung mit Hilfe der Gasgleichung 2 


s—1 1 

p+p _[o+o\ tt -(P+e)® o+0 [T+ 4 

Pot Ps (are) T+ En. Pot Po %+% In: T, 3 
wobei die Gleichungen x=c,/c, und AR=c,— c, benutzt wurden. Während wir in den geo- 
metrischen Zustandsgleichungen mit o,, P,, 7, die Werte der Zustandsgrößen am Boden z =0 
bezeichnet haben, bedeuten die Integrationskonstanten 9,+0,, Po +Ps; Ta + T, in den zu- 
letzt erhaltenen Gleichungen die noch von den Teilchen abhängigen Anfangswerte der Zu- 
standsgrößen. Da wir aber nicht die Lagrangesche, sondern die Eulersche Darstellung 
der Strömung verwenden, müssen diese individuellen Integrationskonstanten eliminiert werden: 
Man erhält durch individuelle logarithmische Differentiation etwa der zweiten der drei Adia- 
baten die gesuchte individuelle Zustandsgleichung in der differentiellen Form (o+o)/(p + p) 

—=x#d(e+o)/d(p-+p) oder 

d@+o) _ 1 d(p+p» 
did I u W, 
p+p a 5 nd De A 
an, R(T-+ T) das Quadrat der Schallgeschwindigkeit bezeichnet. 


wo (c+c)’=x = 
4 


Im Fall einer adiabatischen Atmosphäre mit dem Temperaturfeld T= T, — y*z stimmen 


” oT lg ER du i 
geometrischer Temperaturgradient "Z =’ _- e 7 und individueller Temperaturgradient 
d(T+T) . i \ cds ; 
ae > überein, so daß bei quasistatischem Austausch von Luftteilchen der Schwerpunkt 


der gesamten Luftmasse unverändert bleibt. Der Fall„=y* (d.h. n=x) der adiabatischen 
Atmosphäre stellt also den indifferenten Gleichgewichtszustand dar. Die polytrope Atmo- 
sphäre ist offenbar im stabilen Gleichgewicht, wenn ihr Temperaturgradient y < y* (d.h. x<x), 
und im labilen Gleichgewicht, wenn y >y* (d.h. n >x) ist. Für späteren Gebrauch schreiben 
wir den Zusammenhang zwischen y und n» noch in der folgenden Form 


n, r > 
— —1=(x ni .)= er nt BE TE 5 

2. Die Leewellengleichung. Wegen der zweidimensionalen und stationären Verhältnisse 
sind die fünf Zustandsgrößen 9+0, p+p, T+T, ü+u, w Funktionen der beiden kartesi- 
schen Lagekoordinaten x, z allein, wobei die Größen 0, », T, ü des Grundzustandes durch (3), 
die vorgeschriebenen Bodenwerte p,, T,, 0, =Pp.| RT, und die Anströmgeschwindigkeit a —const 
gegeben sind. Die Zustandsgrößen hängen durch fünf Beziehungen zusammen, nämlich die 
beiden Eulerschen Bewegungsgleichungen, die Kontinuitätsgleichung, die individuelle Zu- 
standsgleichung und die Gasgleichung. Berücksichtigt man, daß der Grundzustand von x un- 
abhängig und die Grundströmung horizontal und konstant ist, so lauten diese hydrodynamisch- 
thermodynamischen Gleichungen der Reihe nach 


E ou DET 1 op 

ru) rn oz o+0o0.x' 

x RIT dm I o(p+p») 

+), —rn - erw z. G; 
‚(@+e) 


du , Ow 
1-0 


+6+0ly, + 32 


KR; 


(ü+ u) = +n 


02 


ED ernannt) 
turn u =(c+e) rt 7 == 9, 
p+tPp. 3 ı FT n 
ver R(T+T). 


Hieraus leitet man nach der Methode der Linearisierung die atmosphärischen Störungs- 
gleichungen für die Störgrößen o, p, T, u, w ab, indem man in diesem System von Grund- 
gleichungen alle in den Störgrößen und deren Ableitungen nicht linearen Glieder vernach- 
lässigt und beachtet, daß wegen p/Jo—= RT unddp/dz= -—- 40 der Grundzustando=p=T= 
=u=mw=( eine strenge Lösung des hydrodynamisch-thermodynamischen Gleichungssystems 
darstellt. Danach erhält man die fünf linearen Störungsgleichungen 
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als Bestimmungsgleichungen für die fünf gesuchten Störgrößen o, p, T, u, w. Da zunächst 
das Aufwindfeld w(x,2) interessiert, ist das nächste Ziel die Herleitung einer Gleichung für w 
allein durch Elimination der übrigen Störgrößen. Hierzu eliminiere man einerseits p aus (7) 
und (9) und darauf o mittels (8) und andererseits p aus (6) und (7) und darauf wieder o 
mittels (8). Auf diesem von W. Tollmien angegebenen Wege erhält man zunächst zwei 
Gleichungen für die beiden Komponenten « und w der Störgeschwindigkeit, nämlich 


du _, [du er 
2 nn > r 
si + En a en a en Fr N Fe 
ww _d /du ,„ dw DEI) du , dm ; 
ee lantge) 93; -#-Dalsctz, (12). 
dx doz\0x oz 02 02) 
AD) zZ ow 


’ an : TE x gu u 
Schreibt man erstere Gleichung in der Form — @? a ?) Fr an )- gw und eliminiert ®*) 


d2 
i aM . . u mw 
man aus dieser und der letzteren Gleichung die Divergenz er: so findet man wegen 


dc/dz=—x»g(n—1)/n als Gleichung für w allein 


f ü? f (x —- 1)nü? 
—N1= - — hi 
re w "mw c? x g0w x»—n (x —n)g? 0 13 
_ = ——- — = = w=ß, . ; 3). 
dx d 2? u nc oz ü? nwo ’ (13) 
1— 0? Se e? 


Für die weitere Untersuchung beschränken wir uns auf den praktisch allein in Betracht 
kommenden. Fall, daß das Quadrat #°/c? der Machschen Zahl innerhalb der Troposphäre 
gegenüber 1 zu vernachlässigen ist. Dann geht (13) in 





x x 
vw ‚am nam (62) f 
02° e® d2 ee ER 











über. Dies ist. die Leewellengleichung, eine lineare partielle Differentialgleichung 
zweiter Ordnung mit nicht konstanten Koeffizienten. Sie lautet nach Einführung der Tempe- 


ratur T=c?’/x R und des Temperaturgefälles y=(n — 1)y/n R wegen der Beziehung (5): 


re EI, 
do RT dat ._@# r 
3. Angenäherte Berücksichtigung der polytropen Schichtung und die Länge der Leewellen. 
Die Leewellengleichung (14) geht im Falle n=1 der isothermen Atmospäre ‚in die lineare 
Differentialgleichung mit konstanten Koeffizienten über: 
mw mw zgydw (x—1)y? 


_ T rn LTE RN ne ar en, D). 
dx ' de © rt we" (15) 


Die Behandlung der allgemeinen Leewellengleichung, welche an der Grenze 2 — ni der poly- 


tropen Atmosphäre singulär ist, führt auf die konfluente hypergeometrische Differentialgleichung. 
Doch interessiert die Lösung der allgemeinen Leewellengleichung nur innerhalb der Tropo- 
sphäre, da nur in dieser die Atmosphäre polytrop, in der Stratosphäre aber isotherm ist. 


se) Für eine ausführliche Darstellung des ganzen Eliminationsprozesses vgl. man die demuächst erscheinende 
3. Auflage des Lehrbuches von H. Koschmieder: Dynamische Meteorologie, Leipzig 1943. 
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Nun weicht das Quadrat der Schallgeschwindigkeit in der Troposphäre nur um höchstens 
10 Prozent von seinem Mittelwert ab. Wenn man daher in (14) c? durch seinen konstanten Mittel- 
wert ersetzt, so wird die strenge Lösung dieser vereinfachten Differentialgleichung im Bereich 
der Troposphäre eine gute Näherung der allgemeinen Leewellengleichung darstellen. Indem 
wir in diese Näherungslösung für c? die ursprüngliche Abhängigkeit von z wieder einführen, 
haben wir den Einfluß der polytropen Schichtung mit einer Genauigkeit berücksichtigt, 
welche in den Grenzen der Linearisierung gelegen ist. Diese Berücksichtigung der poly- 
tropen Schichtung mit einem von Null verschiedenen Temperaturgradienten kommt also darauf 
hinaus, daß in der strengen Lösung des Leewellenproblems für isotherme Schichtung das 
Verhältnis x der spezifise hen Wärmen durch »/n und das Quadrat der Schallgeschwindigkeit 
durch e; {1 eh. ersetzt wird. 
N we Mi 
Es handelt sich also darum, die Leewellengleichung (14) bei konstantem c? zu lösen. 
Machen wir den Ansatz 
RE NE a a EZ nF 


wo H=&e/z9=RT/g die Höhe der homogenen Atmosphäre für die Temperatur T ist, so 
genügt w, (x,2) der Wellengleichung 














vw, =0 mi k= a4 V ” 1 vr (17) 
i r uc n EL N i 

A , a ’ } w_ »’ Wü 
Entsprechend der Voraussetzung über die Machsche Zahl wurde hierbei ee 7 ; durch 


2 -1 ersetzt. Damit ist die Bestimmung der Aufwindkomponente w auf die Te der 


viel untersuchten Wellengleichung unter den noch zu formulierenden Randbedingungen zurück- 
geführt. 
Mit dem Ausdruck (17) für die Wellenzahl k.ist auch die asymptotische Länge /=2n/k und 





n 
identisch mit der von V. Bjerknes’) angegebenen „Transversalfrequenz“ yy/’—I*, wo 


die „Orbitalfrequenz“ 2rü/A=ük der Leewellen gegeben. Die Orbitalfrequenz 2 V * _1 ist 


T=1/n RT= — fe und /*=1/x RT=1/c* den Barotropie- bzw. den Pieozotropiekoeffizi- 


enten bezeichnen. Dieselbe Transversalfrequenz erhält man auch als Kreisfrequenz der quasi- 
statischen kleinen Vertikalschwingungen eines isolierten Luftteilchens um ein mittleres 
Niveau 2 innerhalb der umgebenden, ruhenden und stabil geschichteten Luft. Denn ist £ die 
kleine Vertikalverschiebung des betrachteten Luftteilchens und T=T,— yz die gemeinsame 
Temperatur des Teilchens "und der umgebenden Luft in der Gleichgewichtshöhe 2, dann ist 
T=TN- Y@+2 Shre T—yt die Temperatur der umgebenden Luft in der Höhe z2+{[ und 
Fe ; Y 4 
=} \ -T 
adiabatisch verschobenen Luftteilchen; ferner ist die Vertikalbeschleunigung pro Masseneinheit 
ie A ; . Da d’Z T—-T 
des isolierten Luftteilchens bei Beschränkung auf quasistatische Bewegung durch de I m 
gegeben; führt man hierin für T und T’ ihre Abhängigkeit von £ ein, so ergibt sich bis auf 
kleine Größen zweiter Ordnung 


E* =T-—-y"*"c+--: die Temperatur des von der Höhe z auf diese Höhe 2+{ 


dt gW*—r 
sY—nN,._o. 
df 1 
Aus dieser einfachen Schwingungsgleichung entnimmt man die Kreisfrequenz } 9 (y*— y)/T, die 
wegen (5) mit der Orbitalfrequenz der Leewellen übereinstimmt’). Beachtet man, daß das 


#, V. Bjerknes u. Mitarbeiter: Physikalische "Hydrodynamik. Berlin 1933, S. 328 u. S. 370 ft. 


10) Diese Orbitalfrequenz Vs u* —;), T der Leewellen findet sich zum ersten Mal bei P.Raethjen: Das zwei- 
dimensionale, atmosphärische Stromfeld um ein Hindernis, Z. Flugtechn. Bd. 17 (1926), S. 185 bis 190. Vergl. auch 
Raethjens eingehende Untersuchungen: Zur Vertikalbewegung im atmosphärischen Kontinuum, Meteor. Z. Bd. 46 (1929), 
S. 292 bis 300, S. 329 bis 337, S. 383 bis 303, S. 420 bis 435, und: Stationäre Strömung unter dem Einfluß der Schwere in 
stabil geschichteten Flüssigkeiten und Gasen, insbesondere der Atmosphäre, Meteor. Z. Bd. 48 (1931), S. 288 bis 306. Auf 
denselben Ausdruck für die Orbitalfrequenz der Leewellen kommt F. Linke durch die oben wiedergegebene Überlegung 
in seiner Abhandlung: Zur Vertikalbewegung isolierter Luftimassen, Meteor. Z. Bd. 45 (1928), S. 255 bis 263. Dieselbe 
Frequenz tritt auch noch in anderem Zusammenhang auf. Z.B. hat eine in dieser Zeitschrift demnächst erscheinende 
Arbeit von H. Görtler „Über eine Schwingungserscheinung in Flüssigkeiten mit stabiler Dichteschichtung“ ergeben, 
daß bei kleinen erzwungenen Schwingungen mit einem örtlichen Regularitätsdefekt der Verschiebungsamplituden zugleich 
auch entsprechende Störungen in allen Punkten gewisser Flächen auftreten, wenn die Schwingungsfrequenz unterhalb 
einer kritischen Frequenz liegt, welche gerade gleich der Orbitalfrequenz der Leewellen ist. 
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von Th. Hesselberg eingeführte Stabilitätsmaß E= (y* — y)/ T auch durch die logarithmische 


v9 Bio . n—lz en . 
Ableitung S> | © der durch TO =(p/p,) * oder durch 9=[T, (1 — nn) definierten 


potentiellen Temperatur © gegeben ist, so erhält man nach diesen Umformungen für die Länge 
der Leewellen die folgenden Darstellungen 


2: Sn / un H_ Bun. 2rü Si / fg (18) 
A= Pr: re ung — / 9 oa sv!-r"" (y* — y) . . . . 
IV V I yz 


Insbesondere geht aus dem letzteren Ausdruck hervor, daß die Länge der Leewellen allein von der 
Temperatur, dem Temperaturgradienten und der Anströmgeschwindigkeit abhängt '°*). Setzt man 
in diese Wellenformel (18) z. B. die von J. Küttner‘"’) auf Grund zahlreicher Beobachtungen 
gemessenen Mittelwerte T=248°K, y= 0,006 C/m, @=1lömsec”' ein, welche etwa der 
Wetterlage vom 21. Mai 1935 entsprechen, so ergibt sich die Wellenlänge zu A=7700 m in 
genauer Übereinstimmung mit der Beobachtung. 


4. Formulierung der Randwertaufgabe. Zu dem System der hydrodynamisch-thermo- 
dynamischen Grundgleichungen auf S. 6 treten noch die durch das Strömungsproblem ge- 
gebenen Grenzflächenbedingungen. Da die freie Atmosphäre nur den Boden zur äußeren Be- 
grenzungsfläche hat und wir sonst keine inneren Diskontinuitätsflächen voraussetzen, kommt 
zu den Grundgleichungen nur die kinematische Grenzflächenbedingung hinzu, welche das Ver- 
schwinden der Normalkomponente der Geschwindigkeit am Boden verlangt. Ist 3={,(x) die 
Gleichung der Kontur des Bodenprofils, so lautet diese Bedingung 


(0,60) [tu] R)- 


Die Höhe £,(x) der Bodenerhebung sei für alle & gegenüber der Wellenlänge A klein von 
erster Ordnung. Da für alle & bis auf kleine Größen zweiter Ordnung u (x, £,)=u(x,0) und 
w(x,£,)=m(x,0) und wegen (16) w(x,0)= mw, (x,0) ist, so lautet die kinematische Grenz- 
flächenbedingung in linearisierter Form 





m. 0)=ü& (a) N ee | 





Hierbei wurde von erster Ordnung kleines £, (x), d. h. flach ansteigende Bodenerhebung vor- 
ausgesetzt; jedoch kann diese Randbedingung auch für beliebig steil ansteigende Plateau- 
und Berghänge beibehalten werden: denn wegen des elliptischen Typus der Differential- 
gleichung (14) wird die Steigung w (x,2)/(a+u(x,2)) der Stromlinie im Punkte (x,2) mit 
wachsender Höhe 2 sehr rasch eine für alle x kleine Größe erster Ordnung, so daß (19) schon 
für eine nahe am Boden verlaufende Stromlinie als neue Bodenkontur bis auf kleine Größen 
zweiter Ordnung eine strenge Randbedingung darstellt. Da wir nach den Darlegungen im 
vorigen Abschnitt bei der Berücksichtigung der polytropen Schichtung von der Leewellen- 
gleichung für isotherme Atmosphäre ausgingen, also eine Lösung in unendlich hoher Atmo- 
sphäre zu bestimmen haben, so ist die Bedingung (19) noch durch eine Aussage über das 
Verhalten der Größe w, im Unendlichen zu ergänzen. Nun ist für z2>+ nicht das Ver- 


schwinden der Vertikalkomponente w, sondern nur das Verschwinden der Energie s w? der 


vertikalen Störbewegung pro Volumeneinheit eine physikalisch notwendige Forderung. Da 


‘ 


in der unendlich hohen isothermen Atmosphäre 2 w’» e-:/H (ee ?2H w)’— w} ist, tritt zu (19) 


noch die Forderung 
limw,(&,2)=0 für yatz!>x..... 2.0.0.0. (0) 


als Randbedingung im Unendlichen hinzu. Damit ist die Bestimmung der Aufwindkompo- 
nente w= mw, e*/?rH auf die Lösung der folgenden Randwertaufgabe zurückgeführt: Es ist 


10«) Diese Bemerkung bezieht sich jedoch nur auf diejenigen Leewellen, welche von der störenden Bodenerhebung 
einen hinreichend großen horizontalen Abstand haben, da in der noch zu bestimmenden Lösung der Wellengleichung 


die Länge ?2ru vr 9(x*—y) nur die asymptotische Wellenlänge darstellt. Mit Formel (18) wird also nur die Wellen- 
länge der vom Bodenhindernis entfernteren Leewellen behauptet, so daß der Einwand H. Reinhardts (Beitr. z. Phys. 
d. fr. Atomsphäre, Bd. 26 (1940), S. 33) gegen die Unabhängigkeit der Wellenformel (18) vom Bodenprofil hier urizutreffend 
ist. Die bemängelte Abweichnng der am 28. V.1938 beobachteten Wellenlänge von der dieser Wetterlage entsprechenden 
theoretischen Wellenlänge rührt zu einem Teil auch daher, daß Reinhardt in Formel (18) für 7 fälschlich die Boden- 
temperatur statt die in der betrachteten Höhe herrschende Temperatur einsetzt. Im übrigen bekommt man durch eine 
geringe, noch im Bereich der Beobachtung gelegene Vergrößerung des Temperaturgradienten y eine gute Überein- 
stimmung mit der beobachteten Wellenlänge. 
10b) J. Küttner: a. a. O. insbesondere S. 102 ff. und S. 264. 











sic 
Ge 
du 
x = 


Se 
Int 


üb 
Ini 


















Z. angew. Math. Mech. . Gay = ” ä " 
Bd. Nr. 1 Febr. 1943 Lyra, Theorie der stationären Leewellenströmung in freier Atmosphäre 9 





eine im Inneren der oberen Halbebene 2>0 reguläre Lösung w, der Diffe- 
rentialgleichung Aw,+k’w,=0 zu bestimmen, welche amRand 3=0 die vor- 
geschriebenen Randwerte @£,(x) annimmt und im Unendlichen verschwindet. 


5. Greensche Funktion und Eindeutigkeitsproblem. Die erste Randwertaufgabe läßt 
sich bekanntlich auf das spezielle Problem zurückführen, eine am Rande des betrachteten 
Gebiets & verschwindende „Quellösung“ der Differentialgleichung zu bestimmen, nämlich die 
durch folgende Eigenschaften definierte Greensche Funktion G (x,2;£&,{), wobei der Punkt 

—£,2={, dem Gebiet © angehört: 


1. @ genügt in bezug auf «,z im Gebiet &, ausgenommen im Punkte x=£, 2={£, der 
Differentialgleichung AG + kK?’G=0. 


2. @ wird im singulären Punkt <=£, z={ wie 5 log V(xz — &°?+ (2 —{£)* unendlich, 
d. h. @ ist im Gebiet & von der Form U (x, 2)logy(x — &?+(z2— &)+V(x, 2), wo 
U und V iin © regulär sind und U6,9=5, ist. 


3. @ ist Null längs des Randes von ©. 


- 


Seien nun « (£,Z) und v(£,Z) irgendwelche in einem Gebiet der £Z-Ebene reguläre Funktionen. 
Integriert man die Identität 
Ö 


u(Av+k’v)— v(Au tKu)=75; 


- 


dL 


| dv du 


0v du 
a u u; —- U I 
RI- DF= 0, 0 


über ein Regularitätsbereich von « und v» und wendet auf die rechte Seite den Gaußschen 
Integralsatz an, so erhält man die Greensche Formel der Wellengleichung 
dv du 


dv ou] ,_ s 
ac | u DE v Seas, 


u dE — | DE 

wo das Kurvenintegral über den als analytisch und regulär vorausgesetzten Rand des Gebiets 
erstreckt ist. Um diese Formel auf die Greensche Funktion u=@G und eine reguläre 
Lösung v=m, der Wellengleichung anzuwenden, erstrecken wir das Flächenintegral über 
denjenigen Teil des Gebiets ©, der außerhalb eines ganz in ® enthaltenen Kreises um (x, 2) 
gelegen ist und erstrecken das Kurvenintegral im positiven Sinne über den Rand dieses Teil- 
gebietes. Geht man dann in bekannter Weise mit dem Kreisradius zur Grenze Null über, 
so findet man die Darstellung - 











\ \lutto+ Ro) vidu4 Puldsdt=| 











[ Ga... 06 
w,(2,28) = \ w.lgat JE 


. 


as, 


worin zur Ausführung des Kurvenintegrals nur die Randwerte von w, bekannt zu sein brauchen. 
Speziell für die obere Halbebene ais Gebiet geht diese Formel in 
0 
w.e,)=— \ n.(60 (5), I  — — 
J Slt=0 


—.© 


über. Handelt es sich um ein beschränktes Gebiet ©, dessen Rand aus einer geschlossenen, 
regulären und analytischen ‚Kurve besteht, so hat die Wellengleichung bekanntlich nur ent- 
weder eine endliche Anzahl linear unabhängiger und auf dem Rande verschwindender „Eigen- 
lösungen“, oder sie hat eine und nur eine in & reguläre Lösung, welche am Rand abteilungs- 
weise stetig vorgegebene Werte annimmt. In diesem Falle ist die Greensche Funktion @ 
eindeutig vorhanden, und die Lösung der ersten Randwertaufgabe durch obige Formel ge- 
geben''). Ist aber das Gebiet & unendlich, so ist jeder Wert von k, also auch (17) ein Eigen- 
wert der Wellengleichung, zu welchem es in © reguläre und am Rande verschwindende 
Eigenlösungen gibt'?). Ist insbesondere & die obere Halbebene 2>0, so ist, wenn J, die 
Zylinderfunktion erster Art und »-ter Ordnung bezeichnet, 


= 
De, J,(kr)sinv go 


vs—| 


11) L. Liechtenstein: Journ. f. Math. Bd. 142 (1913), S. 11 bis 13. 

12) A. Sommerfeld: Die Greensche Funktion der Schwingungsgleichung. Jber, dtsch. Math.-Ver. Bd. 21 (1912), 
S. 309 bis 353. Vgl. auch F. Pockels, Über die partielle Differentialgleichung Ju+k?u=t. Leipzig 1891, insbes. 
S. 206 u. S. 305 ff. 
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mit e=rcosp, 2=rsing, eine unendliche Mannigfaltigkeit reeller, überall in & regulärer 
und am Rande verschwindender Lösungen, welche jeder Lösung der ersten Randwertaufgabe 
überlagert werden kann. Im Gegensatz zur entsprechenden Randwertaufgabe der Potential- 
theorie, wo eine reguläre, am Rande und im Unendlichen verschwindende Lösung identisch 
Null sein muß, ist also die Lösung der Wellengleichung durch vorgeschriebene Randwerte 
und das Verschwinden im Unendlichen noch nicht eindeutig bestimmt. Weil aber das hydro- 
dynamische Problem der Leewellen selbstverständlich eindeutig ist, muß in der mathemati- 
schen Analyse dieses Problems eine Bedingung hinzukommen, welche aus der unendlichen 
Mannigfaltigkeit möglicher Lösungen die eine durch die Natur eindeutig realisierte Lösung 
aussondert'?). Es lassen sich indessen verschiedene hinreichende Zusatzbedingungen angeben, 
wodurch die Greensche Funktion eindeutig festgelegt wird. Im folgenden zweiten Teil 
kommen drei wesentlich verschiedene Methoden zur Beantwortung der Eindeutigkeitsfrage 
zur Darstellung. 

Versteht man unter einer Grundlösung F (.,z) der Wellengleichung eine solche Lösung, 


i EN _ ’ i i 
welche in e=z=0 wie , „logVa’ +2 unendlich wird — für die also Fix —&,2—[{) die 
beiden erstgenannten Eigenschaften der Greenschen Funktion besitzt, aber nicht notwendig 
die dritte — so erhalten wir durch die Bildung 

G(2,2;&5,0)=Fla--&2—L)—-Flza-E2+#[0) . ...: 2.2.0.2) 


offenbar eine Greensche Funktion der Wellengleichung für die obere Halbebene 2>0. Aus 
der letzteren Gleichung folgt unmittelbar 
d G(x,2;8&,0) 

0, 





so daß die Formel (21) für die Lösung der allgemeinen im vorigen Abschnitt formulierten 
Randwertaufgabe in 


= 
OF(x &,2) 


w,(2,2)=2 | w,(&, 0) 35 EN PO EEE ER 


— © 


übergeht. Mit Rücksicht auf den Ansatz w= w, e?'?"H und die Randbedingung w, (x,0) = ü £, (x) 
ist also die Lösung der Leewellengleichung für beliebige Bodenkontur £, (x) durch das Integral 











= 0) i ? “ Ö Fix E77 &E,2 . 
w(x,2)=2u er ?rH \&& _ das BE ENTE ur 
oe a 








gegeben. Damit ist das allgemeine Leewellenproblem auf die eindeutige Bestimmung einer 
Grundlösung F(x,2) der Wellengleichung zurückgeführt, 


Nun ist die allgemeinste Grundlösung der Wellengleichung bekanntlich durch 


1 > ; 
F= I N,(k r) u kr) (a „cos vo : b, sın v p) i } E r b E B (25) 

. ; Ren .= 1 { 
mit 2=rcosgp, z=rsing gegeben, worin T N, (kr) wegen 1 N,krn)=2- logr für r>0 


eine spezielle Grundlösung der Wellengleichung ist, und der additiv zugefügte Reihenausdruck 
mit den unbestimmten Koeffizienten a, und b, die allgemeinste überall reguläre Lösung der 
Wellengleichung darstellt. Das Eindeutigkeitsproblem reduziert sich also auf die Frage nach 
der eindeutigen Bestimmung der unendlich vielen unbekannten Koeffizienten a, und b,. Die 
Lösung dieses Eindeutigkeitsproblems, welche im folgenden Abschnitt gegeben wird '?*), lautet 


1 1 


art’ BR u . ei 


A2yy} 17 


13) Hierfür kommt die Sommerfeldsche „Ausstrahlungsbedingung“ nicht in Betracht, da es sich bei unserem 
Strömungsproblem darum handelt, Strömungsquellen im Unendlichen der positiven Halbachse auszuschließen. Bei 
Sommerfeld stellt die abhängige Variable der Wellengleiehung die Amplitude einer eigentlichen Schwingung dar, 
während hier das w, etwas ganz anderes ist, nämlich eine mit dem Faktor e-:/?"H versehene Geschwindigkeits- 
komponente einer Strömung. 


124) Derjenige meteorologische Leser, weleher mathematisch weniger geübt ist, kamı ohne Schaden für das 
Verständnis den Teil II überschlagen. Eine stark elementarisierte Darstellung der Integration der Leewellengleichung 
lindet man bei H. Koschmieder a. a. 0. 








D 


u = 


er 
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für alle v. Nehmen wir dieses Resultat vorweg, so lautet also die Lösung (24) der Leewellen- 
gleichung wegen (25) und (26) für beliebige Bodenkontur 


BRETT a EEE => 1 > cos2r+1)o] ,. 
w=2u et? mu zug) 7Nkr)+— \ Isr,. (kr) +8 dE. 


J 2v—+1 
= v—( 
Dieser Integralausdruck wird in Teil III für spezielle Bodenerhebungen auf eine leicht zu 
tabellierende Form gebracht. 


II. Die allgemeine Lösung der Randwertaufgabe. 


6. Bestimmung der Grundlösung nach der Methode der unendlich vielen Variablen. Be- 
trachtet man die stationäre Strömung (&@ + u, mw) übeı eine Bodenerhebung vom Standpunkt 
eines relativ zur Grundströmung (@,0) ruhenden Koordinatensystems, so handelt es sich um 
diejenige instationäre Strömung, die durch die Bewegung einer mit konstanter Geschwindig- 
keit — @.in ursprünglich ruhiger Atmosphäre fortschreitenden Bodenerhebung hervorgerufen 
wird. Eine physikalische Überlegung lehrt, daß dann die einzelnen Luftteilchen Bewegungen 
um eine mittlere Ruhelage ausführen, während hinter der bewegten Bodenstörung ein Zug 
interner Wellen als Wellengruppe mit der Geschwindigkeit der Bodenstörung fortschreitet. 
Dem entspricht die Tatsache, daß in der vom Standpunkt der Bodenerhebung aus stationären 
Strömung die Störung der ursprünglich geraden horizontalen Stromlinien und damit die Auf: 
windverteilung w (x,2) längs jeder Niveaulinie 2=const. auf der Leeseite in hinreichender 
Entfernung von der Bodenerhebung aperiodisch abklingt. Diese Eigenschaft, welche sich 
offenbar auch auf w, (x,2) und damit auf die Grundlösung F(x,z) überträgt, enthält ins- 
besondere, daß .F(x,0) für alle negativen x hinreichend großen Betrages monoton ist. 

Diese letztere Eigenschaft, daß die Grundlösung für z=0 eine für alle hinreichend 
großen Werte von — x monotone Funktion ist, nehmen wir als diejenige Zusatzbedingung, 
welche aus der unendlichen Mannigfaltigkeit möglicher Grundlösungen 


1 © 
'=TN(k r)+ P; J,(kr)(a,cosv»p+b,sinv®). . .». ». 2...) 
‚=® 
die eine durch die Natur eindeutig bestimmte Grundlösung aussondert. Wir zeigen im folgen- 
den, daß durch diese Monotoniebedingung in der Tat die abzählbar unendlich vielen 
Koeffizienten a, und b, eindeutig bestimmt sind. 

Zunächst sind — ganz unabhängig von der Monotoniebedingung — alle b,=0. Um das 
zu beweisen, setzen wir in Formel (23) (S. 10) für F nur den allgemeinen regulären Bestandteil 
2 J, (kr) (a, cos v p-+b,„sin»g) ein. Dann muß die resultierende Lösung w, auf dem Rande 
r=( y 
z2=0( identisch in & verschwinden, da ja die Randbedingung bereits durch den singulären 


Bestandteil TN,(k r) erfüllt wird. Also ist für alle x das Integral 


© ß) n x J,(kr)(a,cosrp + b,sinvop) 
\ w, (£,0) v„—( 2 


02 


worin r und @ durch 2 — £=rcosp, z2=rsing gegeben sind. Da aber diese Aussage für 
ganz beliebige Randfunktionen w, (£,0) gilt, so folgt nach einem bekannten Schluß, daß dann 
notwendig 


an 


d I J,(kr)(a,cosvp +b,sinvgp) 


für alle x ist. Führt man die Differentiation nach z aus, so findet man für die Ableitung an 
der Stelle z=0, x — £E>0: 


n 


2 vb,J,(kr)—=0 


v‚—| 


für aller >0. Wären nicht alle d,=0, so sei b„ mit mZ1 der erste nicht verschwin- 
dende Koeffizient. Dividiert man die Reihe gliedweise durch r”" und vollzieht den Grenz- 


» \mı 
übergang r>0, so ergibt sich wegen J„(kr) = (3) Imt, daß b)„=0 ist im Widerspruch zur 


Voraussetzung. Daher sind alle b,=0 v=1,2,3, ...), W. z. b. w. Also muß jede Grund- 
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lösung F(x,2) der Wellengleichung, wenn F(x — &,2—{)— F(x —£,z+{£) eine Greensche 
Funktion der oberen Halbebene sein soll, eine zur x&-Achse symmetrische, d. h. in z gerade 
Funktion sein. 

Somit ist die Grundlösung von der Form 


1 PR 
F= 7 N,(kr)+ Y a,J,(kr)cosvo. 
Ta 
‚— U 


worin die noch unendlich vielen unbekannten Koeffizienten a,, a,,a,,... so zu bestimmen 
sind, daß die Monotoniebedingung erfüllt ist. Diese besagt, daß F für =n, d.h. 


gr pe 7 
a T N,(kr)+ >) (— 1)” a,J,(kr) 
es | 
eine für alle hinreichend großen r monotone Funktion ist. Um bei der Bestimmung der 
Koeffizienten a, schwierigere Konvergenzuntersuchungen zu vermeiden, gehen wir ‚zunächst 
von derjenigen Teilmenge aller Grundlösungen F bzw. aller Funktionen f(r) aus, für welche 
bei fest vorgegebenem n>0 alle a, mit v>2n-+1 verschwinden. Wir betrachten also zu- 
nächst den endlichen Summenausdruck 
1 2n+1 
„Ks 1 N,(kr)—+ > (1) a,J,(kr). 
2, 
Um für hinreichend große r ein möglichst aperiodisches Verhalten der Funktion f„(r) zu be- 
kommen, werden wir die Koeffizienten a, so bestimmen, daß die in der asymptotischen 
Entwicklung dieser Funktion enthaltenen Wellen möglichst weitgehend gelöscht werden. 
Setzt man in f„(r) für die 2#»+2 Zylinderfunktionen ihre bekannten asymptotischen Ent- 


u \ 1 ; PR 
wieklungen **) ein und ordnet nach Potenzen von zZ, 50 erhält man die für alle hinreichend 


großen r definierte asymptotische Entwicklung 


7 (— 1)4=! A 2u+l 2 5:3 
V ;krnn= > NBErw Au sin(k r— — a) Bucos|k r— er a) 


“—=U 
mit 
ee 1 ui 
. 1 -1°7-.3°...2u—1) 
n 
+ I (-1P 2» +1? — 1] 2° +1? — 32]... [42 +? —- Qu—1%]a,,;;; 
vv 
n 
Bu = 2 (- 1 42»? — PP] (27? — 3]... @r? —- Qu —1P]a,, 
r—( 
1 n " 
für a>0 und ,=-7+ I (1a... B= 2 (-1”a,,. 
0 v 
Das u-te Glied dieser Entwicklung stellt eine Welle dar, deren Amplituden mit r > 
9 . \ ie 
wie | abklingen. Bei den 2» +2 uns zur Verfügung stehenden Koeffizienten a, erreichen 
wir offenbar die stärkst mögliche Löschung der in dieser Entwicklung enthaltenen Wellen 
dadurch, daß wir die ersten 2» +2 Ausdrücke A,, A,,...,;, Ay, Ba; Bi, --, B„ zum Ver- 
schwinden bringen. Wir erhalten somit ein System von 2n +2 Bestimmungsgleichungen für 
die 2» +2 unbekannten Koeffizenten «a,, @,,..., Aans+.. Dieses System zerfällt offenbar in 
ein System von n+1 linearen Gleichungen für die Koeffizienten mit ungeradem Index und 
in ein ebensolches System für die Koeffizienten mit geradem Index: 
r & 
pi u 1’ a,»;, ee r Pi} ( 1)’ d,y 
v 


I) 
n 


n 
eV +1?-1]a,.=— ZN’ BeN’—1]a,,=0. 


12) Vgl.z.B. Jahnke-Eımde: Funktionentafeln, 2. Aufl., Leipzig 1933, S. 203, oder R. Weyrich:.Die Zylinder- 
Zunktionen und ihre Anwendungen, Leipzig 1937, S. 47. 
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Multipliziert man in den einzelnen Gleichungen die eckigen Klammern miteinander aus und 
formt in naheliegender Weise der Reihe nach die zweite, dritte, ... Gleichung unter Be- 
nutzung der vorangehenden Gleichungen um, so ergibt sich — wie allgemein durch voll- 
ständige Induktion erwiesen wird —, daß die gesuchten Koeffizienten auch den folgenden 
beiden Systemen linearer Gleichungen genügen, welche mit den obigen äquivalent sind: 

n n 

2 (- 1)’ a,,;; -ı 2(- 1)’ a,, 0, 

v v0 


n n 
Ser tltan., 0 I (-1P Era, =0. 


ZN’ +1ra,,,,—0 I (-1 @ry*ra,=0. 

0 0 
Die Determinante jedes Gleichungssystems berechnet sich als Vandermondesche Deter- 
minante'°) zu 

z2arnng2141...@n)! bzw +z2armtmar.113151...2n — 1)! 


Beide Determinanten sind für jedes » von Null verschieden, so daß die 2» +2 Koeffizienten 
Ay, dy..-, Agny, für jedes n eindentig bestimmt sind. Weil das zweite Gleichungssystem 
überdies homogen ist, so sind die Koeffizienten mit geradem Index sämtlich gleich Null: 


I era ne 


’ ’ 


Auch die Lösung des ersten Systems läßt sich in geschlossener Form angeben; man findet 


ee -;_ ae en 4 
Tg, ie m ) n—v, (#). 
Hierbei wurde die Abhängigkeit des Koeffizienten von n durch einen oberen Index zum Aus- 
druck gebracht. Somit ist f„(r) und damit „die »-te Näherung“ 


2v +1 


1 m 
F=gNtkn+Nay, IvulkrlosytDg 2.0... 0) 


der Grundlösung eindeutig bestimmt. 
Nun vollziehen wir den Grenzübergang lim F, für n>x und weisen nach, daß dieser 
Limes die gesuchte, für = bezüglich r monotone Grundlösung darstellt. Zunächst ergibt 


=: ? -_ 
sich wegen der Identität er u. Me ” fer - = I ) ( ")* . ”) aus der „Wallischen 
n n—v n+1 r v 


—1 2 . 
Produktformel* (— r( K )= TE für jedes v>0, daß 
SE. 
. (n) a 2 se > > 
Be LE a2r+1 


(6). 


Der Ausdruck (5) kann wegen aS),,—0 für »>n auch als unendliche Reihe geschrieben 


werden. Diese Reihe hat wegen der Beschränktheit der a, + und der Abschätzung 


(2) 
,kn=\, 1+0(-)| RETTEN 


NRW > krir +! Ä ” a0 
für beschränkte r< R die Majorante K (5 ) (2»+1)!, ist also für jedes n in jedem 
v—=U 
beschränkten Bereich der r g-Ebene gleichmäßig und absolut konvergent. Daher darf in (5) 
der Grenzübergang für n >» gliedweise ausgeführt werden. Also sind die gesucliten Koeffi- 


zienten a,,,, durch den Limes (6) gegeben, und damit die Grundlösung 





en 2 > „cos@r+1)o 
F=7Ndikr)+z J,y+ (kr) 2> +1 


1 ] 











15) Vgl. z. B. G. Kowalewski: Einführung in die Determiuantentheorie, 2. Aufl., Berlin 1925, S. 38, 
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eindeutig bestimmt '®). 
Es bleibt nun noch zu zeigen, daß diese so bestimmte Grundlösung (8) auch wirklich 
die Monotoniebedingung erfüllt. Dies folgt am einfachsten aus der Integraldarstellung 


1 “ ; 
= 2 N, (kr) + 2 \ en reg 4: 
$ 


welche unmittelbar aus der wegen (7) für jedes r bezüglich t gleichmäßig konvergenten Ent- 
wicklung’) des Integranden 


x 


sin (krsin)=2 N J,,,, (kr)sin(2v +1) 
v=0 
durch gliedweise Integration hervorgeht. Aus (9) folgt unter Anwendung der bekannten 


Integraldarstellung **) 


ERSE: 5 
Ben 5 N — krö&int 0. 
N,(kr)=—; \e di+5z, 


ar 


BEER: 4 5:20 A, 


e 
na 


0 
daß insbesondere F für g=n, d.h. 


1 BERTR 
fa=—z, | tens. Be a I 6 
v 


ist. Die Monotoniebedingung, wonach f(r) für alle hinreichend großen r monoton sein soll, 
ist in der Tat erfüllt, da 


fr.) —- f(r,) = = \e*" mn. erkennen 
0 


für O<r,<r, gewiß positiv ist. 

Schließlich merken wir noch an, daß die eindeutig bestimmte Grundlösung (8) nicht nur 
für 2=0 und x<0, sondern auch ve jedes z und hinreichend große Werte von — x monoton 
1 
a2». — -1 
auf der Luvseite («e< 0) und zugleich eine optimale Verstärkung der Wellen auf der Leeseite 
(x >0). Dies ersieht man leicht aus den folgenden einfachen Umformungen der-Grundlösung (9) 
mittels (10): 


ist. Die Koeffizienten a,,,, = bewirken eine stärkst mögliche Löschung der Wellen 


an. © 
F= = \ sin (kr sin Hdt-z- \ ekröint gy 
0 


0 


; 5 1 


2n 





% 
ae \ sin (kr sin Hdt+z- | erkröintgg + I N.(kr) 
ö N 


16) In gewissem Sinne sind die eindeutig bestimmten Koeffizienten (6) Lösungen des Systems von unendlich 
vielen linearen Gleichungen mit unendlich vielen Unbekannten: 


n 1@=0) 
S’(- 17 (2r +1)2P ca +1 - g$ 
v0 =,,2. 3...) 


Diese Gleichungen sindzwar — von der ersten Gleichung abgesehen — für die Lösung (6) alle divergent. Doch sind sie 
„4-summierbar in dem Sinne, daß die Summe der für |x| < 1 konvergenten Reihe 


= © 
SD eWwer +2,40 = a DB -rar+nep-1ar 
v0 v0 


für <>1-—0 gegen t/, oder 0 strebt, je nachdem »=0 oder y>Pist. In diesem Sinne ist die Lösung (6) allerdings 
nieht eindeutig bestimmt, da offenbar jedes Polynom der Form 


Co2r +1! +C,@r +1 +... +Cm@r+12 +1 o9=0,1,2,...) 


eine Lösung des homogenen Systems darstellt. Vermutlich gibt es aber nur eine Lösung, nämlich (6), so daß die 
p-te Reihe (p=0,1,2,...) genau „C2p-summierbar‘ zur angegebenen Summe ist. Die letztere Aussage ist damit 
äquivalent, daß die 2p-fach iterierten arithmetischen Mittel der Teilsummen der p-ten Reihe gegen die angegebene 
Summe konvergieren. 

17) Weyrich:a. a. O. S. 9. 

18) Weyrich:a.a. 0. S. 22 oder 8. 11. 
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Wegen 2 —p>0 für &>-— bei festem z folgert man leicht aus der ersten Darstellung, 
Br . 008 k a 
daß F(x,z) für hinreichend große Werte von — x monoton wie = Er abklingt; wegen 9>0 


sin (k x — n/4) 


für 2&>+% bei festem 2 verhält sich aber F(x,z) wie ‚ also nicht monoton. 


V2rkx 
Insbesondere folgt aus (12) für z=0 die Darstellung 
er 
| ring für z<0, 
an, 
0 
FF, = 
5, 1° ( 5 Noel x) Br: 2 >89, 
0 





Als n-te Näherung für die Grundlösung ist sowohl die »-te Teilsumme der Reihe (8) als auch 
die Summe (5) brauchbar. Erstere approximiert offenbar besser die Umgebung von r=0, 
letztere diejenige vonr=». Z. B. liefert (5) für n=0 die Grundlösung 


1 1 
7 N, (kr) + 7 J, (kr) cosp "). 


Bevor wir zur numerischen Auswertung des gewonnenen Resultats übergehen (Teil III), sollen 
in den beiden folgenden Abschnitten noch zwei weitere Methoden zur eindeutigen Bestimmung 
der Greenschen Funktion erörtert werden. 


7. Bestimmung der Grundlösung nach dem Rayleighschen Limitierungsverfahren.. Die 
vorangehende Methode hat den Mangel, daß sich die Leeseitigkeit der stationären Wellen 
nicht als Folge der Theorie ergibt, sondern vielmehr als physikalisch begründete Voraus- 
setzung in der Monotoniebedingung schon enthalten ist, die zur Behebung der sich einstellen- 
den Mehrdeutigkeit des Problems hinzugekommen war. Daher soll im folgenden noch eine 
andere Methode zur Darstellung kommen, welche auch die Leeseitigkeit der stationären Wellen 
aus den Grundgleichungen der Mechanik herzuleiten gestattet. Diese Methode hat zum ersten- 
mal Rayleigh 1883 a. a. O. zur Berechnung der (durch stationären Oberflächendruck hervor- 
gerufenen) Leewellen auf der freien Oberfläche einer unendlich tiefen inkompressiblen Strömung 
angewandt. Sie beruht auf einem Kunstgriff, welcher darin besteht, daß die Unbestimmtheit des 
Problems zunächst durch Einführung gewisser kleiner Reibungskräfte behoben wird. Von diesen 
Reibungskräften wird angenommen, daß sie der Störgeschwindigkeit (u, w) proportional sind 
und also der Störbewegung einen Widerstand (— „u, — um) mit positivem u entgegensetzen. 
Indem man in der ‘hierdurch eindeutig bestimmten Lösung des Problems mit dem positiven 
#>0 geht, erhält man die Lösung des reibungsfreien Problems. Diese Rayleighsche 
Methode läßt sich, wie H. Stümke in seinem im einleitenden Abschnitt genannten Bericht 
für die isotherme Atmosphäre gezeigt hat, auf den vorliegenden Fall eines kompressiblen 
und polytrop geschichteten Mediums in der folgenden Weise übertragen. 

Durch Einführung der genannten Reibungskräfte gehen die beiden Eulerschen Be- 
wegungsgleichungen (6) und (7) (S. 6) in 

.du_ dp 

07 7 een TA Laie 

m - u. ; 

u u TE. de u. 
über, so daß diese zusammen mit der Kontiuuitätsgleichung (8), der Zustandsgleichung (9) und 
der Gasgleichung (10) (S. 6) das System der atmosphärischen Störungsgleichungen bilden. 
Aus diesen 5 Bestimmungsgleichungen für o, p, T, w und w folgt durch Elimination der vier 
Störgrößen o, p, T, u und Berücksichtigung des Falles «’/® <1 eine Differentialgleichung 
dritter Ordnung für w allein: 


x x * 
o |öw, dm _nIom (2-2) rule en_nTdmi_ 
17 


2 rn. st+32=-52|=0 .. (13). 
ox|dx? 02? ce 02 u? c? 0x" 02 ce 02 (13) 

1%, Diese Funktion wurde in der einleitend genannten Arbeit des Verf. als Grundlösung gegeben (das entgegen- 
gesetzte Vorzeichen rührt daher, daß hier die Greensche Funktion mit anderem Vorzeichen definiert ist). Dem Verf. 


war entgangen, daß die Reihen (14) und (17) a. a. O. S. 201 identisch sind, wie eine Umformung der Glieder und nach- 
folgende Umordnung ergibt. Daher ist insbesondere die a. a. ©. S. 2,2 (Formel 19) angegebene Lösung für w nur das 
erste Glied einer unendlichen Folge partikulärer Lösungen, welche in jedem beschränkten Teilbereich der oberen 
Halbebene gleichmäßig gegen die durch die Natur eindeutig bestimmte Lösung (4) auf S. 23 der vorliegenden Arbeit 
kouvergiert. Die früher mitgeteilte Lösung stellt also lediglich eine erste Näherung für die eindeutig bestimmte 
Lösung dar, von der sie sich aber wegen der starken Konvergenz der Folge (5) in einiger Entfernung vom Störungs- 
zentrum praktisch kaum unterscheidet. 
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Hält man wieder © konstant und setzt ww, '?"#, wo H=&/xg9g=RT]/g die der Tempe- 
ratur T entsprechende Höhe der homogenen Atmosphäre bezeichnet, so genügt w, der linearen 
Differentialgleichung 3. Ordnung 


Ö 
fu (= m. = ar. w,+k? w)+[4w,—r A 
mit % 
nn /#-1 und 1 15 
57] - sion ns MR SR FE Sa SE Fer (15). 


Die hier gestellte Randwertaufgabe, eine im Inneren der oberen Halbebene reguläre 
Lösung w, der Differentialgleichung (14) zu bestimmen, welche am Rande 2=0 die vor- 
geschriebenen Randwerte &L, (x) annimmt, läßt sich auch auf die Bestimmung einer wie in 
Abschnitt 5 zu definierenden Greenschen Funktion zurückführen. Es wird sich im folgenden 
zeigen, daß eine Greensche Funktion bzw. Grundlösung des vorliegenden Problems bei hin- 
reichend kleinem Reibungskoeffizienten « >0 eindeutig vorhanden ist. Geht man mit dieser 
bestimmten Grundlösung zum rechtsseitigen Limes für «0 über, so ergibt sich in der Tat 
gerade diejenige Lösung des reibungsfreien Problems, die. wir nach der Methode des vorigen 
Abschnitts erhalten haben.‘ Das Ziel der folgenden Untersuchung ist also die eindeutige Be- 
stimmung der Grundlösung und damit der Greenschen Funktion für die Differentialgleichung 








Ö 
A me -(Aw, +k:w, I (16) 








auf dem Wege über die eindeutig vorhandene Grundlösung der Differentialgleichung (14). 

Bevor wir an diese Aufgabe herantreten, verschaffen wir uns erst die der Formel (21) 
(S. 9) analoge für die Randwertaufgabe der Differentialgleichung (16) gültige Formel. Seien 
wieder u (£,£) und v(£,£) irgendwelche in einem Gebiet der £{-Ebene reguläre Funktionen. 
Integriert man die Identität 


oP 20 
pr 2 I; 
it uns (do+ko)te; ut get ge 
mi 
2 
Peoge-Jegetulde +k’v) 
und 


u u Dvodu 

PJEac Lo: 
über ein Regularitätsbereich der Funktionen « und » und wendet auf der rechten Seite den 
Gaußschen Integralsatz an, dann erhält man als Greensche Formel der Differential- 


gleichung (16) 
\ \ u 32 (4 e+k’v)+v 3 (du+k?u) 


wo das Kurvenintegral wieder über den als analytisch und regulär vorausgesetzten Rand des 
betreffenden Gebietes erstreckt ist. Wendet man diese Formel auf die noch zu bestimmende 
Greensche Funktion u=G (x, 2;£&,{) und eine reguläre Lösung v=mw,(x,z) der Differential- 
gleichung (16) an, indem man den Integrationsbereich in der in Abschnitt 5 angegebenen 
Weise wählt, und mit dem Radius des Kreises um £=x, © =z gegen Null geht, so findet man 
\ 0%G vr... 
w,(x,2) = \mı de d°- Era | 
Diese Formel, welche speziell für das Gebiet der oberen Halbebene in 








dsdi=\(Pdr— ads), 





oo : »G 
w, (x, 2) = \ w,(£,0) ( ER e)_ „A 


Yn 


(17) 
übergeht, liefert die Lösung der vorgelegten allgemeinen Randwertaufgabe für beliebige ab- 
teilungsweise stetig vorgeschriebene Randwerte von w,, sobald die Greensche Funktion @ 
eindeutig bestimmt ist. Die Bestimmung der Greenschen Funktion führen wir auch wieder 
auf die Bestimmung der Grundlösung F(x,z) zurück, aus der wir durch die Bildung 
G(x,2;&,0)=Flx —£E,2—0)— F(x—£E2+{L) die gesuchte Greensche Funktion für die 
obere Halbebene erhalten. Aus dieser Bildung folgt unmittelbar 


@lm,2;$ eg 0 F(x—£,2) 


dE8L 


= 
ze Oxd2 
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so daß mit Rücksicht auf den Ansatz w=w, e®/?*# und die Randbedingung w (x, 0) = ü £ (x) 
die Lösung der allgemeinen Leewellengleichung (13) für beliebige Bodenkontur £,(x) durch- 
die Integralformel 








+» 
w(e,2)=2ue?"H ET Ta F@—S2) 4: 


Pr d£ IE 1 











gegeben ist. 


Die eigentliche mathematische Aufgabe besteht somit in der eindeutigen Bestimmung 
der Grundlösung F(x,z) der Differentialgleichung (16). 

Da es sich hier nicht mehr um die Wellengleichung, sondern um die allgemeinere 
Differentialgleichung (16) handelt, gehen wir von einer von N. Zeilon?®) angegebenen Defi- 
nition und Integraldarstellung einer Grundlösung für allgemeine lineare partielle Differential- 
gleichungen mit konstanten Koeffizienten aus. Bekanntlich ist für eine Grundlösung F(x, 2) 


der Wellengleichung (55 .)u = =Adu+k’u=0 charakteristisch, daß für eine in einem Ge- 


’dz 
biet & willkürlich gegebene Funktion @ (x,z) das über © erstreckte Integral 


us)=||Fa@-&:2-VPo@Ydedt ....::2:...09 


0 ® 
.. ‚use 2) dar- 
stellt. Benutzt man diese Eigenschaft als Definition einer Grundlösung für eine ganz all- 
gemeine partielle lineare Differentialgleichung mit konstantem Koeffizienten 


0 0 
f (5 3) =0, 
und setzt von der Funktion p(x,2) lediglich die Existenz ihrer Fourierschen Integral- 
darstellung 


in & eine Lösung der zugehörigen inhomogenen Differentialgleichung f (= 


ee \ \\ elle -Htr@-Dln(e t)d&dtdady 
—.® —o & 


0 0 
y. ‚z)u=e len) wegen 


N) Ö BR - ila a BR: . . 
17 rn er 


formal durch den Realteil des Integrals 


voraus, so wird die Differentialgleichung r(; 


o 


1 P> eile(«—-H+rY@—D)] 
ie \ \ I) ren  P@natatace, 


—ıo —o 





gelöst. Hieraus entnimmt man nach (19) für die Grundlösung den Zeilonschen Ansatz 


x 


1 ” eilax +2) 
Fa)=g;R | | uam‘ ady a ar a BE RN (20). 


Für die Differentialgleichung (14) ist 
fua,iy=—ia(l@ +) (@--y’+P), 
also die zugehörige Grundlösung 


Ks} n 


; i(ax+YB)dad 
F(a,2)= Reg \ \ e u. y 
"ale +r-)—Zila+y+P) 





(21). 


Damit ist die Bestimmung der Grundlösung auf die Ermittlung des rechtsseitigen Grenzwertes 
dieses uneigentlichen Doppelintegrals (21) für «>0 zurückgeführt. Ist „ >0, so konvergiert, 
wie 8! ch zeigen läßt, dieses uneigentliche Integral für jedes x und z außer für =2=0. Für 





20) N. Zeilon: Das Fundamentalintegral der allgemeinen partiellen linearen Differentialgleichung mit konstanten 
Koeffizienten, Ark. för Matem. Bd. 9(1911), Nr. 38. Vgl. auch die Darstellung in Frank-v. Mises: Differentialgleichungen 
der Physik, Teil I, 2. Aufl., Braunschweig 1930, S. 862 ft. 
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„0 dagegen hat das Integral keinen Sinn, da der Integrand wegen f,(ia,iy)=—ia(a?+y’—k?) 
in allen Punkten des Kreises a@ +y?=k” und der Geraden a=0 der ay-Ebene singulär ist. 
Wohl aber existiert der rechtsseitige Limes des Integrals für «0, durch welchen der Wert 
des Integrals definiert wird. Die Rayleighsche Methode kommt also auf ein Verfahren 
zur Limitierung divergenter Integrale hinaus. 

Führen wir die Integration zunächst bei festgehaltenem y über a von —co bis + aus, 
so kann das Ergebnis dieser Integration mit Hilfe des Residuensatzes folgendermaßen ge- 
wonnen werden. Da in der Ebene der komplexen a= Rei® der Integrand wegen 


\destrel_g-8- eine 


im Falle &<>0 auf dem unendlich großen Halbkreis der oberen Halbebene (sin © >0), im 
Falle x < 0 aber auf dem unendlich großen Halbkreis der unteren Halbebene (sin 9<0) ver- 
schwindet, kann das Integral 


1 ” gi(ax + 12) 
4m’ | Arüninte 
auch durch komplexe Integration über denjenigen geschlossenen Weg gewonnen werden, 
welcher aus der reellen Achse und dem oberen oder unteren unendlich großen Halbkreis der 
komplexen a-Ebene besteht, je nachdem x positiv oder negativ ist. Danach ist der Wert 
des Integrals (22) gleich der Summe der Residuen des Integranden in denjenigen Nullstellen 
von fu(ia,iy), welche im Inneren der oberen oder unteren Halbebene gelegen sind, je nach- 
dem x positiv oder negativ ist. Da es nur auf den rechtsseitigen Limes des Integrals (22) 
für «>0 ankommt, und die Wurzeln der algebraischen Gleichung fu (ia,iy)=0 stetig von u 
abhängen, genügt es, die Residuen des Integranden für a=0 zu kennen, wobei aber das 
Residuum einer reellen Nullstelle von f,(ia,iyJ)=—ial(a?+y’—k?) zur oberen 
oder unteren Halbebene zu rechnen ist, je nachdem, ob diese Nullstelle für 
kleine positive « in die obere oder untere Halbebene rückt. Um aber das Ver- 
halten der Nullstellen von f,(ia,iy) festzustellen, wenn u stetig von Null auf einen kleinen 
positiven Wert anwächst, genügt es, die Wurzeln von fu(ia,iy)=0 nur bis auf Glieder, 
welche in « von 2. Ordnung sind, zu bestinnmen. Die elementare Rechnung ergibt 





+ 2 
IT ee Zen 20 22217727 
a 
a, (u) = VEZr+G Bro), 
PEREN iu k®+l: 
a,(u)=— Bar oe. 


Hieraus folgt: das Re- 
siduum in a, (0)=0 gehört 
zur unteren oder oberen 
Halbebene, je nachdem 
Iy|<k oder |y|>k ist; 
die Residuen in den (für 
|yI< k reellen) Nullstellen 





An 
—- 





a, (= yk’—y' 
und 

a, ()=— yk’— y? 
gehören beide zur oberen 
Halbebene. Die Wahl des 
Integrationsweges erfolgt 











&%o Rn 
x<o0 
a, M>K 





also — den vier Fällen 
z<=0,|y)|ZSk 
entsprechend — in der 
in Bild 2 angedeuteten 
wu Bild 2. 
Die Rayleighsche Methode, welche in der Definition 
1 K- en ei(«x +72) 
a er. ——d = li \ Pa 7° 2 . . . . % 
im) (a, ) aa, ® Br 
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(22) 
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des für «=0 an sich ieh Integrals (22) besteht, kann auch so aufgefaßt werden: Dem 
Integranden von (23) kann zunächst dadurch ein Sinn beigelegt werden, daß die Gerade 
y=const. der reellen ay-Ebene zur komplexen a-Ebene erweitert wird und die reellen 
Pole des Integranden auf dem Umweg über diese komplexe Ebene umgangen werden. Dabei 
stellt sich insofern eine Unbestimmtheit der Aufgabe ein, als jeder der reellen Pole entweder 
über die obere oder die untere Halbebene umgangen werden kann, so daß sich mehrere Mög- 
lichkeiten für die Wahl des Integrationsweges bieten. Diese Mehrdeutigkeit des Problems 
wird nun durch die Einführung der Rayleighschen Reibungskräfte behoben, indem der 
Integrationsweg in der oben ausgeführten Weise eindeutig festgelegt wird. 

Nun zeigt sich, daß das Integral der in dieser Weise bestimmten Funktion (23) über y 
von — oo bis +0 nicht existiert. Diese Schwierigkeit ist aber leicht beseitigt, wenn man 
von (23) das in derselben Weise ausgewertete Integral 

nn ” ei’z 
”)| f, (ia, „4 


subtrahiert, wodurch offenbar die eine Singularität a—=0 des Integranden behoben wird. 
Durch diese Änderung geht das Grundintegral (20) für a=0 in 


1 e © gilax +72) __ eirz 
\ | rn dalr er waere 


über, was auch wieder ein Grundintegral ist, da ei’*/f,(ia,iy) von x unabhängig, also 


2.2 ei1?2/f,(ia,iy)=0 ist. Im übrigen erhielte man auch ohne diese Abänderung der 
\dx’d2 g 


Grundlösung bei formaler Rechnung dasselbe Resultat für w,, da man nach (18) nur die Ab- 
%”F 

leitung 2 ee e) benötigt, in welcher das additiv hinzugetretene Glied wieder fortfallen 

würde. 


Um endlich das Grundiniegral auszuwerten, bestimmen wir die Residuen des Integranden 
in bezug auf a. Diese lassen sich leicht aus der Umformung 
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ablesen, wobei wir die erste oder zweite Darstellung verwenden, je nachdem |y|< k oder 
Iy|>k ist. 

Unter Beachtung des oben angegebenen Integrationsweges (wobei nunmehr a=0 keine 
singuläre Stelle ist, (@°**—1)/a ist-für a=0 natürlich durch ix zu ersetzen) findet man den 
vier Fällen <=0, |y|=k entsprechend 
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Hieraus folgt endlich durch Integration nach y von —» bis + und Beschränkung auf 
den Realteil 





5, | E12 (er VRR _1)dy (<<0), 
F (x,2) =! 





K 
1 ? k’— y? 
er AA a _ 8 
1 ee rent | le PR —1)dy (x >0). 





Damit ist die gesuchte Grundlösung der Differentialgleichung (16) des reibungsfreien Problems 
eindeutig bestimmt. Nach Formel (18) benötigen wir für die Lösung nur die Ableitung 
2» 
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Zum Schluß dieses Abschnittes soll für diesen Ausdruck noch die einheitliche für alle «=F0 


gültige Stümkesche Integraldarstellung hergeleitet werden. 
Sei <<0, oder, wenn man durch <=rcosp, 2=rsing Polarkoordinaten einführt, 


F<p< 7 und r >0. Macht mıan ferner die Substitution y=kCofr, so wird 
| en ;e’/r-kdy= \ cos (k r sin @ Eoj r) etreosy Sintgr 


vr — 





v 


oo 
—Re| erreine+indr. 


Nun kann nach dem Cauchy schen Integralsatz 
in der (r+io)-Ebene (vgl. Bild 3) der Inte- 
grationsweg o=9, r>0 durch den Weg r=0, 
p<sosa und s=r, r>0 ersetzt werden, da 


' . . n In RN 
coso&Sinr im Streifen Z<0o< 5 negativ ist 





und daher 
jet Since +io)| ag eeascSint,g 
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strebt für >+©. Also wird wegen ®Sinio 48 
=isino, Sin(+in)=— Eintr: En. 


x nt 


| cosy2 pe RAR | ekrSinio;dc-+ \ ekrSin(e+im) gg 
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vr- 


y 
n 


=— | sin(krsino)do+ | er*r&intgr 
" 0 


=-22|7 N, (kn)+an e n(krsint) dt|, 
y 
letzteres nach der Poissonschen Integraldarstellung der Zylinderfunktion N,. 


Sei jetzt e>0, d.h. — bed <op< 5; substituiert man y=kCojr falls y>k und y=kcost 
falls O<y<k, so wird 





88 „ zVagy_g| 0sr2 sinzyk’—y? ,, 
vr ° : dy-2 ve dy 
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[.) ; 0 
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Wegen (25) haben wir also die für +0 gültige Darstellur 


n/2 


F?(<—£,2)_ d[1 _ Es 
en DSG Nelkn)+z, | sinkrsingdt Te 








$ 
Der Vergleich mit (9) (S. 14) ergibt, daß die 2. Ableitung 9° F (x —£,2)/jdxd2z der Grund- 
lösung der Differentialgleichung n (Aw, +k’w,)=0 identisch ist mit der im 6. Abschnitt 


erhaltenen 1. Ableitung OF(=—£,2)/öz der Grundlösung der Differentialgleichung 
Aw, +k’w,=0, so daß wegen der Formeln (24) (S. 10) und (18) die nach beiden Methoden 
erhaltenen Lösungen für w, übereinstimmen. 


8. Bestimmung der Grundlösung durch Aussonderung mittels der Monotoniebedingung. Am 
einfachsten gewinnt man folgendermaßen die Grundlösung unseres Problems: Man geht von 
einer Grundlösung der Wellengleichung Aw, + k’w, = in der Zeilonschen Integraldarstellung 


© » 
eilax +72) 


1 

Fer, )=ReTn \ \ Pa Gel ae ae a (27) 
aus; zupächst wird diesem Integral dadurch ein Sinn beigelegt, daß bei festem y die Inte- 
gration hinsichtlich « unter Umgehung der reellen Nullstellen von k’ — a’ — y? auf dem Um- 
weg über die komplexe a-Ebene ausgeführt wird (vgl. S. 18/19); sodann wird die Wahl 
des hierbei noch unbestimmt bleibenden Integrationsweges durch die Monotoniebedingung 
eindeutig festgelegt. 

Zunächst lautet nach leichter Umformung das Zeilonsche Fundamentalintegral 








k © 
a eiax cosyz 1 IE ER 
F(x<,2)= \ | \ al yo) de dy 
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Die beiden inneren Integrale hinsichtlich a werten wir wieder mit Hilfe des Residuumsatzes aus: 
nach der Überlegung auf S. 18/19 sind die inneren Integrale gleich der Summe der Residuen 
des Integranden in denjenigen seiner Pole, welche im Inneren der oberen oder unteren Halb- 
ebene gelegen sind, je nachdem x >0 oder << 0 ist. Dabei kann jetzt aber jede der beiden 
rellen Pole a =yk?— y?, a = — yYk?—y? des ersten Integranden zur oberen oder zur unteren 
Ebene gerechnet werden. Daher ergeben sich vier Möglichkeiten für die Wahl (des Integrations- 
weges: entweder rechnet man beide Pole zur oberen Halbebene (erster Weg), oder beide zur 
unteren Halbebene (zweiter Weg), oder den einen zur unteren und den anderen zur oberen 
Halbebene (dritter bzw. vierter Weg). Bei Wahl des ersten Integrationsweges ergibt sich fast 
unmittelbar aus (28) 
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Wegen der Formeln (25) und (26) ist daher 
1 ıT 
F@&o)=7N(kn+z, \ sin (krsin t) dt. 
$ 
Auf dem zweıten Integrationsweg findet man 
k E) 
1 ee 1 \ cos y2 en 
PR dy--3—- | ——e*!?-kdy (z<0), 
zı k? — y? 2 ? __72 . 
; V y ar 
F,(x2,2)=| 
1 COS y 2 vu 
_ ser kg (x >0) 
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also =F, (— x, z); daher ist . 


1 
Fu@,)= N, (kr) +5; \ sin(krsint)dt. 
2—y 


Auf dem dritten und vierten Integrationsweg gelangt man zum selben Ergebnis, nämlich 











k E 
1 ( cosyz-sinx=yk?-- y? Ya 1 cos SEIT, pe 
\ yR—y ar. e’/r-Mdy (2<0O), 
) 
F, (2,2) = k 
1 ( eosyz-sin zyk?- y’ 3 1; cos y2 Ku 
a — dy _—_— eV Bl, 
3 | yk? —_ y? Vy’ u u d) (x >0), 
v 


also 2F,=F,+F,. Daher ist 
F, (x, 2) = N, (kr). 


Je zwei dieser drei Grundlösungen sind linear ‘unabhängig, so daß man die unendliche 
Mannigfaltigkeit ; 
1 er 
F(@,2)=cF,+(1—-o)F,=—N,(k N)+5- \ sin (kr sin t)dt 

$ 
von Grundlösungen erhält, wo c beliebig reell ist. Es gilt aus dieser mittels der Monotonie- 
bedingung eine Grundlösung eindeutig auszusondern. Zu diesem Zweck schreiben wir nach 
(11) (S. 14) die Grundlösung F für z2=0 und x<0 in der Form 





F(&,0)=ecF,(&,0)+(1—c)F,(&,0) 


oo 


rt N, (kr) 54° \ erruntgg, 
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Aus dieser Darstellung liest man unmittelbar ab, daß F(x,0)=f(r) für hinreichend große r 
dann und nur dann monoton ausfällt, wenn c=1 ist. Also ist F,(x,2) die gesuchte Grund- 
lösung in Übereinstimmung mit dem Ergebnis (9) (S. 14) der Methode des 6. Abschnittes. 


III. Meteorologische Anwendungen der Theorie. 


9. Die Aufwindverteilung über spezielle Bodenerhebungen. Die Lösung des allgemeinen 
Leewellenproblems, d. h. die Aufwindverteilung w (x,2) in einer polytrop geschichteten atmo- 
sphärischen Strömung über dem Boden mit der beliebig vorgeschriebenen Kontur 2=L, (x) 
ist nach den übereinstimmenden Ergebnissen der vorstehenden drei Abschnitte durch die 
Integralformel 


+» © 
1 1 2 1)o 
w(x,2)=2üe*?"H \voslman+t ) FARB 271 An dEl u, 
v‚—=U 


— © 
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mit —- £=rcosp, z=rsing gegeben, worin & die Anströmgeschwindigkeit, H=c’/xg=RT]jg 
die Höhe we zur Temperatur T gehörigen homogenen Atmosphäre bezeichnen und die Wellen- 


länge eue durch 
2nüc BIN: er A 
i= m I In ee 
g L RAR Vs -h1) @) 
n . 


—1g 
n R 





gegeben ist. Dabei ist y- - = und in z der vertikale bzw. der adia- 


batische Temperaturgradient. 

Die einfachste Bodenerhebung, aus der jede andere Bodenform durch hinreichend feine 
Superposition beliebig approximiert werden kann, ist offenbar das Hochplateau mit vertikalem 
Hang. Sei h die Plateauhöhe, und der Plateauhang in x=0 gelegen, so ist die Kontur des 
Plateauprofils durch 
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(x) 0 für z<0) 

'd = == 3 

” h für <>0 


gegeben. Aus der allgemeinen Integralformel (1) folgt durch partielle Integration wegen der 
Randwerte (3) und des Verschwindens der Ableitung > [-..] im Unendlichen für die Auf- 
windverteilung wp über dem Plateau (3): 





+» © 
un a 0” u 1 cos(2r’+1)o 
wp(x,2)=—2üe:/?"H \ Az T N, (k AR Y)ı..0n 1 Fa: 
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SE ı U 008 (2v+1)@ 
»/2a„H__|_ N (kr Bee y En 
2ühe ja Men+z I Jar (kr) 2, +1 


‚u 
In dem letzteren Ausdruck führen wir noch die Differentiation nach 2 aus und erhalten 
schließlich unter Anwendung bekannter Differentiationsformeln und Funktionalgleichungen 
für Zylinderfunktionen ?‘) als Aufwindverteilung über einem Hochplateau der Höhe h den ein- 
fachen Ausdruck 
z 2/2 1 Be 
wp(x,2)=—2ühker/?rH R N,(kr)sing+ r- 


v‚—_ 





Jr (kr)sin2vo x (4). 


Diese Darstellung wird der numerischen Berechnung von w (x, 2) zugrunde gelegt. Der Reihen- 
ausdruck ist leicht zu tabellieren, da auch für die Besselschen Funktionen höherer Ordnung 
Tafeln ??) vorliegen und der durch Abbrechen der Reihen entstehende Fehler wegen der Be- 


ziehung } (kr]2)?” 1 
Jsln= ai v er Q) 
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Bild 4. Aufwindverteilung längs äquidistanter Niveauflächen bei plateauförmiger Bodenerhebung. 


mit wachsender Gliederzahl rasch abnimmt und leicht abgeschätzt werden kann. In Bild 4 
ist die Verteilung der Aufwindkomponente w längs äquidistanter Höhen dimensionslos auf- 
getragen. Um das Verhalten von w (x,2) bei festem z und großen Werten von — x allgemein 
zu kennen, gehen wir won der ersten Darstellung (12) (S. 14) aus, wonach die Aufwind- 
verteilung durch 
’ ann 9 | 1 s . 1 — kr &int 
w=2ühe?r# — | — \ sin(krsint)dt — ;— \ er""Stdgt 
d2 |2n n 

0 0 


21) Jahnke-Emde: a. a. O0. S. 212 f. oder Weyrich: a. a. O. S. 72 f. 
22) Jahnke-Emde: a.a.O.S. 242 fl. oder Rep. of the Brit. Assoc. for the Adv. of Science, London 1915, S. 29 bis 32. 








Z. angew. Math. Mech. 





24 Lyra, Theorie der stationären Leewellenströmung in freier Atmosphäre Bd. 23 Nr. 1 Febr. 1943 
o i es ; 2 2 
gegeben ist. Aus dieser Darstellung folgt für «> — © bei festem 2 wegen np = = | >0— 
ı 


die asymptotische Beziehung 
w(x,2) = nd e/’?"Hgsinkz, 
n|®| 

wonach auf der Luvseite in den Schichten O<kz<n, n<kz<2n, 2n<k2<B3n,... ab- 
wechselnd Aufwind und Abwind von sehr geringer Stärke herrscht. Der Aufwind der untersten 
Schicht ist der Hangaufwind, der nur in Hangnähe eine beträchtliche Stärke erreicht. Auf 
der Leeseite treten in Abständen von einer Wellenlänge Aufwindfelder auf, welche jedoch 
erst bei bergförmigen Bodenerhebungen stärker ausgeprägt sind. Die einfachste Bodenform 
ist diejenige mit rechteckförmigem Profil. Man erhält offenbar die Aufwindverteilung über 
einem solchen Bergrücken der Höhe h und Breite b, wenn man diejenigen beiden Aufwind- 
felder voneinander subtrahiert, welche von zwei um die horizontale Strecke b gegeneinander 
versetzten Hochplateaus derselben Höhe erzeugt werden: 


2X 
z: 


2X 


»|a 





x 
Bild 5. Aufwindverteilung längs äquidistanter Niveauflächen.bei einer Bodenerhebung; 
mit rechteckförmigem Profil der Breite Ä 2. 
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Bild 6. Aufwindverteilung längs äquidistanter Niveauflächen bei einer Bodenerhebung 
mit rechteckförmigem Profil der Breite A. 
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b b 
we + 3,2) mp (2-2) N ee 


Man bekommt dabei eine optimale Verstärkung der von den beiden vertikalen Berghängen 
herrührenden entfernteren Leewellen, wenn die Bergbreite gleich der halben Wellenlänge 
le. 


7-7 ist, und eine stärkste Löschung der entfernteren Leewellen, wenn die Bergbreite ge- 


u ‚_2n u ) 
rade eine Wellenlänge 4 = z beträgt. Es ist bemerkenswert, daß sich gerade in letzterem 


Fall das erste Aufwindfeld besonders stark ausbildet (vgl. Bild 5 und 6). Die drei Bilder 4 | 
bis 6 sind bezüglich des Längenmaßstabes dreifach überhölit. Dagegen stellt Bild 7 die Auf- 
windverteilung durch Linien konstanter Aufwindstärke unverzerrt dar. 


1. Aufwindfeld 


Hangaufwind 








2.Feld 3. Feld 
7. Aufwindfeld 
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Bild 7. Darstellung der Aufwindverteilung nach Bild 5 und Bild 6 durch Linien gleichen Aufwindes. 


Die durch den Punkt (x, 2) = (— %, 2,) hindurchgehende Stromlinie ist durch die Gleichung 


x 


; 1 
2=2,+{(x.2) mit a2) \ A 
gegeben. Für 2,>0 konvergiert diese Stromlinie gegen die Kontur der beliebig abteilungsweise 
stetig angenommenen Bodenerhebung: es ist lim Z(x,2,)={,(x) an allen Stetigkeitsstellen. 
20>)0 


In den Bildern 8 und 9 ist der Verlauf der nach (4) berechneten Stromlinien wiedergegeben, 
welcher den Aufwindfeldern der Bilder 5 und 6 entspricht. Aus Bild 9 ist zu ersehen, daß 
für die Bergbreite von einer ganzen Wellenlänge nur eine Leewelle auftritt. Die sich aus 
Formel (4) ergebenden maximalen ‘Werte für w innerhalb der ersten Aufwindgebiete liegen 
im Bereich der beobachteten Aufwindstärken. 

Die Formeln (4) und (4a) liefern trotz der Linearisierung der Randbedingung (S. 8) die 
richtige Verteilung des Hangaufwindes. Denn setzt man in den Ausdruck (4) die für kleine 
r gültige bekannte asymptotische Beziehung 


k r 2v 
%) 81 
” ApR (k r) == &»)! und N, (k r) = — = Er 
ein, so erhält man als asymptotische Darstellung der Aufwindverteilung (4) für die Umgebung 
des Störzentrums r =0: 
hei .. 
y—uhsinp_ü 2 
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Bild 8. Stromlinien um einen Höhenzug mit rechteckförmigem Profil der Breite A/2. 
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Bild 9. Stromlinien um einen Höhenzug mit rechteckförmigem Profil der Breite A. 


Daraus ergibt sich für die Linien gleichen Aufwindes in Hangnähe: 
üh \? üh \’ 
er ( 7 sen) " Pa 
‘was mit den Ergebnissen der bekannten durch die Erfahrung bestätigten Theorie des Auf- 
windes an Plateau- und Berghängen übereinstimmt). Die vorliegende allgemeine Aufwind- 
theorie entlıält also die Theorie der Hangaufwinde und der bodennahen Strömung um Hügel 
als Spezialfall, und beschreibt sowohl die kleinräumige Hangströmung als auch die großräumige 


atmosphärische Leewellenströmung durch die eine Formel (4) bzw. (1) für beliebig vor- 
gegebene Bodenformen. 


28) Vgl. z. B. H. Koschmieder a. a. 0. $ 34. 
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10. Druckverteilung am Boden und Rotoren. Mit der Lösung für die Aufwindverteilung w 
erhält man die übrigen vier Störgrößen u, p, o und T durch einfache Integrationen. Eliminiert 


0 _[d ) 
man 2 aus (6) und (9) (S. 6) und ersetzt as®+n 2 durch —o +) zufolge (8), so 
tindet man für die Horizontalkomponente der Störgeschwindigkeit 

x 
RE: PER... N 
u(x,2)= \ ( w— 3e)dx. a RS ET 
. Dann ergibt sich aus (6) (S. 6) für den Stördruck p= ii u oder 
BE 
__wp([(9, oO w . 
p(x,2)= len - 5.) de Nee a 2? Fe 


Schließlich entnimmt man o und T aus (7) und (10) (S. 6). 


In dem Fall, daß das Profil der Bodenstörung auf den isolierten Punkt =z=0 kon- 
zentriert ist, läßt sich das Integral für « und damit auch für p, o und T durch stark kon- 
vergente Reihen darstellen, welche nach Zylinderfunktionen fortschreiten. Zunächst erhält 
man nach dem vorigen Abschnitt wegen der Linearität der Leewellengleichung die Aufwind- 
verteilung über einem Bergrücken mit rechteckförmigem Profil der Breite 5b und Höhe h durch 
die Superposition (4a) zweier von gleich hohen Plateaus erzeugten Aufwindfelder. Ist F=hb 
die Profilfläche, so lautet das Aufwindfeld für das Bergprofil 


el + 2 ) ( ® 
b b p'Üp x+5,2)—- mp 5,2 
wp %+75,2)— wp e—- Zı2 =, — 7 a re a 


worin wp durch Formel (4) gegeben ist. Der Grenzübergang b>0 liefert als Aufwindverteilung 
über einer auf die y- Achse konzentrierten Bodenerhebung mit der infinitesimalen Profilfläche F 
den von h unabhängigen Ausdruck 
ae F Di) wp ( 
"T > EL ER TE EEE 


Danach lassen sich in der Tat « und p nach (6) und (7) und damit auch o und T geschlossen 
durch den Ausdruck (4) für wp darstellen. 


Auch für die Verschiebungen £ der gestörten Stromlinien über einem schmalen Berg- 
rücken mit der Profilfläche F ergibt sich nach (5) und (9) die von h unabhängige Darstellung 


F 
We) pr, 2). N ee en a 
Danach ist der Stördruck (7) für schmale Bodenerhebungen durch 
P(x,2)_ ürFrl(g Od mp 


gegeben. 
Von praktischem Interesse ist die Druckverteilung am Boden 2=0. Aus (11) folgt 
wegen wp (x,0)=0 für <+0 und Formel (4) für die relative Störung des Bodendrucks 





DEM, lt anne 
u SEE TH 32\7 N,(kr)sinp + z 1 J,,(kr)sin2» o) RR 
„1 
. . Ku h) 419 ß) 
Nach Ausführung der Differentiation = — | oder =— — ‚Je nachdem x 
d2l.-0 r |dplp=0 r dply=a 














positiv oder negativ ist, ergibt sich schließlich der folgende einfache Ausdruck für die relative 
Druckstörung am Boden: 





SL. SELL ö y EN 
Di er gH «| |4 N,(kle))+ (&v)? — Fr (k En) nr 101 a 2 
v=1 











worin das positive oder negative Vorzeichen gilt, je nachdem & >0 oder << 0 ist. 
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Bild 10. Verlauf der relativen Druckstörung am Boden bei schmalem Bodenhindernis. 


Dieser leicht zu tabellierende -Ausdruck ist in Bild 10 für die Profilfläche F= 3.10° m? 
dargestellt. Aus diesem Barogramm geht hervor, daß am Boden horizontale Druckunterschiede 
von 5mb, imb und 0,6 m b unter den ersten drei Aufwindfeldern bestehen, welche wegen 


des Faktors = mit wachsender Entfernung vom Bodenhindernis stark ab- 


nehmen. P. Raethjen hat zuerst darauf hingewiesen, daß sich aus solchen abnehmenden 
Unterschieden des Bodendrucks die Erklärung für die Entstehung der sog. „Rotoren“ ergibt. 
In der Tat bildet sich wegen der Luftreibung in Lee der Bodenerhebung zunächst eine bodennahe 
Schicht mit kleiner Horizontalgeschwindigkeit aus. Da hierdurch die vorwiegend von höheren 
Schichten dem Boden aufgeprägte Druckverteilung (12) nicht verändert wird, so werden die über 
den Stellen maximalen Bodendrucks befindlichen Luftmassen in Richtung des vorangehenden 
kleineren Druckminimums, d. h. entgegen der Richtung der Grundströmung beschleunigt, 
Wenn diese vor jedem Druckmaximum beginnende Luftströmung- über die Stelle maximalen 
Drucks hinaus in das Gebiet des Druckanstiegs eintritt, so muß sie zum Stillstand kommen 
und aus Kontinuitätsgründen nach oben ausweichen. Damit gelangt die Strömung in die 
höheren Schichten mit größerer Horizontalgeschwindigkeit und wird also in Richtung der 
Grundströmung mitgeführt, um sich jeweils über den Stellen maximalen Bodendrucks zu 
einem Wirbel zu schließen. Danach stellen die „Rotoren“ nur eine Begleiterscheinung der 
Leewellenströmung dar, die sich ohne eine Voraussetzung von Grenzflächen lediglich aus der 
theoretischen Druckverteilung erklärt. 436 
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Charakteristikenverfahren für Kugel- und Zylinderwellen 
reibungsloser Gase. 
Von Robert Sauer in Aachen. 


Es wird ein auf der Charakteristikentheorie beruhendes rechnerisch- 
zeichnerisches Näherungsverfahren entwickelt für die nichtstationäre isen- 
tropische Gasströmung mit Kugel- oder Zylindersymmetrie. 


1. Einleitung. 


Bei Vernachlässigung der Reibung und Wärmeleitung und Voraussetzung konstanter 
Entropie führt sowolıl die ebene wie auch die kugel- oder zylindersymmetrische nichtstationäre 
Gasströmung auf eine Potentialgleichung vom hyperbolischen Typus. Man kann sich infolge- 
dessen des Charakteristikenverfahrens von J. Massau') bedienen, was in verschiedenen 
Arbeiten geschehen ist. Das weitestreichende Verfahren dieser Art für sehr allgemeine 
Voraussetzungen stammt von G. Guderley. Das ebene Problem wurde insbesondere von 
R. Sauer?) und F. Schultz-Grunow.), das kugel- und zylindersymmetrische Problem 
neuerdings ebenfalls von F. Schultz-Grunow*) behandelt. 

In der vorliegenden Abhandlung wird ein weiteres Charakteristikenverfahren für 
kugel- oder zylindersymmetrische Gaswellen entwickelt, das einem rechnerisch-zeichnerischen 
Näherungsverfahren zur Konstruktion achsensymmetrischer räumlicher Überschallströmungen °) 
nachgebildet ist. Der Kürze wegen wird das Verfahren hier nur unter den oben angegebenen 
engen Voraussetzungen und unter Beschränkung auf vollkommene Gase mit konstanten 
spezifischen Wärmen besprochen. Es läßt sich jedoch leicht nach verschiedenen Richtungen 
verallgemeinern. 

Die rechnerisch-zeichnerischen Näherungsverfahren für die nichtstationäre Gasströmung 
sind von erheblicher praktischer Bedeutung, da man nur wenige strenge Partikulärlösungen ®) 
des Problems kennt und diese Lösungen zur Erfüllung der Randbedingungen bei praktischen 
Aufgaben selten ausreichen. 


2. Potentialgleichung und Charakteristiken. 


Die Strömungsgeschwindigkeit. « ist eine Funktion der ÖOrtskoordinate x und der 
Zeitkoordinate ft; x bedeutet bei der ebenen Strömung den Abstand eines veränderlichen 
Querschnitts von einem festen Bezugsquerschnitt, bei der zylinder- oder kugelsymmetrischen 
Strömung den Abstand von der Zylinderachse bzw. dem Kugelmittelpunkt. 


Als Grundgleichungen hat man die Eulersche Bewegungsgleichung 


1 
1: + uu.—— > Px 


die Kontinuitätsgleichung 


4 
9 +_OUx+ 0x U vu u—0. 
und die Adiabatengleichung der vollkommenen Gase 


re 


Hierin ist o die Dichte, p der Druck und x das Verhältnis c,/c, der spezifischen Wärmen. 
Die Zeiger t, x bedeuten partielle Differentiation, der Zeiger Null bezieht sich stets auf einen 
festen Bezugszustand.. A=0 kennzeichnet das ebene, A—1 das zylindersymmetrische und 
i=2 das kugelsymmetrische Problem. 


Durch den Ansatz 


p 
“” (dp Je 
M=1=-57- — 


0 »—1 
»o 


. Massan: Ene. der math. Wiss. II, 3, 1, S. 162. 
. Sauer: Ing.-Arch. Bd. 13 (1942), S. 7989 und Ing.-Arch. Bd. 14 (1943), im Druck. 
. Sehultz-Grunow: Forsch. Ing.-Wes. Bd. 13 (1942), S. 125 134. . 
. Schultz-Grunow: Ing.-Areh. Bd. 14 (1943), im Druck. 
Sauer: 'Theor. Einführung in die Gasdynamik, Berlin 1943, $ 21 und Luftf.-Forschg. Bd. 19 (1942), S. 148 132. 
{. Bechert: Ann. Phys. 5, 39 (1941), S. 169 202 und H. Marx: Ann. Phys. 5, 41 (1941), S. 61/88. Eine Lösung 
von grundsätzlicher Bedeutung gab neuerdings G. Guderley in Luftf.-Forschg. Bd. 19 (1942), S. 302,312. 
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mit der durch 


E 


gegebenen Schallgeschwindigkeit ce ergibt sich aus (1) bis (3) bekanntermaßen die Potential- 
gleichung 











4 
(u) px u par pyartrzdum0 a re A 





Nach (4) und (5) ist 
x»—1 
U=@x, e=0'—- 59’ —D)p, 
die Potentialgleichung (6) ist also eine nichtlineare Differentialgleichung zweiter Ordnung. 
Die zu (6) gehörige Charakteristikengleichung 
(? — u)d®? +2udtde— de’=0 
liefert die beiden reellen Charakteristikenscharen 


dx Br 
dd »+re a U ER (A), 














die Potentialgleichung (6) ist also vom hyperbolischen Typus. Da _c die Schallgeschwindigkeit 
relativ zu dem mit der Geschwindigkeit'u strömenden Gas bedeutet, sind die Charakteristiken 
(7) die Wegzeitlinien für die Fortpflanzung infinitesimaler Störungen. 

3. Transformation und geometrische Deutung der Potentialgleichung. 


Wir ordnen jedem Punkt _A des in Frage kommenden Bereichs der x, t-Ebene ein 
schiefwinkliges &, n-Koordinatensystem zu, dessen Achsen in Richtung der Charakteristiken (7) 
verlaufen (Bild 1). Aus den Transformationsgleichungen 


z=£cosa+tncosß, t=!sina+nsinß, etga=u+c, ctgß=u—c 
erhält man dann die Differentiationsregeln 
WRITE WERE | FRE. SARSRUR 
Jg way,tsinay;  Pruend # V Edel Y £ 
durch deren zweimalige Anwendung 
Pen Pxx 608 acosß + Yxı(sinacos + cosasin B)+Y::Sinasinf, 


u. EN 
sinasinß xx 


folgt. Hiermit geht die Potentialgleichung (6) über in 


cotg a cotg P + Pxı (cotg a+ cotg P)+ Prr = Pxx lu — ce) +2 u px: + pri 








A j i 
P%n=— *usinasin VEN INT LUST | 








Natürlich gilt gemäß der Herleitung die einfache Beziehung (8) nur in dem betrachteten Punkt A 


Zur geometrischen Deutung von (8) untersuchen wir nunmehr die Strömung gleichzeitig 
in der x, t-Wegzeitebene und in der u, g-Zustandsebene, wobei wir die Komponenten u=9, 
und q=9. zum: Zustandsvektor iw zusammenfassen. Die Projektionen des Zustandsvektors w 
auf die &- und n-Achse (Bild 2) sind 

m=ueosa+tgqsina=g;, n=ucosß+qgsindß=gn. 


Die Potentialgleichung (8) ist daher gleichwertig mit den beiden Beziehungen 





dm usina sin Bdn 
Ä (9), 
dn=— eu sina sin ßd£ 
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Bild 1 (links oben). 
Berührsysteme des Charakteristikennetzes. 


Bild 2 (oben). 
Projektionen des Zustandsvektors. 


Bild 3 (links). 
Geometrische Deutung der Potentialgleichung. 


die man bei Verzicht auf die hier benützte geometrische Herleitung auch aus dert Verträg- 
lichkeitsbedingungen der Charakteristikentheorie gewinnen kann. 

In (9) ist dm der Zuwachs, den die Projektion m. des Zustandsvektors w im festen 
&,n-System des Punktes A erfährt beim Übergang zu einem Nachbarpunkt C auf der n-Achse 
mit der Koordinate dn (Bild 3); ebenso bedeutet dn den Zuwachs der Projektion n beim 
Übergang zu einem Nachbarpunkt B auf der &-Achse mit der Koordinate d£. 


4. Vergleich des kugel- und zylindersymmetrischen Problems mit dem ebenen Problem. 


Durch u=u(z,t), q=g(x,t) wird die x,t-Wegzeitebene auf die w, g.Zustandsebene 
abgebildet. Wir setzen voraus, daß diese Abbildung für die in Frage kommenden Bereiche 
punktweise eineindeutig sei. Dann entspricht dem Kurvennetz der x, t-Charakteristiken 
der Wegzeitebene in der Zustandsebene ein Kurvennetz von u, q-Charakteristiken; beide 
Netze sind durch die Gl. (7) und (9) miteinander verknüpft. 

Beim ebenen Problem spezialisieren sich nun (8) und (9) mit A=0 zu 


In=0; dm=0, dn=0. 


Dies besagt, daß hier die Komponente m beim Fortschreiten längs der n-Achse und die 
Komponente n beim Fortschreiten längs der &-Achse ungeändert bleibt, daß also die Linien- 
elemente der x, t- und «, gq-Charakteristikennetze paarweise aufeinander senkrecht stehen. Eine 
nähere Untersuchung’) zeigt, daß das «, q-Netz von den Randbedingungen unabhängig ist und 
daher ein für allemal fest vorliegt. Das x,t-Netz ändert sich mit den Randbedingungen, 
bleibt aber stets in orthogonal-reziproker Beziehung zu dem festen u, g-Netz und läßt sich 
auf Grund dieser Beziehung sehr einfach näherungsweise konstruieren. 

Beim kugel- und zylindersymmetrischen Problem mit A=2 bzw. A=1 geht die ortho- 
gonale Beziehung zwischen den beiden Charakteristikennetzen verloren. Beide Netze sind 
hier mit den Randbedingungen veränderlich, müssen also jeweils beide konstruiert werden. 


5. Näherungskonstruktion der Charakteristikennetze. 


Wir ersetzen die beiden Charakteristikennetze durch Streckenzugnetze und haben dann 
die Aufgabe, ausgehend von den vorliegenden Anfangs- und Randbedingungen die beiden 
Netze Schritt für Schritt zu konstruieren, indem wir die Differentialgleichungen (9) als 
Differenzengleichungen umdeuten und in geeigneter Weise auf die Streckenzugnetze anwenden. 
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In den zu konstruieren- 
den Streckenzugnetzen sollen 
die Knotenpunkte einander 
zugeordnet sein. Wir be- 
zeichnen die entsprechenden 
Knotenpunkte mit P,, P/ 
und mit ] ik» M;x die Mittel- 
punkte entsprechender Seiten 
P;,Pı, P/’Py (Bild 4. 2;x 
, bedeutet den Wert der 
a Funktion 


4 
L 


- 


Pr = 
am 


c’usinasin 





für das Punktepaar M;x, M;'x; 
durch M;; ist x festgelegt, 
M;'; bestimmt u,q und da- 
Bild 4. Näherungskonstruktion der Charakteristikennetze. mit gemäß (4) und (7) auch 
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u“"—(»—1)g und cotga=u+c, cotgß=u-—ec. 


REEL VOR 
e=c, 


Die Zuordnung der beiden Streckenzugnetze wird nun durch folgende Vorschriften 
definiert: 


a) Von den beiden durch einen Punkt M;', des u, gq-Netzes festgelegten Winkeln a, ß 
bestimmt der eine die Richtung der entsprechenden Seite P; P, des x,t-Netzes von der 
Länge d£;, bzw. dn;x, der andere liefert die Richtung für die Projektion dn;; bzw. dm;x 
der Seite P;’ P%’ des u, q-Netzes. f 


b) Die Größen d&;x, dy:r, dm;x, dn;; und @;, sind verknüpft durch die als Diffe- 
renzengleichungen gedeuteten Beziehungen (9). 


Auf Grund dieser Vorschriften lassen sich die beiden Charakteristikennetze im Anschluß 
an vorgegebene Anfangs- und Randbedingungen schrittweise Knotenpunkt für Knotenpunkt 
konstruieren. Dabei hat man im wesentlichen immer wieder folgende Aufgabe durchzuführen: 


Vorgegeben sind zwei sich entsprechende Punktepaare P,, P,'; P,,.P;,', gesucht ist das 
Nachbarpaar P,, P,'. Bei der Lösung dieser Aufgabe geht man zweckmäßig iterativ vor. 
Zunächst schätzt man den Punkt P,’ und ermittelt hieraus die Größen a,,, Qys, Bus: Pas» 
Us, Uys, Cs, Cys, durch die dann auch der Punkt P, und die Längen d£,,, dn,, festliegen. 
Das Einsetzen von Q,,, Ay, d&,, dn,, in (9) liefert die Projektionen dn,,, dm,,, mittels 
deren sich ein verbesserter Punkt P,’ konstruieren läßt usw. 


Gegenüber dem ebenen Problem?) besteht die Mehrarbeit bei dem hier für Kugel- und 
Zylinderwellen' beschriebenen Verfahren einerseits in der zusätzlichen Rechenarbeit für die 
Ausdrücke 42;, und andererseits in dem Auftragen der Projektionen. dm;x, dn;z. Außerdem 
verursacht die Iteration einen erhöhten Zeitaufwand, der sich jedoch praktisch als unerheblich 
erweist, da man gewöhnlich P,’ leicht schätzen kann und daher meist mit dem ersten Ite- 
rationsschritt zum Ziel kommt. 


6. Zusammenfassung. 


Ausgehend von der Potentialgleichung der kugel- oder zylindersymmetrischen nicht- 
stationären isentropischen Gasströmung ohne Reibung und Wärmeleitung wurde ein 
rechnerisch-zeichnerisches Charakteristikenverfahren entwickelt, durch das sich im Anschluß 
an vorgegebene Anfangs- und Randbedingungen die Charakteristikennetze in der Wegzeitebene 
und der Zustandsebene näherungsweise konstruieren lassen. Während beim ebenen Problem 
das Netz in der Zustandsebene ein für allemal festliegt und nur das Netz in der Wegzeitebene 
von den Anfangs- und Randbedingungen abhängt, sind beim kugel- und zylindersymmetrischen 
Problem beide Charakteristikennetze mit den Anfangs- und Randbedingungen veränderlich. 
Der mathematische Grund für dieses verschiedenartige Verhalten liegt darin, daß nur bei der 
ebenen Strömung die Potentialgleichung durch eine Legerdre-Transformation linearisiert 
werden kann. 451 














(A 


genü 
eine 


entw 
digk: 
Nähe 
Gleic 


wo d 


aus ( 
denk 
hinte 


Von { 
versce 


wo F 
einers 


Der \ 


diesen 
Kurve 





Li u | 


@ 9 


vo5Bson 


n 
ß 
e 


ıe 


Tr 
rt 








1 


XUM 


2. . Math. . r r f 
BEN. Febr 100 Krahn, Die Janzen-Rayleighsche zweite Näherung 33 





Die Janzen-Rayleighsche zweite Näherung der kompressiblen 
Strömung um ein beliebiges Profil. 
Von E. Krahn in Göttingen. 
(Aus dem Institut für theoretische Aerodynamik der Aerodynamischen Versuchsanstalt Göttingen.) 


Es wird gezeigt, daß man die Geschwindigkeit der kompressiblen Unter- 
schallströmung an einem Profil in zweiter Näherung nach Janzen-Ray- 
leigh berechnen kann; wenn man die Ableitung der Abbildungsfunktion 
des Profils auf einen Kreis und die Geschwindigkeitswerte der inkompres- 
siblen Strömung nur an der Profilkontur kennt. 


Das Potential 9 der kompressiblen ebenen adiabatischen Strömung, das der Gleichung 


1 
Prxt Yun (PrPıxt2PxPuPYzut Pr yon) - u et re 


genügt, wo c die Schallgeschwindigkeit ist, kann man nach Janzen und Rayleigh') in 
eine Reihe nach Potenzen der Machschen Zahl 


U 
c 


M= 
entwickeln, wo U die Anströmgeschwindigkeit und c,„ die ihr entsprechende Schallgeschwin- 
digkeit ist, und erhält zur Bestimmung der zweiten Näherung 9, von @, wenn man als erste 
Näherung das Potential @, der inkompressiblen Strömung um dasselbe Profil betrachtet, .die 
Gleichung 


Mm: 
Pızzt Pıyu= Te (Pox Poxxt290x Puv Poxu + You Povn) - en ° 


wo die Indices x und y partielle Ableitungen nach den Ortskoordinaten bedeuten. 
In der x y-Ebene sei ein Profil gegeben, das durch die Abbildung 
a 


e=ertiy- [Ott +++ ER 


aus dem Kreise mit dem Radius R um den Nullpunkt der {=£-Hin-Ebene entstanden zu 
denken ist. Die Bogenlänge der Profilkontur, gerechnet von einem festen Punkt, z. B. dem 
hinteren Staupunkt, bezeichnen wir mit s, die Richtung der inneren Normalen mit n. 
Die inkompressible Strömung setzen wir als bekannt voraus. 
Y, genügt dabei folgenden Bedingungen 
a En a Er 
2. am Profil ist 9,n=0, 
3. im Unendlichfernen ist 9,x=U, 9%,=0 Ü’+y>o). 


4), 


Von , verlangen wir, daß am Profil 9,„=0 ist und im Unendlichen die Ableitungen von o, 
verschwinden. 
Macht man für @, den Ansatz 
M? (> (2) 
a Sum 5. ae BE ee Ba Due a ae ee u Te Ze ze Zee (5), 
wo 9%” und 9%” zwei Funktionen sind, die der ‚homogenen Gl. (4) genügen, so hat man 
einerseits die Gl. (2), andererseits 


Pıxzt Yıyy» = ET! BE - + ee 
Der Vergleich der rechten Seiten zeigt, daß z. B. 
N a ea Deere Ei 
Pr 
= Inu tn da + gr EEE 


diesen Bedingungen genügen. Die Integration ist im Außengebiet des Profils längs einer 
Kurve oder am Profil von einem Punkte P, bis zu einem Punkte P,(x,y) durchzuführen. 


1) O. Janzen: Phys. Z. Bd. 14 (1913), S. 639; Lord Rayleigh: Phil. Mag. Bd. 32 (1916), S. 1. 
3 
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Y Uns interessiert die Lösung am Profil, daher betrachten wir. 

ds weiter nur den Fall, daß P, und P, auf der Profilkontur liegen 


T 
Ei und die Integration längs dieser durchgeführt wird. 

z Wir wollen hier von einer Zirkulation um das Profil ab- 

eo nis sehen; in diesem Falle wird das Integral um das Profil ein- 


s deutig. Der Fall mit Zirkulation läßt sich auch behandeln, wie 
Bild 1. an anderer Stelle ausführlich dargestellt werden soll?). 








Bezeichnen wir den Winkel zwischen der s- und «-Richtung am Profil mit r, so läßt 
sich das Integral (8) für am Profil gelegene Punkte so anschreiben: 


1 i 
5 YxtYy)sinrds. ne 5 


a _ 
a 


Se nu 


und die Ableitungen in Tangenten- und Normalenrichtung 


E= 2 
ps = (Pix + Py) sin ar A u a aa 


1 
Pen 5 (Pax + Pay) 08 1 Pe FRE Na N 


Bei der Berechnung der Ableitungen von p, beachten wir die Gleichungen 


Poxn = Youys \ (12) 
Poyun = —Voxs J 
und können dann schreiben 
M: 1 ee ’ 
Pu=n! gt )sinrtaun | gie t ei) sin rd! ra Di 
0 
s 
M° 1 Q 2 1 , 2 ng 
Yın= 77% —Pou 9 (PoxtYoy)CcoSst — YPoxs Sei teisinds)! Ru... 


v0 


Um 9,„ am Profil zum Verschwinden zu bringen, addieren wir zu p, eine Lösung der homogenen 


Gleichung 
a EEE 


Die Summe genügt natürlich wieder der inhomogenen Gl. (2). Die Lösung der Potential- 
gleichung (15) bestimmen wir, indem wir mit Hilfe der Abbildungsfunktion (3) zum Kreise 
übergehen, dort die Potentialfunktion bestimmen und dann wieder zum Profil zurückkehren. 

Zu dem Zwecke multiplizieren wir (14) mit |f’(£)| und entwickeln das Produkt in eine 
trigonometrische Reihe nach dem Kreiswinkel ©. (Den Punkten des Profiles entsprechen die 
Punkte des Kreises £ = R ei®.) 


FON dsinkO+crcoskd '. .........(6). 
£ k 


Das ist eine Radialgeschwindigkeit am Kreise, der eine Tangentialgeschwindigkeit 


ESEL en: 


“ 
entspricht. Diese müssen wir jetzt auf das Profil übertragen und haben in der Summe von 
(13) und dem durch |f’(£)) dividierten Ausdruck (17) die gesuchte Lösung unserer Aufgabe. 
Bezeichnen wir hier m# 9, die den Randbedingungen genügende Lösung von (2), so haben wir 
m: N © 
Y9s = U? Far a (P x + Pu) sınr+ Pys \ > (P5 x + Pu) sind s\ 


: 1 
+ 2 [3x 008 k O-+exsin KO, - ra) 





2) Erscheint demnächst in der „Luftfahrtforschung‘*. 
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Damit ist die gesuchte Lösung der Aufgabe — die II. Näherung der Geschwindigkeit der 
kompressiblen Strömung am Profil — ausgedrückt durch 1. die von der Profilform abhängige 
Größe r, 2. die Ableitung der Abbildungsfunktion |f’(£)|, wobei deren Werte nur für die 
Punkte der Profilkontur erforderlich sind, und 3. die Geschwindigkeit ,, der inkompressiblen 
Strömung am Profil, denn die in (13) und (14) vorkommenden Werte @,x und 9,, sind die 
Projektionen von @,; in Richtung der x- und y-Achse, lassen sich also durch 9,;, cos r und 
sin r ausdrücken. 

Man sieht, daß die Lösung der inhomogenen Gl. (2) am Rande des vorgegebenen Profils 
sich auf die Lösung der homogenen Gl. (4) am Rande des Profils zurückführen läßt, ohne 
daß man die Werte im Außenraum zu kennen braucht. 472 


Berechnung der Auslenkung und Spannungen beim Kippen 


des geschlossenen Kreisringes. 
‘Von Karl Federhofer in Graz. 


Die bekannten Untersuchungen über die Kippstabilität des Kreisringes 
beschränken sich auf die Berechnung des für den Eintritt des Kippens 
maßgebenden kritischen Druckes, sie lassen aber die beim Kippvorgang 
eintretende Formänderung und damit auch die Spannungen im Ringe 
unbestimmt. Die vorliegende Arbeit bringt die bisher fehlende Berechnung 
der größten Auslenkung und der Spannungen des Ringes. 


Zur Beurteilung der Kippsicherheit des geschlossenen Kreisringes ist die Kenntnis der 
Größe jener gleichförmig über den Ringumfang verteilten radialen Drucke notwendig, bei 
deren Überschreitung das stabile reine Druckgleichgewicht des eben bleibenden Ringes nicht 
mehr bestehen kann, so daß der Ring aus seiner Ebene kippt und einen Achter bildet. Der 
für die Instabilität maßgebende kleinste kritische Druck ist bereits bekannt. A. Hencky‘') 
berechnete ihn unter der Annahme, daß die Drücke auch nach der Verzerrung zum Kreis- 
mittelpunkt gerichtet bleiben, S. Timoshenko?) nimmt an, daß die Drücke ihre ursprüng- 
liche Richtung bei der Kippung beibehalten, also lediglich parallel verschoben werden, die 
Rechnung des Verfassers?) fußt auf der Annahme, daß die Drücke nach der eingetretenen 
Formänderung in den Richtungen der Hauptnormalen der zu einer Raumkurve verzerrten 
Ringachse wirken. Bezeichnet A die Biegungssteifigkeit des Ringquerschnittes für die in die 
Ringebene fallende Querschnitts-Hauptachse, a den Ringhalbmesser, so beträgt der kleinste 


kritische Druck p=ah, wo der Beiwert a abhängig ist von der jeweiligen Annahme über 

die Richtung der Außendrücke nach der Formänderung des Ringes; er ergibt sich in der 
1: 

i+4’ 1739’ 3, worin A das Verhältnis der Bie- 

gungssteifigkeit A zur Drillingssteifigkeit C angibt*). 


obigen Reihenfolge dieser Annahmen mit 


Da das Knicken des Ringes in seiner Ebene durch einen kleinsten Außendruck ®) =: 


bewirkt wird, wo B die Biegungssteifigkeit für die zur Ringebene senkrechte Querschnitts- 
hauptachse ist, so ist für Ringe, bei denen A< B ist, die Kippgefahr größer als die Gefahr 
der ebenen Einbeulung. Über die nach Überschreitung des kritischen Außendruckes ent- 
stehende Auslenkung des Kreisringes aus seiner Ebene geben die angeführten Untersuchungen 
keinen Aufschluß; denn zur Berechnung des kritischen Außendruckes genügt es, in den 
elastischen Grundgleichungen nur die kleinen Glieder erster Ordnung in den Formänderungs- 
größen zu berücksichtigen, womit die Grundgleichungen linear’ und homogen werden. Sie 
versagen daher bei der Ermittlung der beim Kippvorgange auftretenden Formänderung, die 
aus diesen vereinfachten Grundgleichungen nur bis auf einen Zahlenfaktor bestimmt ist. 

Zweck der folgenden Untersuchung ist die zahlenmäßige Ermittlung der größten 
Auslenkung des gekippten Kreisringes und die Berechnung der nach Überschreiten der 
Kippgrenze auftretenden zusätzlichen Spannungen im Ringe. Die Untersuchung wird sowohl 
für die vorhin angeführte Belastungsannahme von Timoshenko, — die die kleinste Kipplast 
ergibt — sowie für die Annahme von Hencky, die eine um ein Drittel größere Kipplast 
liefert, durchgeführt. 


1) A. Hencky: Z. angew. Math. Mech. Bd. 1 (1921), S. 454. 
2)S. Timoshenko: Z. angew. Math. Mech. Bd. 3 (1923), S. 358. 
3) K. Federhofer: Eisenbau Bd. 12 (1921), S. 293. 
4) Vgl. auch E. L. Nicolai: Z. angew. Math. Mech. Bd. 3 (1923), S. 297. 
5) J. Boussinesq: Comptes rendus Ac. Sc. Bd. 97 (1883), S. 843. 
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1. Die Gleichungen des Gleichgewichtes. 


Die Zentrallinie des Kreisringes vom Halbmesser a geht beim Kippen in eine Raumkurve 
über mit den Krümmungskomponenten x,x’ und der Verwindung r; für ein Bogenelement ds 
der deformierten Ringachse gelten mit X, Y,Z als den Teilen der äußeren Belastung je 
Längeneinheit des Ringumfanges in den Richtungen der Hauptnormalen, Binormalen und 
Tangente folgende Gleichgewichtsgleichungen *) 


y IN’ . A 
ırd+XeNe..... de Ga -TatY =-Nt.... (a), 
dr b2 dH Pr ’ < 
ds — Nx+Z=--N'x. .. . (1b), ds Gx =—— Wu. 
‚ 1@G 
T+N =Hx-Gr. .. (lie), Ga tHr—N=Gr “nn. (2). 


N, N’, T bedeuten die Spannkraftkomponenten in den vorhin angegebenen Achsrichtungen, 
G und @’ die Biegungsmomente für die Haupt- und Binormale, H das um die Tangente 
drehende Drillungsmoment. 


Die Spannungsmomente sind mit den Krümmungskomponenten und der Verwindung 
durch die Gleichungen verknüpft 


G=Alx —%,), G'=B(x —x,), H=C(-—1,), 


worin die Werte mit dem Zeiger 0 sich auf den ursprünglichen Zustand des Kreisringes 
beziehen; daher ist zu setzen 


0, a u=d 


und es bleibt 


1 
G=Ax ... @), E=Bl@-.). 2: @. Beer. 2: 
Im Ringe herrscht infolge des konstanten Außendruckes von Haus aus eine konstante 
Ringspannkraft T7,=— pa, so daß die Spannkraft T nach Eintritt des Kippens gesetzt 
werden kann 


a Eu a | 
wo T, den kleinen Zuwachs während des Kippens angibt. 


Mit dem Auskippen des Ringes ist eine Störung der Kreisgestalt der Ringmittellinie 
verbunden, die in erster Näherung durch die unendlich kleine Verschiebung v in der 
Y-Richtung (Ausweichen aus der Ringebene) und durch die kleine Drehung £ der Ring- 
elemente gegen die Ringebene beschrieben wird. Die zu dieser Formänderung gehörigen 
N’,H,@G sind klein I. Ordnung, somit sind die in dem Gleichungssatze (1) rechts vom 
Gleichheitszeichen stehenden Ausdrücke klein von II. Ordnung, woraus folgt, daß auch die 
beim Auskippen noch auftretenden Verschiebungen in der Ringebene (u radial, positiv zum 
Kreismittelpunkte, w in der Kreistangente, positiv mit wachsender Bogenlänge s) klein von 
II. Ordnung sind. Von der gleichen Größenordnung sind daher auch die mit dieser ebenen 
Biegung verbundenen Spannkräfte N, T, und das Biegungsmoment @’, womit aber die im 
Gleichungssatze (2) rechts vom Gleichheitszeichen stehenden Ausdrücke als klein III. Ordnung 
festgestellt sind. Werden diese vernachlässigt, so vereinfacht sich System (2) in das den 
eingangs erwähnten Arbeiten zugrunde gelegte Gleichungssystem, aus dem in bekannter Weise 
zwar der für den Eintritt des Kippens maßgebende Außendruck p;, aber nicht die Größe 
der Auslenkung berechnet werden kann. Entscheidend für die Bestimmung der Formänderung 
beim Kippvorgange ist somit die Beibehaltung aller jener Glieder in (2), die in vo und £ und 
deren Ableitungen klein von III. Ordnung sind. 


2. Die Belastungskomponenten X, Y, Z. 


Falla. Die Radialdrucke p bleiben parallel der Ringebene. Der Punkt P (Bild 1) des 
Kreises vom Halbmesser PO =a ist infolge der Verschiebungen u, w in der Kreisebene nach 
(P) und infolge der Auslenkung v (senkrecht zur Kreisebene) nach P’ gelangt; voraus- 


#6) A. E. H. Love: Lehrbuch der Elastizität. Deutsch von A. Timpe. Leipzig u. Berlin 1907, S. 456. 





wä 


lief 
vor 


hin 


ergil 
III. 





Z. angew. Math. Mech. < En 
Bd. Nr. ı Febr. 1943 Federhofer, Das Kippen des geschlossenen Kreisringes 37 





setzungsgemäß wirkt dann der Außendruck p in P’ parallel zu (P) O. Er hat in den 
Richtungen der Kreistangente und der Kreisnormalen PO die Komponenten p sin y und 


pcosy, wo siny= =, cosy=1 (genau einschließlich der Glieder klein III. Ordnung, da u und w 


bereits klein Il. Ordnung sind). 


Da sich die Kreistangente beim Übergang von P nach P’ in der Ringebene um y= 
2 — und aus der Ringebene heraus um s dreht, und da ferner die x 2-Ebene in P’ mit 


der Ringebene den kleinen Winkel ß einschließt, so gibt die Druckkomponente psiny in den 
x, y, 2-Richtungen die Teile 
a N dv ; 
psiny-siny, psiny-as: psiny-cosy, 
während die Druckkomponente pcosy die Teile 
pcosy, — pcosy-sinß, p<os y - sin y 


liefert. Damit folgt für die Belastungskomponenten bei Vernachlässigung der Glieder klein 
von höherer als III. Ordnung 


Z=p». a a ee ET (7a), 
j wdv 
= psinß+P 7 ds RE I 
w du r 
i=-I— tpv=r7 N a he a Eee 


Fallb. Die Radialdrücke bleiben auch nach dem Auskippen zum Kreismittelpunkt 0 


Bild 1 (links). 


hin gerichtet. 








Bild 2 (reehts). 


Der Außendruck p hat in den Richtungen der x,, Y,, 2,- Achsen (Bild 2) die Komponenten 
pcosy-cos®, —psin®, — psiny-cos®. 

Diese Komponenten geben in den Richtungen der x, y, z-Achsen die Teile 

pcosycos#, —psind-sinß,. psinycosd-y, 

—pcosycosd-sinß, — psin®, psin y cos 0-$°, 
r dv R 
Pcosycosd-y, —psind- az —psinycos®. 
Wegen 
in9—4.(14%) sin y 
sind — = + a)’ ‚=57: 


vo? 
cosal— 5: cosy—=1, 


ergıbt sich nach Zusammensetzung dieser Teile und Streichung aller Glieder von höherer als 
III. Ordnung 
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v? ©; 
x=p(1-9,)-pZsinß. a a ee (7d), 
rel )nnnt (rt) Ä 
u sind —p — a PA Er (Te), 
du v dv 


Z=Pgs Pads 
3. Die Krümmungskomponenten x, x’ und die Verwindung r. 


Um die Differentialgleichungen zur Berechnung der Verschiebungen u, v, w und des 
Winkels 8 zu erhalten, müssen die Ausdrücke für x, »’, r als Funktionen der ersteren ent- 
wickelt werden. Für x’ gilt die bekannte Beziehung 


a 


1] 
wobei yagkı®, da « und w klein II. Ordnung sind und x’ in den Gl. (2) im Produkte 


mit @ und ÄH (klein I. Ordnung) vorkommt, so ist die Beschränkung auf die Glieder klein 
II. Ordnung in x’ und y hinreichend’). Hingegen müssen x und r bis einschließlich der 


Glieder klein III. Ordnung in v und £ entwickelt werden, da sie im Produkte mit A auftreten. 


Nehmen wir zunächst ® gleich Null und 3=+0 an, so bilden die ursprünglich in die 
Ringebene fallenden Hauptachsen aller Querschnitte den Mantel eines Kreiskegels, dessen 


Krümmung für die Punkte der Ringmittellinie gleich Tr, (#- 2 ist. Faßt man nun 


die Auslenkung v als Funktion der Bogenlänge s auf, so ist = —=sin D (wo ® Tangenten- 
winkel) und die Krümmung ist gleich „En. Wegen n = 008 . r folgt für die Krümmung 


in der YZ-Ebene — sc a oder bei Beschränkung auf die Glieder klein III. Ordnung: 


_ (1 + + vl ) vll, 
Damit ergibt sich 


sin I. 
PR (1+4 5 00°) on ee ae de er a 
(Hochgestellte römische Zahlen bedeuten Ableitungen nach der Bogenlänge s.) Dem Winkel 8 
entspricht ein Drall vom Betrage ze. Nehmen wir aber $ und = gleich Null und für ein 


Ringelement v=+0, S® = Const an, so bilden die zum Element gehörigen ursprünglich in die 
Ringebene fallenden Hauptachsen ein kleines Stück einer Schraubenfläche (Steigungswinkel 
tga= vl), deren Windung gleich — ist, wofür bis zur III. Ordnung genau za — vl?) vl 
gesetzt werden kann. 
Somit ergibt sich für die Verwindung 
a el a a a 3 


Dem Zuwachs T, zur Ringspannkraft (Gl. (6)), der klein II. Ordnung ist, entspricht eine 
zusätzliche Ringdehnung e vom Betrage ey (F Ringquerschnitt), die auch gleich sein 


muß ze, wenn ds und do die Längen eines Bogenelementes vor und nach dem Kippen 


bezeichnen, 
Somit ist 
T, 
1+ EF 


Ey 
Qa 


(12). 


?) Vgl. die Gl. (15a) und (15b). 





Irma 


Er: 
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Die Richtungskosinusse des Bogenelementes do mit den Verschiebungsrichtungen u, v, w 
betragen 
du wds dv ds dw wuds 
do  ado’ do’ do 'do ado' 


Da die Summe ihrer Quadrate gleich Eins sein muß, so ergibt sich 


do\’ /du  w\’ /dv\’ dw u\ 
(a) ts) +Ha) ++ 2). 
oder bis auf kleine Glieder III. Ordnung genau 


() =1+vl?+2 (m -*) 


und demnach in Verbindung mit (12): 


T 
vl? u 2 (m 4) = 25F . . . . . . . . : . . » . (13). 


Mit T,=0 folgt hieraus als Bedingung für eine dehnungslose Ringachse 
or +2 (mi #)=0 ae NET 
Die Benutzung dieser Gleichung, welche die unmittelbare Berechnung von u aus w 
und v gestattet, vereinfacht die weitere Rechnung erheblich; sie hat zur Folge, daß alle 


A 
Glieder, die den Faktor Fo enthalten, gegenüber den davon freien Gliedern vernachlässigt 


werden, was für dünne Ringe, bei denen der Trägheitshalbmesser v5 klein gegenüber 


dem Ringhalbmesser a ist, zu rechtfertigen ist. Mit dieser Vereinfachung wird 


v=—+a(mi+oton) NE N te Re 


und damit 


= (IHM) + a [wit (ot)? + ot ol] ee 0 EN 


4. Die Grundgleichung zur Berechnung der Auslenkungen v. 


Fall a. Die Grundgleichung für die Auslenkung v ergibt sich aus dem Gleichungs- 
system (2, a bis c) auf folgende Art: Wir tragen dort für T und Y und x’ die ihnen gemäß 
Gl. (6), (7b) und (15b) entsprechenden Werte ein und schaffen jene Glieder, die klein von 
höherer als der I. Ordnung. sind, auf die rechten Gleichungsseiten und bezeichnen sie ab- 
kürzend mit 7‘, /’,, 7,; damit verwandelt sich das System (2) in 


. (16a), 
. (16b), 


. (16c) 


T=N(4+ + T,( oe). 2 


de 


n=-@(£ on) +6 yi ae 1 0 5 


a 


n=6 (+0) Ay. De a a re een. 
Ersichtlich sind alle /' klein von III. Ordnung. 
Wird (16c) differenziert und dabei auf (16a und b) Bedacht genommen, so ergibt sich 


+ LG-paui=Ni+1,-a=T. a a 
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Schreiben wir die Gl. (10) und (11) wegen sin ß .s-E in der Form 


u. (£ on) Ba er a rate, a Me 
=(#+ ze), TEN EEE a EIER: a 
wo 3 
y= F + ol BE (21a), 
klein III. Ordnung 
„zo . (21b) 
so folgt aus (19) und (3) 
E-® +ou+y, 22), 
und aus (20) und (5) H 1 
PZN 


oder wegen (16c) 
ER EEE ur 
C a 2 a Yı- 
Wird hieraus #! mit Hilfe von (22) beseitigt und Gl. (18) beachtet, so ergibt sich für das 
nn G der Ausdruck 
Ir 
2 a+9=a + (1422) + (4r- GT. ‚rltapl). et 


worin p 
au > ee ee a .._ 


Hiermit führt Gl. (18) zur gesuchten Grundgleichung für die Auslenkung v (s) in der Form 
a 1+4 d? 1\/a 1 
a ter av + li a) Hr-eralEl-Fr-Ntan)=e 
oder mit ds=adp: 
d’v pa’\ d'v pa’\ d’v 
Ta + rl nu ame ee  \ 
_ worin 2) klein III. Ordnung in v und deren Ableitungen ist. Vernachlässigt man die 
kleinen Glieder von höherer als I. Ordnung, das heißt, wird 2()==0 gesetzt, so vereinfacht 
sich Gl. (25) in die von S. Timoshenko angegebene Grundgleichung, deren Lösung, da v (p) 
eine periodische Funktion von g mit der Periode 2r sein muß, mit 
v(p)=csinmp 
angesetzt werden kann, wo m eine ganze Zahl ist. Das kleinste m ergibt sich mit m =2 
(m =1 kennzeichnet eine Drehung des unverformten Ringes um einen Durchmesser). Wird 
die Lösung v (%)=csin29 in die durch 2(p)=0 homogen gewordene Grundgleichung (25) 
eingetragen, so ergibt sich daraus der kleinste kritische Druck für den Eintritt des Kippens zu 
94 
Pk = a’ (A+4) (26). 
Die Auslenkung v (9) =.csin2% und die Verdrehung £ bleiben wegen 2(p)=0 unbestimmt. 
Das Näherungsverfahren zur Ermittlung der Konstanten c wird in (5) entwickelt. 


Fall b. Mit Benutzung des in Gl. (7b) angegebenen Wertes für die Belastungs- 
komponente Y und Beachtung der Formeln für T und x’ verwändelt sich das System (2) in 

















dN’ 

ds -pasi_p=T, a A en en et € . (27a), 
dH 6 

ds a a a N 
. a 


mit 


ee ee Sl 





Vor 


wor 
der 


odeı 


wo 
Gl. 


=C: 


Vers 
den 
sicht 


hierı 


aus | 


und 


bestiı 


und ı 
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Die Funktionen /', und 7, bleiben in Übereinstimmung mit den früheren G]. (17b) 
und (17e). 
An die Stelle der Gl. (18) tritt nun 


Gu+LG-pari-p®=T. TE ee 4 


Für das Biegungsmoment G findet man nach dem im Falle (a) genauer beschriebenen 
Vorgange den Ausdruck 
G ‚vu u a 1 
riEDELIDEN Er Hg r-Sh- ten), 


woraus sich in Verbindung mit (29) die gesuchte Grundgleichung für die Auslenkung v (s) in 
der Form ergibt 


, 1 pa’ 1 
u Ehe Er, 
av +, (e+ A )» +5 
1+4 d? 1\/a 1 
=ZaT-(ta)lat- Gr ntan)=0, 
oder mit ds=adpg: 
d’v 


dv pa°\ d'v „pa 
a +le+ amtrı+a-n . 


d’v a? 
gro)... . 





wo %2(9) klein III. Ordnung in v und dessen Ableitungen ist. Mit 2(g)=0 vereinfacht sich 
Gl. (30) in die von A. Hencky angegebene homogene Grundgleichung; die Lösung v (p) 
=csin29 führt hier zur kleinsten Kippbelastung 


12 A 


Pk >= "At (31). 


5. Berechnung der Konstanten c in v(p)=csin2p. 


Falla. Wir berechnen zunächst die mit der Verschiebung v»—=csin2%9 verbundenen 
Verschiebungen « und w in der Ringebene aus den Gl. (la bis c), wobei es genügt, die in 
den rechten Gleichungsseiten vorkommenden Werte N’, x, r als klein I. Ordnung zu berück- 
sichtigen. 


Man erhält aus Gl. (23) 


hiermit aus Gl. (22) 


aus Gl. (20) 


Damit sind die Ausdrücke 
eN’r=f, 
—aNx=f, 
a(Hx— Gr)=axı(C— A)=f, 


bestimmt und es gehen die Gl. (1,a bisc) wegen 


du u 
X=p, Z=p 7, (Gl. 7a und e), “= (145% 


dy 


und nach Streichung von N de 


als klein IV. Ordnung über in 
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+n- -. a 
TANTE re ee FE 5° 
Tu +Na BL ao ri A 


Wird aus der ersten und zweiten dieser Gleichungen T, beseitigt und in das Eliminations- 
ergebnis der aus der dritten Gleichung berechnete Wert von N eingetragen, wobei wegen 


Gl. (4) und (8) -4% zu Setzen ist, so kommt 








En -gllbnbE 
Zufolge Gl. (15a) ist aber 
avan+gr—4[Ö)sinte a A 
und gemäß Gl. (14) 
42T os 2p Va Be Re 


u 
dp 
womit sich die folgende Grundgleichung für die tangentielle Verschiebung w (p) ergibt 














dp* dy* dp 


Werden mit Hilfe der vorhin berechneten x, r und N’ die Funktionen f,, f,, f, be- 
stimmt. und in vorstehende Gleichung eingetragen, so läßt sich bei Beachtung der Gl. (26) 
nach längerer Rechnung die Grundgleichung in der Form darstellen 


dy* 








d’w 9 A\d'w 18 A\d’w 
dg* +let ga atli are = 900 0 (© -)sin 4 9. ur 
worin der dimensionslose Beiwert o durch 
Bi 27 (A+2) A 
IT ar B Pe 


bestimmt ist. 

Infolge Berücksichtigung der kleinen Glieder II. Ordnung in den Gl. (labisc) ist in 
der Grundgleichung (36) eine Störungsfunktion von gleicher Größenordnung (c?) aufgetreten; 
ihr entspricht die Lösung 
ne Te et Me 


wenn », wie aus der Eintragung dieses Ansatzes in (36) ersichtlich, den Wert 


a 
A 


erhält. Der dimensionslose Beiwert £ ist nur von den Steifigkeitsverhältnissen 2 und = c 


abhängig gemäß 
(40). 





tm 4 
Br 1 0 | 
"TBA+NB 
Mit der tangentialen Verschiebung w sind nach den Gl. (34) und (35) auch die Ver- 
drehung y der Kreistangente und die radiale Verschiebung u bestimmt; es ergibt sich hierfür 


ay= al -)e- Dein te, 


u(p)= (1-752)eos4P]. 





Aus G’= Bart folgt für das Biegungsmoment @’ 


@ =168(5 ) ED eos4p 


E> (es Be) En (122) Er -0(Z )sin 4p[1 )-5 (5: - 2-45) 
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A 316 A — 176 
oder wegen -1=- 5 mars 


\ ce? 3164 — 176 
G=—A a? Bür: Firz cos 4 q a Be a ee ee (41). 
Hiermit berechnet sich die Schubkraft N aus (33c) zu 


4 © 91A+49 
Nı- Ar ur FE er — ° 
Die Gl.- (33a) ermöglicht dann die Berechnung des Zuwachses T, zur Ringspannkraft — pa; 
er ergibt sich mit 
5A A > A 3161 — 176 
= - -(& Arte ) a 106 + 784 - -95 a 3 cos4p . (83). 
Nachdem nun jene Spannkräfte T,, N und das Spannungsmoment @’, die der mit dem Kippen 
verbundenen ebenen Biegung des Ringes zugehören, berechnet sind, können die Funk- 
tionen I‘, 77, T, aus den Gl. (17 a—c), sodann /’ aus Gl. (18) und schließlich die Störungs- 
funktion 2(g) in Gl. (25) ermittelt werden; nach Durchführung dieser längeren" Rechnung, 
die im Einzelnen nicht wiedergegeben wird und in der in allen Gliedern, die den Druck p ent- 
halten, hierfür genau genug p=px (Gl. 26) gesetzt worden ist, entsteht schließlich 


2(W=(5)lo,sin2r+o,sinsg] . a 


4013 + 1244 A + 366 4? +52 4° 46 
= Gars EN RER Te. 
139181 + 106195 A + 8717 4? — 785 4° 
nn A+4) 
Die Grundgleichung (25) läßt sich durch den Ansatz 


saeeraieareiae . . .,.. 5 rer 





. (45b). 


integrieren. Die Eintragung dieses Ansatzes in (25) und Umformung des Gliedes linker 
Hand liefert 
ge * 32 (384 + 143) 
e(2 -1)sin2p— aa = 5 Hk m sin69= lo, sin29+w,sin69]. 


Durch Gleichsetzen der Koeffizienten von sin 29 ergibt sich für das gesuchte c der Wert 
5 p. ie 
= N ee are ter ae 
Ve (47) 


Die weitere Koeffizientenvergleichung führt zur ia der Verbesserung c,; man erhält 


(+9)o, fe}. 
1152 88i +18) \a) '° 


Für einen schmalen Rechteckquerschnitt bA(b >h) dessen Breite 5 mit der Richtung 


des äußeren Druckes p zusammenfällt, berechnet sich das Steifigkeitsverhältnis A aus 


A _ EJ, 2(m-H1)J, . i 
ch a gj wegen G= mn! zu beahiage yo ei; mit m als Poissonscher Zahl. Da 


der Drillungswiderstand J für das sehr schmale Rechteck J -4 bh’ gesetzt werden kann 
und das Trägheitsmoment J, des Querschnittes für die zu b parallele Hauptachse gleich 


N bh’ ist, so folgt 





=— 


(48). 


ar oder mit u = 065 
2m 3 Reh 


Mit diesem A findet man 
aus Gl. (45a) w,= 496, 
aus Gl. (45b) = 2105, 
suCL) ,=— 008 (<) e 


und es lautet daher gemäß (46) die gesuchte Lösung für die Auslenkung v (p) 


2 
o(p=e| sin 29 — 005 (2) sin 6] ; 
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Hieraus entnehmen wir schließlich als hinreichende Näherung für das maximale Maß c der 
Auslenkung beim Kippen des Kreisringes mit schmalem Rechteckquerschnitt®) 


® EEE. p 


max _ 7 Pk -. (50). 





Man erkennt, daß die Auslenkung für p=p; verschwindet und daß bereits kleine Über- 
schreitungen der Kipplast zu großen Auslenkungen v„,, Anlaß geben. Ist z.B. p nur um 
1 v.H. größer als p,, dann beträgt v„., bereits 8,6 v. H. des Ringhalbmessers, bei p = 1'02 px 
ist Oyax gleich 12 v. H. von a. \ 

Die größte Änderung des Ringhalbmessers a berechnet sich aus 


uW)=Z 1+(1-152)00s49) (Seite 42) 





c? 
ZU Uni y wobei Z£ genau genug gleich 1 gesetzt worden ist. Mit Beachtung des 
Wertes von c cn. 50) folgt endlich 


un = 0788 0 (2 6 1) 


1 
Ta 1:01 Pr, Umax 10 ’max » 


Hiernach ist, wein 


und für 
: u 102 px, Unmax 013 Unax- 


Fall b. Mit dem die Grundgleichung (30) unter der Annahme 2(9)=0 befriedigenden 
Ansatz v(p)=csin2%_ ergeben sich bei Beachtung des in (31) angegebenen ersten kritischen 
Druckes px für x,ß,t und N’ dieselben allgemeinen Ausdrücke wie im Falle a, natürlich 
hat die Konstante c darin einen anderen noch zu ermittelnden Zahlenwert. Die in 5 (Fall a) 
eingeführten Funktionen f,, f,, f, sind im Falle b wegen der geänderten Belastungskomponenten 
X,Z (Gl. 7d und f) zu ersetzen durch 


(= tan c+palgs+"L), 


v dv 


f,=-aN’s+p— a dp’ 


h, =a(Hx — G?). 


Die Gl. (33a—c) bleiben erhalten, die aus ihnen in gleicher Weise wie im Falle a abzu- 
leitende Grundgleichung für die tangentiale Verschiebung ®w unterscheidet sich von Gl. (36) 
nur in einigen Beiwerten und lautet 
d*w 12 A\d'w 24 4 
dg‘ ++ at li i+iB 





ne =9600 (2 -) sin 4 9 2 se 





worin der dimensionslose Beiwert o abweichend von (37) durch 


_j_ A152+1% 
LFI OBHOÄa+N 


bestimmt ist. Der Gl. (50) entspricht die Lösung 
w(p)=»sin4p, 





2 
wo wie im Falle a: ‚=-5[&)% während {= er 4 hiermit ist auch die Ver- 
" BüÜ+4B 
drehung y der Kreistangente und mit dieser das Biegungsmoment G’ bestimmt, für das sich 
ergibt 
c? 829 1 — 644 


@=—A a Tara cost pP. 


8) Formal gleiche Ergebnisse gelten beim Kippen eines-geraden Stabes (K. Federhofer, Z. angew. Math. Mech. 
Bd. 6 (1926), S. 43), beim Knicken des Stabes (v. Mises, Z. angew. Math. Mech. Bd. 4 (1924), S. 435) und bei der Einbeulung 
eines Kreisringes (K. Federhofer, Verhandlg. d. 2. Intern. Kongr. f. Techn. Mechanik, Zürich 1926, S. 388). 
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Die Gleichgewichtsgleichungen (33c und a) liefern sodann für die Schubkraft 
4A c® 1544 +31 


N=- ga lrep 


sindg 
und für den Zuwachs T, der Ringspannkraft — pa: 
s e\ ,3254+4) =) A GE - . 
= -(&) (A +4)? (= 7T5(A+4)° ee 5 +4 
Nach Berechnung der Funktionen 7‘, (Gl. 28a), 7‘, und 7’, (Gl. 17b u. c) gewinnt man schlie&- 


lich für die Störungsfunktion 2(g) der Gl. (30) einen formal mit (44) übereinstimmenden 
Ausdruck, worin aber die Beiwerte ®, und w, (abweichend von 45a, b) ‚bestimmt sind durch 


A829 A — 644 . 
’ cos4p. 


2534 — 3714412324 754° 
(+4) , 
177600 + 128424 A + 13733 2° — 671 4° 
(+4) i 
Die Integration der Grundgleichung (30) mit Hilfe des Ansatzes 


W,= 


oa 


ve(p)=cesin?2@+ e,sin6g 
liefert nach Koeffizientenvergleichung 
Pr Ga / p 
 Vof Pr 
GrYyo,_ (<) 
3360 134+8)\a) ° 


Für den Kreisring mit schmalem Rechteckquerschnitt folgen daraus wegen A=065 die 
Sonderwerte 


und 


= 


2 
0) 744 wc > C 
0,=235, w, = 2653, ce, = —- 0% (2) c, 


so daß sich gemäß (51) die Auslenkung ® (p) mit 
2 
v(p)=e| sin 29 — 0'06 (&) sin6g | 
und die größte Auslenkung beim Kippen hinreichend genau mit 


= #198 Bu; 
Pk 


) 
® max 





ergibt.. 
Die größte Änderung des Ringhalbmessers «a berechnet sich aus 

16 c? 

Umax = 15 a 


bei Beachtung des Ergebnisses (52) zu 


Una = 1536 a (2 ı) 
Pk 


6. Berechnung der größten Druckinanspruchnahme des gekippten Ringes. 


Nun sollen noch die Größtwerte der Biegungsmomente G, @’ nnd der Ringkraft T, die 
im gekippten Kreisring bei kleiner Überschreitung des kritischen Druckes Pr wirksam 
sind, zusammengestellt und aus ihnen die größte Druckinanspruchnahme des Ringes be- 
rechnet werden. 


Falla. Aus den Gl. (6) und (43) folgt, daß die größte Ringkraft 7,., an den Stellen 
r=l@n+1) T (n=0,1,2,3) auftritt mit dem Betrage 
c? A h 


= Bo — Pr Bär 4) 1350 A + 784 (7) 


3161 — em 
A+4 


Setzen wir für den bei Überschreitung des kritischen Druckes px» wirksamen Außendruck 


em: . BE re ET 5 
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worin » >1 ist, so läßt sich der Ausdruck für Ta, bei Beachtung der Gl. (26) und (50) und 
mit A= 065 vereinfachen in 


REER pa 14 (1-—)(1167 00155) ee 


An den Stellen der Größtwerte von T erhalten gemäß den Gl. (41) und (32a) auch die 
Biegungsmomente @’ und @ ihre Größtwerte, und zwar 
c? 316 4 — 176 c 12 
G’ max = a? 3(A+4)? ’ Gns=Airy 
wofür wir mit Benutzung von (26), (50) und mit A=065 setzen können 


Gay =002 pa (1). BI We a ae ed a 


Gas = eew-ı a A 


Die größte Druckinanspruchnahme o„., des gekippten Ringes mit Rechteckquerschnitt b-h 
beträgt daher 


= (1-—)[v1e7 - 0:04 „ron (1-2) HB ueit Br 





Nun berechnet sich aber der kritische Druck p; aus Gl. (26) mit A -7 bh’ und A=0%65 zu 


b h\’ 
Pr= gg # (2) 
woraus für die in Gl. (57) als erster Faktor rechts stehende Beanspruchung des Ringes bei 


reiner Zusammendrückung 


ER 5 
sel): - 22222220. 
folgt. Tragen wir diesen Wert in Gl. (57) ein und bedenken, daß diese nur für kleine Über- 
schreitungen der Kipplast, also für wenig von Eins abweichende Werte »(>1) und für einen 
schmalen Rechteckquerschnitt (kb) gültig- ist, so kann Gl. (57) ersetzt werden durch die 
einfachere & 


1 ER 
on = 5 8) 515 Vr-i+v(n) ES a Re are 
Diese Formel gilt bis zum Erreichen der Proportionalitätsgrenze o, des Ringmaterials und 


liefert mit ons, = op für den Größtwert des Verhätnisses = hinreichend genau 


(2) = 12 0 
ya TW-—i E' 


Für Stahl folgt hieraus mit op = 2500 —, kg E= 2200000 mn: 


em?’ 


ML 


(eV 
TUR 733 33 yr — —y’ 


so daß die obere Gültigkeitsgrenze der Formel (59) 


für »=101 bei (&) —=0014 


max 


für »=10 Sei (=) — 0010 
max 


liegt. Der Vergleich der Formeln (59) und (58) zeigt, daß das Anwachsen vou 0a, bei Über- 
schreiten der Kippgrenze außerordentlich rasch erfolgt. 
Fall b. Mit den hier gültigen Zwischenergebnissen 
2 
Tan palı+{i tler +01% )|; 
Gar = 0089 p a’ (1 er =), 


2 
Ga, =124 8 v1 
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entsteht für die größte Inanspruchnahme des gekippten Ringes die Formel 


a FAZ N‘ 744 a rer 
P 1+(1- „(2176 + 023 )+ „vrzil, 


Fe bh , b v 


die aus den vorhin angeführten Gründen durch die einfachere 


Ba a 
Omas = 4.68 E() 44 yv—1+v ) 
ersetzt werden kann. Bei Erreichen der Proportionalitätsgrenze op = o 


o S 0,625 Op 
a = vv —ıE 


ergibt sich 


max 


h i i # j 
wonach der Grenzwert (=) fast um die Hälfte kleiner ist als im Falle a. 


max 


Zusammenfassender Bericht. 


Die konfluente hypergeometrische Funktion 
mit besonderer Berücksichtigung ihrer Bedeutung für die 
Integration der Wellengleichung in den Koordinaten eines 
Rotationsparaboloides. 
Von Herbert Buchholz in Berlin. 


(Mitteilung aus dem Zentrallaboratorium für Fernmeldetechnik der AEG.) 


In dem vorliegenden Bericht werden die wichtigsten Eigenschaften der 
beiden, gewöhnlich mit M;,m (2) und W;,m (2) bezeichneten konfluenten 
hypergeometrischen Funktionen besprochen und dabei eine Zusammen- 
stellung der hauptsächlich für sie geltenden Formeln gegeben, um vor 
"allem dem Physiker und dem Ingenieur das Arbeiten mit diesen beiden 
noch wenig bekannten, physikalisch aber doch recht wichtigen Funktionen 
zu erleichtern. Im besonderen wird auf ihre Bedeutung für die Integration 
der Wellengleichung in den Koordinaten eines Rotationsparaboloides hin- 
gewiesen. Die im Zusammenhang damit erforderlichen Reihenentwicklungen 
beider Funktionen bei rein imaginären Werten von k und z werden im 
Laufe der Arbeit hergeleitet. Bei derartigen Werten von k und z besitzt 
die Funktion M;;,m (2) für ein festes z unendlich viele einfache Nullstellen 
in bezug auf k. Das aus dieser Eiyenschaft ableitbare, einem bestimmten 
Rotationsparaboloid zugeordnete System von Eigenfuhktionen ist orthogonal. 
Fir den Ausdruck exp (ikR)|R einer fortschreitenden Kugelwelle wird 
eine Integraldarstellung angegeben. Ein Schrifttumsverzeichnis über die 
wichtigsten Arbeiten beschliefit den Bericht. 


1. Der Zugang zur konfluenten hypergeometrischen Funktion 
von der mathematischen und physikalischen Seite her. 


Die gewöhnliche hypergeometrische Funktion von Gauß, die durch die für alle |x)<1 


ı 


absolut konvergente Reihe (1.1a)') definiert ist, stellt bekanntlich die in der Umgebung von 


1) In der hier gewählten, von Barnes stammenden Bezeichnungsweise der hypergeometrischen Funktion gibt 
der vordere Zeiger die Zahıl der Zählerparameter und der hintere Zeiger die Zahl der Nennerparameter an. Diese 
zusätzliche Angabe empfiehlt sich immer dann, wenn man es gleichzeitig auch noch mit der verallgemeiuerten hypergeo- 
metrischen Funktion 


& 

grand (Eon tapdn xu 
u la dad td me! 
n=V 


pFa ter, any rer ap; da, da, rer dus a) 


zu tun hat. Von der sie definierenden Reihe gilt allgemein, daß sie in |x| unbeschränkt konvergiert für p<g und 
nirgends konvergiert für y24+?2. Für p=g +1 hat sie den Konvergenzradius 1. 
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‚F,(a,0,;ß.;2)= (a, is ja mit (a),—=1 und (a„=a(a+1)---(a+n—1) (l.la), 


n= 


x-.(1—x)- Hu - (,+a,+1)x]- eh 


n. 
x=0 reguläre Partikularlösung der hypergeometrischen Differentialgl. (1.1) dar mit den drei 
außerwesentlich singulären Stellen 0, und 1. Zu ihnen gehören als charäkteristische 
Exponenten für die Stelle =0 die Werte O0 und 1— ,, für die Stelle >»a, und a, und 
für die Stelle 10 und $,—a,—a,. Soll die dritte singuläre Stelle nicht bei <=1, sondern 
bei 2=c liegen, so braucht man in den obigen Gleichungen lediglich z=z/ce zu setzen. 
Die charakteristischen Exponenten werden dadurch nicht geändert. 

Bei der konfluenten hypergeometrischen Funktion handelt es sich um den Fall, daß die 
beiden singulären Stellen 2=c und z2>x» der durch die Substitution e—=z/c modifizierten 
Funktion (1.1) „konfluentieren“, d. h. zusammenströmen oder verwachsen. Die unmittelbare 
Durchführung des Grenzüberganges c>x an der Gl. (1.1) würde jedoch die Voraussetzung 
in sich schließen, daß auch nach der Konfluenz der unendlich ferne Punkt der z-Ebene noch 
eine außerwesentlich singuläre Stelle mit endlich großen charakteristischen Exponenten a, 
und a, bleibt. Die genauere Betrachtung dieser Verhältnisse an Hand der Riemannschen 
P-Funktion zeigt aber, daß in diesem Falle der unendlich ferne Punkt eine wesentlich 
singuläre Stelle wird, in dessen Umgebung die Funktion dann nicht mehr eine polare, sondern 
eine Laurent-Entwicklung besitzt. Das bedeutet aber nichts anderes, als daß bei der Konfluenz 
von c und © zum mindesten einer der beiden Exponenten a, oder a, mit ce zugleich unendlich 
groß wird. Setzt man dementsprechend in (1.1) und (1.1a) außer z=z/c auch noch etwa 
a,=c und führt erst darnach den Grenzübergang c >» durch, so entsteht für die konfluente 


PA . 
Y= ıF, (a,; ß,; 2) un; | a 


hypergeometrische Funktion die in z unbeschränkt konvergente Reihenentwicklung (1.2a), und 
an der aus (1.1) auf diese Weise entstehenden Ars (1.2) läßt sich ohne Schwierigkeit 


d’ı 
2 aa+(lı—e Ns ER RN (1.2) 


zeigen, daß sie durch die Funktion ‚F', (a,, _ tatsächlich befriedigt wird. Es stellt darnach 
diese Funktion die an der Stelle 2=0 reguläre Lösung der Gl. (12) dar. Die zweite, für 


y,=s!-fı. Fl, —-B, +1;2—-B;e) -» - -» -» . ... (12b) 


beliebige Werte von ß, von der ersten unabhängige Partikularlösung der Gl. (1.2) gibt die 
Gl. (1.2b) an. 

Der weitere Ausbau der Theorie der konfluenten hypergeometrischen Funktion ist nun 
von Whittaker [3]?) nicht unmittelbar an die Differentialgl. (1.2) angeschlossen worden, 
sondern er legte dieser Theorie aus Gründen der größeren Symmetrie, die dadurch die 
Lösungsintegrale in Abhängigkeit von den Parametern bekommen, die Differentialgl. (1.3) 


d?F 
1er re .........48 
zugrunde. Sie geht aus (1.2) durch die Substitution (1.3a) hervor, so daß zwischen den neuen 
Parametern p und »°) der Gl. (1.3) und den alten Parametern a,, £, der durch die Gl. (1.3a, ß) 


@ .. 248 
y=e’-z2 ° -Fk) . : (13a) ,=5.—-r . . (130) ß=1+p. . (1.3) 


beschriebene Zusammenhang besteht. Die grundlegenden Eigenschaften der Differentialgl. (1.2) 
sind aber durch die Substitution (1.3a) nicht geändert worden: Auchı für die Differentialgl. (1.3) 
ist nach wie vor der Nullpunkt der z-Ebene eine außerwesentliche und der unendlich ferne 
Punkt eine isolierte wesentliche Singularität.. Damit ist auch F(z) eine konfluente hyper- 
geometrische Funktion. 


2) Die Zahlen in eckigen Klammern verweisen auf das Schriftitumsverzeichnis am Schluß des zweiten Teiles 


dieser Arbeit. 
3) Für gewöhnlich werden diese Parameter mit k und m bezeichnet, so daß k=r und m=pj2 ist. 
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Um zu zeigen, daß diese Art von Funktionen auch für die angewandte Mathematik [5] 
von großer Bedeutung sind, wollen wir die Partikularlösungen der Wellengleichung in den 
Koordinaten eines Rotationsparaboloides zu bestimmen suchen. Nun macht es bekanntlich 
keine Schwierigkeit [1] die Gestalt der Wellengleichung 
in diesen neuen Koordinaten anzugeben, sobald erst ein- 
mal der Zusammenhang zwischen den Koordinaten (x, 4, 2) 
eines rechtwinkligen und rechtshändigen Bezugssystems 
und den krummlinigen und gleichlalls rechtshändigen 
Koordinaten einer Schar sich rechtwinklig schneidender, 
konfokaler Rotationsparaboloide gegeben ist. Als ein für 
unseren Zweck geeignetes paraboloidisches Bezugssystem 
wählen wir gemäß Bild 1 das System der beiden Scharen 
orthogonaler Rotationsparaboloide, das in dem zu (x, 9, 2) 
gehörenden Zylinderkoordinatensystem der (0,9,2) durch 
das Gleichungspaar (1.4a,b) dargestellt wird. Aus ihm 


ei c—d .. ... . (le) 

e=4in-n+9)..... . . (l4b 

ergibt sich zwischen den Koordinaten &, y,2 und &,n,9@ _G&ldieserAnnabel: 6l.dieser Parabel: 

der folgende Zusammenhang: g’480 (80-2) 9°-4n0 [Me*2) 
&o" 60" 


Bild ı. 








z=0:.00sp—=2YEn-cosp . . . (14a) 


y=o-sinp=2/En-sing . .  (1.4B) e=e—p (1.47). 


Über den Variabilitätsbereich von & und n gelte stets, daß OS£,n< ist. In der rechten 
Raumhälfte ist dann durchweg &>n und in der linken &<n. Die Punkte Z=n liegen in 
der xy-Ebene. Wegen anderer Einzelheiten wird auf das Bild verwiesen. Für die Ent- 
fernung R eines Punktes (0, 9,2) vom Ursprung besteht die einfache Gl. (1.4). Schließlich 


BRETTEN : rar nt, BR 


berechnen sich noch für die Längen der Bogenelemente dsz,d sn, dsp auf solchen Kurvenbögen, 
auf denen sich jeweils nur die Koordinaten £, n oder g allein ändern, die Formeln (1.4c,d, e). 


E 1 
a =(°4?) . +n 
5 7 
Damit sind dann auch die Ausdrücke für die drei Öberflächenelemente und das Volumen- 
element gegeben. Wir verzichten darauf, sie explizite anzuschreiben, sondern bringen statt 
dessen sogleich die Beziehung (1.5) für den Differentialoperator A, die nach einer bekannten 


BE : TE ER e[, Ö Ö Ö £+n ®% . 

Baht Jet tr ET BEt) (see) +3.le ” - 2En 3 N 
Umrechnungsformel [2] sofort angeschrieben werden kann, sobald erst einmal die drei Bogen- 
elemente bekannt sind. 

Bei den weitaus meisten physikalischen Aufgaben wird sich die Abhängigkeit, die die 
das erregende Feld beschreibende und unter dem Zeichen A auftretende Größe von der 
Koordinate @ zeigt, durch eine Fouriersche Reihe darstellen lassen. Es bedeutet also 
praktisch keine Beschränkung der Allgemeinheit, wenn von dieser skalaren Größe ® ein 
gesetzmäßiger Aufbau gemäß der Gl. (1.6) angenommen wird. 


cos 

sin j 

Damit nimmt dann aber die Wellengleichung, der die Funktion v (£,7) zu genügen hat, die 
Gestalt von Gl. (1.7) an mit k? in cm”? als Wellenzahl. 

vo dv ”v dv E+n 

u. a Bde u 28.8. 

era tan PO gen 

Sie ist A umgekehrt proportional der Wellenlänge mit k=2n[/A der erzeugten Schwingung. 


Um die Integration der Wellengleichung vorzubereiten, machen wir in bekannter Weise 
den Ansatz (1.8) in Gestalt eines Produktes aus zwei Faktoren, von denen jeder nur allein von 


Be ea 


4 


) ".dn. (1.44) dsp=2-(En)'r:dp. (1.4e). 


dE£. (1.4c) dy=( 


D(E,n,9)=vif,n) - pyY Bu u... ..0 A 


v+k(l+Mv=0...... (MN. 
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€ oder n abhängt. Diese Maßnahme führt bei der Funktion f, mit C als einer von & und 7 
unabhängigen Konstanten zu der selbstadjungierten Differentialgleichung: 


4 
OH + EP 0).n9=0. 


Die ganz entsprechend aufgebaute Differentialgleichung für die Funktion f, (n) enthält lediglich 
an Stelle der Konstanten C die Konstante — C. Wir ersetzen nun in der oben angeschriebenen 
Gleichung für f,(£) die unabhängige Variable & durch die neue unabhängige Veränderliche 
u=+2ik-£. Nach Division durch « nimmt damit die Differentialgleichung für f, (w) die 
Form an: ö ' a i a 
u N 
+ (+85 . + a er 
Für f,(v) gilt dann mit v an Stelle von « genau dieselbe Gleichung, wenn statt der 
Variablen n die Variable v= #2ik.n benutzt wird. Wir haben es demnach in beiden Fällen 
mit einer und derselben Differentialgleichung zu tun. Um sie auf die geeignetste Form zu 
bringen, führen wir noch zwei weitere Substitutionen durch. Die erste betrifft lediglich den 
Ersatz der von £ und n unabhängigen Konstanten € durch die gleichfalls von £ und n un- 
abhängige Konstante » gemäß der Gl. (1.9). Die zweite ersetzt die Funktion f(w) durch die 


Cl(#z2ik=v ... (9 Wu "s-F,.lw) .. . (1.10) 
Funktion F(u) von Gl. (1.10) und stellt damit für die beiden Komponenten f,(£) und f,(n) 


#F@),j 1,» , (A-plt 
da +| ır;r 2° 


Fe)=F,(w)=F,(+2ik-d). . (1.lla) F()=F,w)=F,(F2ik-n). . (1.11b) 


der Wellenfunktion (1.8) den Zusammenhang mit der durch die Differentialgleichung (1.11) 
definierten Funktion F (2) her. 

Mit der Gl. (1.11) sind wir damit auch von der physikalischen Seite her an die konfluente 
hypergeometrische Funktion herangeführt worden. Die beiden voneinander unabhängigen 
Partikularlösungen dieser Gleichung, von deren Lösungseigenschaften sogleich zu sprechen 
sein. wird, pflegt man durch die Funktionszeichen 


M, ‚Dia (2) und W, » DI2 (2) 


wiederzugeben. Demgemäß werden wir die entsprechenden beiden unabhängigen Partikular- 
integrale der Gl. (1.8a) durch die Gl. (1.12a,b) definieren 


m? (o=(nj22)'r-M,.n() . - 120) mr —(al2e): Wy,ple) - - (1.126) 


Zwei mögliche Systeme unabhängiger Partikularlösungen der Wellengleichung (1.7) werden 
dann z. B. durch die Gl. (1.132, b) beschrieben, in denen entweder die oberen oder die unteren 


v,P(En)un,P(+2ik-)-m,Pr(FRik-M. .... 0... (13a) 
v,P(En)=m,P(+2ik-8-wP(F2ik-n). ..... 0. (18h) 


Vorzeichen zusammengehören. Der Zeiger » darf in der Gleichung ohne Gefärdung ihrer: Eigen- 
schaft als Partikularlösung jeden beliebigen reellen oder komplexen Wert annehmen. Selbst- 
verständlich können in der Gl. (1.13) die Argumente £ und n auch ihre Stellungen vertauschen. 

Nach diesen vorbereitenden Bemerkungen soll zunächst die ausführlichere Erörterung 
der verschiedenen Partikularlösungen der Gl. (1.11) erfolgen. Im Hinblick auf die bereits 
vorhandene Literatur [3] werden wir uns dabei im allgemeinen sehr kurz fassen können. 
Da jedoch die bisherigen Untersuchungen über diese Funktionen kaum jemals auf eine Nutz- 
anwendung bei der Lösung physikalisch-technischer Aufgaben abgezielt haben, so sind hier 
noch eine Reihe von Bemerkungen nachzutragen, um das Auffinden eines geeigneten Lösungs- 
ansatzes zu erleichtern oder um es durch die Herstellung einer nach dem Produkt zweier 
Wellenfunktionen fortschreitenden Reihe für die physikalisch so bedeutsame Dipolstrahlung 
überhaupt erst zu ermöglichen. Außerdem wird die bisher noch gar nicht behandelte Frage 
nach den Nullstellen der Funktion m,® (z) bei rein imaginären Werten von » und 2 behandelt 
werden. Diese Nullstellen bestimmen nämlich die Eigenwerte und Eigenfunktionen, die maß- 
gebend sind für die Gestaltung der einzelnen im Innern eines hohlen Rotationsparaboloides 
möglichen, freien Teilwellen. Um die Wellenfunktionen für solche Werte von » und z auch 
der erforderlichen numerischen Berechnung zugänglich zu machen, werden im Laufe der 
Arbeit auch die dazu notwendigen Reihenentwicklungen aufgestellt. 


I F@=0. re 2 
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2. Die verschiedenen Partikularlösungen der Differentialgleichung (1.11). 


Für die weiteren Betrachtungen wird es sich als sehr nützlich erweisen, neben der mit 
der grundlegenden Gl. (1.11) identischen Differentialgleichung (2.1) auch noch über die be- 


rW+l-4 ER L eu Da a, Me 


sonderen Varianten dieser Gleichung zu verfügen, die aus ihr durch verschiedene Änderungen 
in den Variablen entstehen. Die meisten der auf diese Weise abgeänderten Gleichungs- 


u a u’ 2 „ Me 5 ©; 3 


formen ergeben sich durch Spezialisierung der Werte von a, ß und y in der Definitionsgleichung 
für die Funktion y von Gl. (2.1a), die die Differentialgleichung (2.1b) erfüllt. 


#y_ ho, +2PHtr =) .4Yy, I ,2ehtahr N, BB+D+BEO—) 
13-|2a + Frl, + - + ZT 
(2.1b). 

1— 
+0: |- rat y-o 


Für die am Ausgange unserer Betrachtungen stehende Differentialgleichung (2.1) ist bei 
beliebigen Werten von p,» und z, wie man durch Einsetzen leicht erkennt, eine erste Gruppe 
von zwei voneinander unabhängigen Partikularlösungen durch das Gleichungspaar (2.2a, b) 


1+ 
M,„,pıa (2) = 3 .. N 4 rl - w. 1+p; e) . . . . . r (2.2a) 


et ‚[1—-p . 
M,-pm(@)=z2?2 -e I „mi: „0 


gegeben. In der Tat geht aus der Differentialgleichung selbst hervor, da sie p nur im Quadrat 
enthält, daß, wenn M,,p),(2) eine Lösung ist, dies auch gewiß von M,,_», (2) zutrifft. Für 
Rep > — 1 verschwindet mit z=0 die erste der beiden lösungen und für die zweite gilt das 
für Rep<-+i. Da für p=0 die beiden Lösungen überhaupt miteinander identisch werden, 
so stellt demnach das obige Gleichungspaar für alle nichtnegativ ganzzahligen Werte von p 
in Wirklichkeit nur eine einzige Lösung mit der Eigenschaft dar, für einen verschwindenden 
Wert von z endlich zu bleiben. 

Die Funktion ‚F\, in (2.2a) erfüllt nun, wie wir schon wissen, die Differentialgleichung 
(2.1b) für ein a=+1/2, y=1 und $=—(1+p)/2. Ersetzt man in dieser Gleichung 2 und » 
durch —z und —», .so entsteht dieselbe Differentialgleichung mit a=—1/2. An Hand der 
Substitutionsgleichung geht daraus hervor, daß sich die beiden Definitionsgleichungen (2.2a, b) 
für die Funktion M,,p)(2) auch in der Form der Gl. (2.3a,b) schreiben lassen müssen. 


1+» 
M,nn (=: ? DE A 2 1+p; 2). a 


DL... 


= 1—p 
M,,-pıs (2) = 2 2 en RP +»; 1—p; 2). ee ir 


Da die beiden Gleichungspaare (2.2) und (2.3) durch bloßen Vorzeichenwechsel von » und 2 
ineinander übergehen, so muß demnach für die Funktion M,,,,,(z) die Relation (2.4a) oder 
(2.4b) bestehen. 


M,,pn (ee) erwil ri. M le 2 222. (24a) 
M_,.on (a: et")—etwrdail2. Mm) - >» 2.2... (4b). 


Nach der Definitionsgleichung (2.2a, b) ist diese Funktion, von der multiplikativen Verzweigung 
infolge des Faktors 2" +P’? abgesehen, ganz und transzendent in 2. Nur in dem besonderen 
Falle, ww» =n+(p+1)/2 und » nichtnegativ ganzzahlig ist, nimmt ihr letzter Faktor den 
Charakter eines Polynoms an. Man bezeichnet diese Polynome nach Laguerre und definiert 
sie in ihrem Zusammenhang mit der M-Funktion durch die Gl. (2.4). 


T@a+pr+!1) Me DR =, 
LP (@=7 In+1) (p+0 ‘2 BET. "Mu+2*, 3 (2)= 


T’(na+p+1) 
Tp+D) Ia+1) 





‚FRl-n;1+p; eo 
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Außer den Funktionen M,,»j; (2) hat Whittaker noch eine zweite Gruppe von zwei 
voneinander unabhängigen Partikularlösungen in die Theorie der Differentialgleichung (2.1) 
eingeführt. Da für einen beliebigen Wert von p die durch die Gl. (2.2) definierten Funktionen 
bereits ein vollständiges Lösungssystem von (2.1) darstellen, so muß sich auf Grund der 
allgemeinen Theorie jede weitere Partikularlösung von (2.1) linear und homogen durch die 
bereits bekannten beiden ersten ausdrücken lassen. Das ist nach den beiden Definitionsgl. 


I(— Rs 
W,, Pia (@)=- 1 Sp M,, pla (@)+- + P)_ a M,, -plI2 (2) nn... 
r( „) r($rP_, 


5 Gh 
p=#0,1,2... jarc(+ 2)|< 37/2 
I (— I 
WW,» (-2)= TU + x trat M_„pel-2) . . . (2.5b) 
rt) Fe) 


(2.5a,b) auch in der Tat der Fall. Gegenüber den beiden M-Funktionen haben die soeben 
eingeführten Funktionen W vor allem den Vorzug der größeren Allgemeingültigkeit. Zwar 
scheinen auch sie zunächst dem Sonderfall eines nichtnegativ ganzzahligen Wertes von p 
nicht gerecht zu werden. Die genauere Betrachtung zeigt aber, daß ihre Definition nach den 
obigen Gleichungen in diesem Falle wenigstens im Sinne einer Grenzwertgleichung richtig 
bleibt, denn für ein p>2 z.B. wird zwar das erste Glied in (2.5a) unendlich groß. Dasselbe 
gilt aber nach Gl. (2.2b) wegen des Nennerparameters 1—p auch für das zweite Glied. 
Es stellt sich also für solche Werte’ von p die obige Lösung der Differentialgleichung (2.1) 
in der unbestimmten Form © — © dar, und man kann, wenn man es wünscht, von hier aus 
lediglich mit Hjlfe der einfachen, für solche Fälle gültigen Regel sogar bis zu dem expliziten 
Ausdruck für W,,,ı, vorstoßen. 

Um die Vorstellung zu festigen, geben wir diesen Ausdruck nachstehend in Gl. (2.6) 
wieder. Er gilt für jedes p=0,1,2---. Für ein p=0 fällt dabei in (2.6) die erste, endliche 


er 
ea ee - 3.00 (3 4) 
er Fr —j [ ..,. 70.8 


W,s(@)= I 0 8 
u FE Re 2 





=u 


2 


» r/[p+t1 a 
Hiper) en, |Pa+9+Ya+i4 p)- v5? P4a-v)-Ine| 
1=0 


jarc2|< 37/2, p=0, 1, 2 .--, absolut und unbeschränkt konvergent. 





Summe ganz heraus. Bequemer läßt sich die Beziehung (2.6) aus der Integraldarstellung 
der W-Funktion nach Barnes gewinnen. Als charakteristisch ist demnach für diese 
Funktion bei ganzzahligen Werten von p das Auftreten eines logarithmenbehafteten Gliedes 
anzusehen. Auch das entspricht durchaus der allgemeinen Lösungstheorie. 

Die einzige Ausnahme von dieser Regel bildet allein der Fall eines »=n+(p+1)/2 
mit nichtnegativ ganzzahligen Werten von n. Bei Annäherung von » an diesen Wert gehen 
zwar die beiden Gammafunktionen vor der geschweiften Klammer gegen unendlich. Das 
gleiche trifft aber auch zu von der Gamma- und Psi-Funktion in der zweiten Summe mit 
dem Argument (p + 1)/2 +4 — », wenigstens solange als A kleiner ist als n+1. Die Funktion 
W,,p., entartet daher unter dieser Voraussetzung in ein Polynom n-ten Grades, und wir 
werden noch sehen, daß es sich auch hier um das Laguerre-Polynom von Gl. (2.4) handelt. 
Der symmetrische Aufbau der beiden Gl. (2.5a,b) zeigt im übrigen, daß für alle Werte von 
s,v und p stets die Beziehung (2.7) bestelit. Kommt p auch in » vor, so darf dieses p den 
Vorzeichenwechsel natürlich nicht mitmachen. 


W; v,‚+pı2 (z) = W; »,-Pia (2) A . ee a . . .. (2.7). 


Da für ein p=0, 1,2... jede der beiden Funktionen W in der Nachbarschaft von 
#=0 die Eigenschaften beider M-Funktionen auf sich vereinigt und für ein p=0, 1,2. 
zum mindesten ein Glied mit: In (2) auftritt, so besitzt die zuletzt eingeführte Gruppe der 
Partikularlösungen für p+0, 1 oder v#n+(p+1]j2 die Eigenschaft, für 2>0 über alle 
Grenzen zu wachsen. Davon abgesehen stellt sie aber für jedes p ein in sich selbst voll- 
ständiges Lösungssystem dar, da sich die beiden Funktionen (2.5a) und (2.5b) in der Nähe 
von 2=0 oder z2> ganz verschiedenartig verhalten. Neben der Gl. (1.13) gibt daher auclhı 


v,P(En)=n_,P(F2ik-E-w4,P(F2ik-n) (v beliebig) . . . . . (28) 
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noch die durch Gl. (2.8) wiedergegebene Kombination von Wellenfunktionen eine mögliche Lösung 
der Wellengleichung an, in der wiederum die gleichgestellten Vorzeichen zusammengehören. 
Nach der allgemeinen LSsungstheorie muß es natürlich außer den beiden Beziehungen 
(2.5a,b), die die Funktionen der zweiten Lösungsgruppe durch die der ersten darstellen, 
auch noch eine Beziehung geben, die das Umgekehrte leistet. Diese Relation kann rein formal 
aus den angegebenen beiden Gleichungen dadurch gewonnen werden, daß man in (2.5b) die 
M-Funktionen mit dem Argument —z mittels. der Gl. (2.4b) durch M-Funktionen mit dem 
Argument +z ersetzt und sodann das neu entstandene Gleichungspaar nach einer der 
beiden M-Funktionen auflöst. Auf diese Weise entsteht die sehr wichtige Beziehung*) (2.9). 
M,,p12 @)= As: „2 .etrir. W_,psla-et")+ 
T'\: — ») . r 
(-3212< arc2< +5) + 
Den. 
;  W,,pır(2) 
B+7,,) 


Ist darin entgegen der besonders angeführten Voraussetzung der Bereich von arc z auf das 
Intervall -n/2< arez< +3n]/2 beschränkt, so ist in (2.9) das Vorzeichen von i im Exponenten 
der drei e-Funktionen von + in — umzuwandeln. 





3. Die allgemeinen Integraldarstellungen für die Funktion M,,,.. (2). 


Den bequemsten Zugang zu der ersten und ältesten Form der Integraldarstellungen für 
die beiden hier in Rede stehenden Funktionen M und W, die im wesentlichen von Whit- 
taker selbst stammt, gewährt die bekannte Methode von Laplace. Ihre Anwendung auf 
die Lösung der Differentialgleichung (2.1) bietet zugleich die einfachste Möglichkeit, den 
inneren Zusammenhang und die gemeinsame Wurzel dieser beiden Funktionen bloßzulegen. 
Geht man gemäß der Vorschrift dieser Methode mit dem Ansatz (3.1) in die Differential- 


P»—1 p—1 


v I. 


B 
Er ang 
get Ferre: n -w—1) ° ER ea ae a ar De 
A 


gleichung (2.1b) für die besonderen Werte a=1/2, $=(1+p)/2 und y=A=1 ein, so zeigt 
eine einfache mit einer partiellen Integration verknüpfte Rechnung, daß der gewählte Ansatz 
tatsächlich das Gewünschte leistet, falls in den beiden Endpunkten des Integrationsweges der 
v-Ebene für den Ausdruck 
e+!,, o+i 
e-?’.» 2 :w—1) 2 


— + 


die Differenz der dort herrschenden Werte verschwindet. Dabei ist es bekanntlich gleich- 
gültig, ob dieses Verschwinden eintritt, weil die beiden Glieder der Differenz einander gleich 
werden oder weil jedes Glied für sich zu Null wird. 

Der allgemeinsten dieser Bedingung entsprechenden Lösung, bei der das Verschwinden 
der Differenz aus dem zuerst genannten Grunde erfolgt, läßt sich durch eine geeignete Ver- 
fügung über die willkürlich wählbare Konstante die Gestalt von Gl. (3.2) geben. Das darin 


Pp+D „{l- ee WEHEN 
Mg reger). 2 .e ( 
(14.04,1-,0-) 
Be, 3... 


re, T -(#—1) 2 -dv O<v;<1 


VA 


p—1 


2 +) 





vorkommende Integ-al hat, um eine Beschränkung der für die Parameter » und p zulässıgen 
Werte möglichst weitgehend zu vermeiden, als Integrationsweg eine Doppelschleife, auf der 
im Falle der Gl. (3.2) zu Beginn des Weges zwischen den Punkten O0 und larcv=0 und 
arc(—1)=—n zu setzen ist. Mit Ausnahme der negativ ganzzahligen Werte für p dürfen 
sonst in der Tat » und p in (3.2) jeden beliebigen reellen oder komplexen Wert haben. Für 
geeignete Werte von » und p kann das Doppelschleifenintegral (3.2) auch als ein nur den 

%) In dem bekannten Buch von Whittaker-Watson [3] ist in der Angabe dieser Beziehung ein Druckfehler 


unterlaufen. Vornehmlich aus diesem Grunde wird der Herleitung der richtigen Beziehung aus den Integraldarstellungen 
im Abschnitt 4 ein breiterer Raum gewidmet, als es sonst nötig wäre. 
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Punkt 0 oder den Punkt +1 umlaufendes, einfaches Schleifenintegral geschrieben werden. 
Ist z. B. Rev+p/2) >—1/2, so kann die Darstellung der Funktion M,,,;, unter denselben 
Anfangsbedingungen für arcv und arc(v— 1) wie oben auch durch das einfache Schleifen- 
integral der Gl. (3.2a) erfolgen, da unter dieser Voraussetzung der Integrationsweg in den 


1—p (+) 





Bu }) 


re + +) sn \ at 2" 
‚2ni. re? 


NM,» (= ‘2 a y w—1) ? don Bu: (3.2a) 


+») ' 


Nullpunkt hineinlaufen kann. Ist umgekehrt die Bedingung Re(p/2— r) > — 1/2 erfüllt, so 
läßt sich aus der Doppelschleife eine vom Punkt 1 ausgeliende einfache Schleife machen, die 
lediglich den Nullpunkt der v-Ebene umschlingt. Treffen schließlich die beiden letzten Un- 
gleichungen gleichzeitig zu, so entsteht das gewöhnliche Linienintegral (3.3), das sich noch 


1 
+1 
T(p+V T.rn.\ 


Be LIFE Sn ö 


p—1 v rn 
1—u) ? u 2? .e"?.du (3.3), 





be | 
+1 
P(p+V I 703 en um 
M,,pn. (2) = ru aan GE Ba ı ne 5 (5) ni (3.3a) 
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symmetrischer in der Gestalt der Gl. (3.3a) schreiben läßt. Diese letzte Darstellungsform 
wird sich als besonders geeignet erweisen, wenn » und z rein imaginär sind. _ 

Auf Grund der Tatsache, daß nach Gl. (2.2a) die Funktion M,,»,(2) mit der ganzen 
transzendenten Funktion ‚F, in Zusammenhang steht, können zur Darstellung von M,, p1s 
natürlich auch solche Integrale dienen, die im wesentlichen die Funktion ‚F, beschreiben. 
Das bekannteste unter ihnen, das in der Tat im vorliegendem Falle mit großem Vorteil be- 
nutzt wird, findet sich bei Watson [4] und ist in Gl. (3.4) wiedergegeben. Paßt man hierin 


R + 
ven euren) > 


0 


die Konstanten b, « und y der Gl. (2.3a) an, so entsteht die zunächst nur unter der Be- 
schränkung Re (pl2+») >— 1/2 gültige, neue Integraldarstellung (3.5). Auch hier läßt sich 
die die Gültigkeit einschränkende Voraussetzung durch Übergang zu Kurvenintegralen erheb- 


T’(p+1) 


M,n.(2)=2- er rer. erw. au 2) .du 41 2 
res ) ; 


+» 


lich mildern. Das kann einmal dadurch geschehen, daß man nach dem Ersatz von u? durch « 
den Integrationsweg vom Punkte — © aus im positiven. Sinne um den Nullpunkt der «-Ebene 
herumführt. Bemerkenswerter ist jedoclı eine andere Art der Wegführung des Integrals, die von 
Meijer [6] angegeben worden ist. Er läßt mit dem gleichen Integranden wie in (3.5) von — 
her den Integrationsweg aus dem linken Winkelraum zwischen den beiden Hauptdiagonalen 
der «-Ebene oberhalb des Nullpunktes nach dem unendlich fernen Punkt des rechten Winkel- 
raums verlaufen. , Auf diese Weise entsteht dann für M,,»),(2) die mit Ausnahme der Fälle 
p=—1,—2,—3,... und2»—p=1,3,5,... für alle Werte von 2,» und p gültige Integral- 
darstellung (3.6). In ihr bedeutet « einen Phasenwinkel mit |arcu|<r/4. Beim Bestehen 
der zu (3.5) gehörigen 


@.+expiu 


Ma pH ger) 2. error. | em, u: 2) wrdu 6) 


© :exp(i(”+ u)) 


Gültigkeitsbedingung läßt sich natürlich das Integral (3.6) auf das einfache Linienintegral (3.5) 
zurückführen. Von A. Erdelyi [7] ist die Darstellung (3.6) dazu benutzt worden, um das 
Verhalten der Funktion M,,»,, für große Werte von Re» festzustellen. 
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4. Die allgemeinen Integraldarstellungen für die Funktion W,,,,. (2). 


Gemäß den Ausführungen zu Beginn des vorigen Abschnitts müssen sich zwei weitere 
Partikularlösungen der Differentialgleichung (2.1) durch die beiden nachstehenden Integrale 
darstellen lassen. Damit darin der Integrand eine in v eindeutige Funktion ist, sind bei be- 


(+) 

Be . p—1,, Br 
f,@)= a I ?. \ he -(—1) 2 AIR . (4.1a), 
».exp(it,) 

+) 
. Di .. Az 


2 o.e 2. \o-re., s .w—]1) 2 er (4.1) 


© .exp(it,) 


liebigen Werten von p und » die beiden in Bild 2 angedeuteten, von den Punkten O0 und 1 
ausgehenden Verzweigungsschnitte vorzusehen, die von den zugehörigen schleifenförmigen 
Integrationswegen umschlungen werden. Ferner bedeute in den obigen Integralen r, den Phasen- 
winkel von v» oder v— 1 zu Beginn des Weges. Dann ist die Bedingung (3.1a) gewiß erfüllt, 
solange die Ungleichung — rj?2 < r, + arc (2) <+ [2 besteht. Zugleich ist damit die Konvergenz 
der Integrale gesichert. Verlaufen also wie in Bild 2 die Verzweigungsschnitte in der oberen 
v-Halbebene, so bleiben die angegebenen Integraldarstellungen gültig für alle arc(z) mit 
—3n/2<arc(2)<+n/2. Verlaufen sie in der unteren v-Halbebene, so muß — n/2< arc (e) 
<+3n]/2 sein. 








u-1 











+” 
Ru Bild 3. Die Phasenwinkel von v und u zählen 
stets von der positiv reellen Achse aus, und 
zwar positivim mathematisch positiven Sinne. 





Mit der Wahl der Konstanten 1, =— 7’ b - _— exp (a i (v - —) wird in der heute 


üblichen Bezeichnungsweise die Funktion f, (2) identisch mit der Funktion W,,,1(2). Es gilt 
also bei einem geeignet gewählten |r,| < nr im ganzen Winkelbereich jarcze| < 3 r/2 die 
Beziehung‘ 


+r —_—+r 
Wo@)=-—g,7 — 5’ -\er2.v 2" .w-1) ? do (Al). 


(1+) 
r(v+"5®) | ei 1— p—1 e—! 
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Die Darstellung (4.1) versagt scheinbar, wenn »+(1— p)/2=0, —1,—2,... wird. Dieser Um- 
stand ist jedoch nicht wesentlich, denn für ein Re(+» — p/2)< + 1/2 darf der Integrations- 
weg in (4.1) sowieso in den Nullpunkt der «-Ebene hineinverlegt und aus der Schleife ein 
geradliniger Integrationsweg gemacht werden. Die auf diese Weise entstehenden besonderen 
Integraldarstellungen sind nach naheliegenden Änderungen der Integrationsvariablen in den 
Gl. (4.1a,b) wiedergegeben. 
@:exp(itg) 
j p—1 p>1 
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In dem zweiten Integral f,(z) setzen wir v=—u-+1 und also v— 1=-.u. Dann geht 
mit dieser Substitution, da in Rücksicht auf die zu Bild 3 angegebene positive Zählrichtung 
von u und u—1 

arc(-u+1)=arce(w)=are(w— 1)+r und arc(— uw)=arce( — 1)=are(u)+n 
zu wählen ist, das Integral (4.18) in die Gestalt der Gl. (4.2’) über. Das hierin vorkommende 
Integral entspricht aber gemäß Gl. (4.1) hinsichtlich der Führung seines Weges und des Anf- 
(+) 


p—1 B;.,, 


—r 
2, \eou.u 2 .w—1) 2 
© .+eXplitg—im) 


baus seines Integranden, dem Integranden für die Funktion -W_,,p, (— 2). Dabei ist offenbar — 2 
auf Grund der oben angebenen Konvergenzbereiche der* beiden Integrale f,(z) und f, (2) 
mit 2-e+*i zu identifizieren. Mit C,=(, ist daher 
BL 
Tr (5? pP +) 
L.et?rivr+rip—)/2,W 


f. @)= rise u ) - 


Kann wegen Re (v+p/2) > — 1/2 der Integrationsweg in den Nullpunkt hineingeführt werden, 
so läßt sich das Integral (4.1) für f,(z) in ein einfaches Linienintegral überführen, und man 
erhält dann die’Gl. (4.2a), in der » durchweg die Phase r, hat. 

p+! © exp(iza) 


(z .e*) 2 
p+1 
r( B) +r j 


Zum Beweise der wichtigen Beziehung (2.9) möge nun für ein z mit 32/2 <arc(e2)<+n/2 
der Schleifenweg in Gl. (3.2a) etwa in der Mitte zwischen den Punkten 0 und 1 in die obere 
v-Halbebene hineingezogen und in der durch den Phasenwinkel r, festgelegten Richtung ins Un- 
endliche geführt werden. Wegen-der Voraussetzung über arc (z) ist das ohne Gefährdung der 
Konvergenz des Integrales möglich. Nach der Durchführung dieses Prozesses haben wir es 
dann in der v-Ebene mit zwei getrennten Integrationswegen zu tun. Der eine dieser Wege 
besteht aus einer Schleife, die den Punkt +1 im mathematisch positiven Sinne umschlingt. 
Infolge der für Gl. (3.2a) geltenden Voraussetzung hat dabei am Anfang der Schleife arc (v—1) 
einen Phasenwinkel, der gegenüber dem im Integral für f, (2) um — 2x kleiner ist. Von dem 
dadurch erforderlich werdenden Faktor abgesehen stellt aber dieses Integral gemäß der 
Gl. (4.1) die Funktion W,,,), dar. Das zweite der oben erhaltenen Integrale besteht aus 
zwei in den Nullpunkt der v-Ebene einmündenden Wegen, auf denen sich wiederum die 
Phasenwinkel von 9» — 1 untereinander um den Winkel 2x unterscheiden. Nach Gl. (4.2a) 
ist dieses Integral durch die Funktion W_,,p), (— 2) darstellbar. Die Zusammensetzung dieser 
beiden Teilergebnisse im Sinne der an der einfachen Schleife vorgenommenen Zerlegung führt 
unmittelbar zu der Gl. (2.9), die damit für den Winkelbereich —3r/2< arc (2)< + r/2 unter 
der Voraussetzung eines Re (»+p/2) >—1/2 bewiesen ist. Nach dem Prinzip der analyti- 
schen Fortsetzung überträgt sich dieses Ergebnis sofort auch auf den Fall beliebiger Werte 
von » und p. Geht man von vornherein von der Darstellung (3.2) für M,,»;; durch eine 
Doppelschleife aus, so läßt sich auch schon der Beweis selbst von der obigen Beschränkung 
hinsichtlich der » und p befreien. Ganz analog gelingt der Nachweis für die Richtigkeit der 
Gl. (2.9) innerhalb des Winkelbereichs — n/2<arc (2) <3r/2, indem man die Schleife in 
Gl. (3.2) in die untere v-Halbebene hineinzielt. Für ein z mit Jarc(2)|< r/2 siud daher 
zwei Deutungen möglich. 

Eine andere allgemeine Folgerung aus Gl. (4.1) betrifft den Sonderfall »=n + (p+1)/2 
mit positiv ganzzahligen oder verschwindenden Werten von n. In diesem Fall kann nämlich 
der Integrationsweg von (4.1) zu einer in sich geschlossenen, den Punkt »—=1 umschlingen- 
den Kurve gemacht werden. Die Auswertung des Integrals führt durch Entwicklung von t+1 
nach Potenzen von t=v — 1 wieder zum Laguerre-Polynom. Damit ist für diese Polynome 

(0+) 
LP («@)=(—)*- +1 FREE. R ri. HhrtHr.et2.dt. . (4.3) 
y T(n+h 42,5 Ani (tI<l) n=1,2,3...pbel. 
die Integraldarstellung (4.3) entstanden, und es hat sich nebenbei ergeben, daß sie für ganz- 
zahlige Werte von n die nachstehend angeführten beiden Darstellungen durch die Wellen- 
funktionen m und w gestatten. 


“u. 2 


met. er 


Bi... B?,. 
.ei.\ er:.0 3 .w_n 3 de. . „(42a 


W_,,p1(2- er") — 
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Ip +n+ 1) -pI2 (p) i 
BYE ee FT 
-pI2 ’ ET Rn (4.3 4, b). 


In (@)= (@]a)" 


Dr (—)" : (2Ir)' 2. Tal .ez?. Wr p4 (z) 


An weiteren Integralen für die Funktion W,,,,, (2) kommen wir darnach zu denen, die 
den Beziehungen :(3.5) und (3.6) für M,,p,, (2) entsprechen. Hier knüpft man nach °A. Er- 
delyi [8] am besten an das Integral (3.6) an. Hat p zunächst einen ganz beliebigen Wert, 
so ist es jedenfalls statthaft, die Definitionsgleichung (2.5a) für W,,,,, (2) anzuwenden. Geht 
man hierin mit der Gl. (36) ein, so treten in dem vereinigten Integral die Besselschen 
Funktionen J, und J_, gerade in derjenigen Kombination auf, in der sie die Hankelsche 
Funktion definieren. Es entsteht somit die Integraldarstellung (4.4). An ihr läßt sich auf 
Grund einer bekannten Eigenschaft der Hankelschen Funktionen erneut die Richtigkeit 
von Gl. (2.7) nachweisen. 

; @.expui 
W. En 24 il) u (1) 1 2» 
v,213(2) = 2": e eu .H, (2u:z'2)-w’-.du Jarcul<ajt.‘. (4.4). 
exp(t+u)i 

Zwei ganz andersartige Integraldarstellungen stammen von E. Barnes und dem Ver- 
fasser. Zu der ersten [3] gelangt man, indem man in der Gl]. (4.1a) das Binom (1+1)”+®-»": 
durch das Mellinintegral (4.5) ersetzt. Wegen der absoluten Konvergenz der Integrale ist 

+oi 

Er *; r| 1—p 1 ’ 1—p\ 7. 8 

(i+2) i re r(s-r+52)1 (-8):#:ds Jjaret|<n (45) 

— moi 
in dem entstehenden Doppelintegral die Vertauschung der Reihenfolge der Integrationen zu- 
lässig. Das innere Integral läßt sich dann nach einer bekannten Formel durch Gamma- 
funktionen ausdrücken. Man gelangt so unmittelbar zu der Formel (4.6). Der Integrations- 
weg ist in diesem und in den meisten ähnlichen Integralen dieser Art so zu führen, daß die 


pP.ezı 


1— 1 
rr+ 2) r(r4 52) 





W,,n, ()= 


+@i (4.6) 


—.»ı;ı 





beiden nach rechts und links ins Unendliche gehenden Polketten der Integranden, nötigenfalls 
durch Bildung von Schleifen, auf verschiedene Seiten des Weges zu liegen kommen. Die 
Beziehung (4.6) erlaubt auf besonders bequeme Weise die unmittelbare Verifikation der’ 
Gl. (2.5) und (2.6). Auch die Eutwicklung von W,,») (2) für große Werte von z läßt sich 
sehr leicht [3] aus ihr gewinnen. Wegen der vor dem Integralzeichen stehenden Gamma- 
funktionen bricht die aus (4.6) hervorgehende Reilienentwicklung bereits nach einer endlichen 
Zahl von Gliedern ab, wenn »— (1+p)/2 einer nichtnegativen ganzen Zahl gleich wird 
(s. Gl. (4.3)). 

Die vom Verfasser a. O. benutzte Integraldarstellung vom Mellin-Typus nimmt ihren Aus- 
gang von der Gl. (4.1b), und sie entsteht, wenn -man den darin vorkommenden Quotienten 
((+1)/(t—1))” nach Potenzen von t?— 1 entwickelt. Das gelingt am einfachsten mit Hilfe 
der Darstellung dieser Entwicklung durch das in Gl. (4.7) angegebene Mellinintegral für die 


t-1? (2-1)? +17” (e 1) u 
ST anal 








Po+Dr{r+5) 
(@2v+1) 


=. (1). 





rl + 55274 1;1 -®) 
Re(t)>0, |1—t|<1 

+i® 
T(—s)T(s+») r(s+r+5) 


T6+%» +1) —.(P —1)’*°.ds 
lare (? —1)|<n 
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hypergeometrische Funktion von Gauß. Geht man mit dieser Beziehung wiederum in die 
Gl. (4.1b) ein und nimmt die auch hier zulässige Vertauschung der Reihenfolge der Inte- 
grationen vor, so läßt sich das innere Integral durch die Hankelsche Funktion mit rein 
imaginärem Argument ausdrücken. Im ganzen ergibt sich auf diese Weise die Gl. (4.8). Löst 


Är- zUi2. e"ipit 


2. et) 


+wi 


ze 


W,,2:()=— 


p+l 
jarev|<r 





3 Ps Tor(lo+g)T(s+ Er pe 


1 (N 
Te r+l )- H „„6eld) ds 


2 





larc2]< nr 


— wi 


man hierin die Hankelsche Funktion in die Summe zweier Besselschen Funktionen auf, 
so läßt sich das Integral mit dem Faktor J,;}7,, sowohl durch den unendlich fernen Halb- 
kreis der rechten als auch durclı den der linken Halbebene zu einem geschlossenen Weg er- 
gänzen. Für den Faktor J_s_,j, ist dies nur möglich in der linken Halbebene. Auf diese 
Weise gelangt man zu sehr merkwürdigen Reihenentwicklungen, über die anderorts berichtet 
werden soll. Da für »=0 die beiden auseinander strebenden Polketten von /'(s) und /'(— s) 
mit ihrem ersten Pol s=0 zusammenfallen, so muß dieser Fall, da er die Hindurchführung 
des Integrationsweges zwischen den ersten Polen der beiden Ketten nicht mehr gestattet, 
aus dem Gültigkeitsbereich der Gl. (4.8) ausgeschlossen werden. Schiebt man aber den Inte- 
grationsweg noch vor dem Grenzübergang über den Pol s=0 hinweg und geht erst darnach 
zur Grenze über, so entsteht die bekannte Beziehung: 


WW. (2) = e (n2)"r . e"tr/t. Hl), 9) ER Da | | 
Schließlich erwähnen wir noch die von C. S. Meijer [9] aufgefundene Integraldarstellung 
(4.9), die dadurch bemerkenswert ist, daß sie im Integranden eine hypergeometrische Funktion 


(+) 
Wun@=TU-a:err..; | ‚ueri. Ze ;a;)-du ci 


2ni 
__ 
lare2|< nr 


mit einem willkürlichen Nennerparameter enthält, der nur verschieden von —1,—2, —3,... 
sein muß. Setzt man z. B. in dieser Gleichung mit nichtnegativ ganzzahligen Werten von n 
und p den Zeiger v=n-+(p+ 1)/2, so reduziert sich die linke Gleichungsseite nach (4.3) im 
wesentlichen auf das Polynom von Laguerre. Mit a=.tP_, entsteht dann rechts nach 
Übergang zu dem obigen Wert von » die schon früher dafür angegebene Integraldarstellung (4.3). 
Für a „122 — » erhält man nach demselben Übergang die der früheren zwar gleichwertige, 
in ihrem Aufbau jedoch verschiedene Integraldarstellung (4.10), die sich natürlich auch auf 
direktem Wege gewinnen läßt. 


(+) 


—1e)P- . \rr P-1. (I HNR.eet.di RS 


urn. 
n+p=0,1,2--- 


Tm+1) (4.10). 


LP (d)= 


Auf die große Fülle weiterer Integraldarstellungen, die vor allem von C.S. Meijer [6] 
aufgestellt worden sind, kann im Rahmen dieser Arbeit nicht eingegangen werden. Wir 
möchten hier nur noch auf zwei besonderes bemerkenswerte Ergebnisse dieser Untersuchungen 
hinweisen: Das ist die Erweiterung der Laplace-Transformation auf die beiden konfluenten 
hypergeometrischen Funktionen, die C. S. Meijer in der unter [6] zuletzt erwähnten Arbeit 
gelungen ist, und auf die Darstellung der Funktion W,, ,,, (2)/2 als die Lösung einer homogenen 
Integralgleichung mit den Grenzen 0 und & durch J. L. Sharıma [20]°). 


5) Die Abschnitte 5 bis 12 werden im nächsten Heft dieser Zeitschrift erscheinen. 
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KLEINE MITTEILUNGEN 


Numerische Berechnung von 


x 
y@a)=e-*|etat. 
0 


Zur Auswertung einer theoretischen Arbeit über 
Sekundär-Elektronen-Emission!) war die Tabu- 
lierung der genannten Funktion erforderlich. 

x 


Zwischen Null und 2 sind die Werte von | et? dt 
f) 
angegeben in der Sammlung von Jahnke-Emde?), 
ihre Berechnung kann z. B. auf dem von H. G. 
Dawson?) vorgeschlagenen Weg erfolgen. Für 
Argumente über 2 sind uns keine Tabellen bekannt 
geworden; da die Formeln von Dawson hierfür 
schwer zu handhaben sind, haben wir einen anderen 
Rechnungsgang eingeschlagen, der im folgenden 
kurz skizziert wird. Der Faktor e-*? ist hinzu- 
gefügt, weil einmal die gewünschte Funktion in 
dieser Form in der Anwendung auftrat, zum anderen, 
weil das mit €-*? multiplizierte Integral für jedes 
z beschränkt bleibt, was in der Tabulierung etwas 
bequemer ist. 
x 


Setzt man y=e-“?\el’dt, so gilt für y die 
0 
Differentialgleichung: 

Betr i: .: 


Da % eine ungerade Funktion von z und 4 (0) =0 
ist, setzen wir heuristisch die Potenzreihe an 


y=42+ 32°+ ,2°+---- 
daraus 


= .+34,2°+504,2%°-+---- . (8) 
ferner 


By 2,2? +2a,20°%+.--- 
(3) und (4) addiert gibt: 


yY+22y=a+(3a+2a,)2?+(59,+2a,)cr+--- 
=1 (wegen (1)). 


(2); 


(4); 





Danach bestimmen sich die Koeffizienten zu 
2 
= —- —4,;, 


3 
Man erhält also 


=]; 4=——G 


usw. 


2 
y#2—-2-°4z 


22: 22° 


22',232° ar 
a ae 
Wegen ihres einfachen Bildungsgesetzes, nach dem 
sich jedes folgende Glied einfach aus dem vorher- 
gehenden ableitet, ist diese für jedes z konvergente 
Reihe auch für die numerische Berechnung, wenig- 
stens für z<2 gut brauchbar. Für größere z ist 
eine andere Entwicklung vorzuziehen : 

Durch partielle Integration des unbestimmten 
Integrals | et? dt erhält man 


1) H.Schlechtweg, demnächst. 


2) E. Jahnke und F. Emde, Funktionentafeln mit 
Formeln und Kurven, Berlin 1938, S. 32. . 


s) H.G. Dawson, On the Numerical Value of | ex? dx 


0 
Proc. Lond. Mathem. Soc. Bd. 29 (1898), S. 519/22. 


et? et? dt 
f* Tr. +) 


ei? et? 
Be u 


3.et? dt 
———— usw. 


n 
BR VolElı 2. 


L, ki22k+1.g2k+1 


(2n+1)! el?dt 
a n!22n+1.fen+t2 





Wegen der negativen Potenzen von £ läßt sich 
die Grenze O0 nicht ohne weiteres einsetzen. Man 
x 
nimmt daher zunächst | et?dt und läßt v nach 0 
gehen: 


x 


+ et? .dt 


n! 22n+1.f2u+2 





u 


In das zuletzt stehende Integral wird für et? die 
au 

Reihe > er 
u! 


eingesetzt; wegen deren gleich- 


us 
mäßiger Konvergenz ist die gliedweise Integration 
erlaubt und man erhält, wenn man auch für eu? 


a» 
die Reihe > E einführt: 
= 


x n 


fe di= 


u =U 


(2k)! ex? 
kiQ2k+1 g2k+i 
n e2 
7 (2)! u2r-2k-1 
k!22k+ı ERRE \ 
=U v=U 
z?u—!n—1 
u'(2u—-2n-1) 


2n+N)! 


( ! 
+1 Br 


u= 


= 
“R u2 u—2n—1 
Y B ‚u! 2u-2n-1)' 


u=Uu 


(Zn-+1)! 
1 Ratı 





Es läßt sich zeigen, daß sich die 2. und 4. Summe 


für jede vorkommende negative Potenz von u 


identisch weghebt. ‚Für positive Potenzen von u ist 
dies zwar nicht der Fall, aber auch ohne Interesse, 
da die entsprechenden Glieder beim Grenzübergang 
u>0 beliebig klein werden. Soll die Behauptung 
zutreffen, so muß also für eine willkürlich heraus- 
gegriffene negative Potenz von u, z. B. u?e-1, 
gelten: 

—20o —-1=2u —2n—- 1=2v—2k—1 (11) 
oder 

vn=n—o und v=k—o; 





Kleine Mitteilungen est: 











0,1450 
0,1408 
0,1367 
0,1329 


S ©» @® <J 9 


_ 


0,1293 


rt 


0,1260 
0,1228 
0,1197 
0,1168 
0,1141 


oe oo o sea a a om 
Sn 


N 
- 


0,1115 
0,1090 
0,1066 


DD 


0,1043 


oc 2 so u mW 


0,1021 


2 


- 


0,1000 
0,0980 
0,0961 
0,0943 
0,0925 
0,0908 
0,0891 
0,0876 
0,0860 
0,0845 


SO © oo A 9 U a8 


so 


2 
3 
4 
5 
6 
7 
8 
9 
0 


2, 


_ 


[> 


-_ 
SO © JS u 2». vw 


— 


per 


0,0831 
0,0817 
0,0804 
0,0791 
0,0779 
0,0767 
0,0755 
0,0744 
0,0733 
0,0722 


1) 
ur DO Een 


oe eo [3290 a » 
®& 


SO ©. &@ SI 9» WW 8 8 
_ 
o 


so 
© ®8 ® SS 9 a B»8 


a m 0 m - 
Rs 



































Z. angew. Math. Mech. 
Bd. 23 Nr. 1 Febr. 1943 


Kleine Mitteilungen 61 





also für ein beliebiges fest gewähltes n und e<Sn 
soll sein: 


@n+1! rt 
n! (n—o)! 22rt1 (20 +1) 


n 
 ;; @M)! 
" ki kt (kg)! 
=Uu 


Diese Beziehung schreibt sich für n=m-+ 1: 


(2m+3)! 
(m +1! (m-e+N! 2m+3(20+1) 


m+1 


(2 k)! 
+ ZRH (ko)! 
=U 


und für n=m: 














@m+1)! \ 
m! (m — og)! 2m+1(20-+1) 


m 
(2 k)! 
Zn ZRH (kg)! 
=U 


(14) von (13) subtrahiert, gibt nach Vereinfachung 
die leicht nachzuprüfende Identität: 


ee 15) 
.4=(2m+29)!(2e-+1) un 


Da die Beziehung (12) für n=0 richtig ist, gilt 
sie durch vollständige Induktion für jedes n, womit 
die Behauptung bewiesen ist. 

Man erhält . also nach Zufügen des Faktors e”*? 
zu (10): 





(14) 








x 


y-e-x [er dt 
{0 


n 
5 (2k)! 1 
Be ki 22%k+1' z2k+1 
m 
(2 n+ 1)! e-x? “ z2u—2n—1 
m anti Le 
u= 


Die Entwicklung besteht also aus zwei Teilen, 
einer bei dem Glied mit dem Index n abgebrochenen 
(bei Fortsetzung semikonvergenten) fallenden 
Potenzreihe in z, deren einzelne Glieder eine ge- 
wisse Ähnlichkeit mit den reziproken Gliedern 
der Reihe (7) aufweisen, und einem Restglied, das 
seinerseits aus einer mit e *? multiplizierten, für 
jedes z konvergenten Potenzreihe bes:eht, deren 
erste n+1 Glieder negativ sind. 

An sich ist natürlich der Wert von % nicht davon 
abhängig, weiches n gewählt wird, oder m. a. W. 
wie viele Glieder der fallenden Potenzreihe man 
nimmt, da deren (aus lauter positiven Summanden 
bestehende) Summe durch das Restglied auf jeden 
Fall auf den wahren Wert von % ergänzt wird. 
Es ist indessen zweckmäßig, die semikonvergente 
Reihe mit dem kleinsten Glied abzubrechen, dann 
ist der absolute Betrag des Restes kleiner als 
dieses. Hierzu hat man für n die nächste unter z? 
liegende ganze Zahl zu wählen. 

Zur Abschätzung des hei günstiger Wahl von n 
zu erwartenden Restes kann man in Formel (16) 
für das letzte Glied 











(2 n)! 
nn Mntl.z2n+l 





an Stelle der Fakultät die Stirlingsche Näherung 


om. „2 +12 
„Ye 
e 


einführen und erhält mit n=r?: 
ex? 1 
z-V2 xz-.yY2 


der Rest ist also für 


An .10-0,43. x? 





(18); 


= 2 3 4 5 
kleiner als 6,7-10”°; 3,1-10 5; 2,4-10-®; 2,9-.10=2, 


Die mit dem kleinsten Glied abgebrochene halb- 
konvergente Reihe (16) liefert demnach olıne Be- 
rücksichtigung des Restes den Wert von y bei e—=2 
auf etwa 2, bei 2—=3 auf 4, bei z—=4 auf 7 und 
bei 2=5 auf 11 Dezimalstellen genau. 

Zur Aufstellung der Zahlentafel wurden die 
Funktionswerte für die durch 0,1 teilbaren x un- 
mittelbar auf 8 Dezimalen berechnet, und zwar 
bis 2=2 nach Formel (7), für ze >2 nach Formel 
(16); für x >4 brauchte das Restglied nicht mehr 
berücksichtigt zu werden. 


Essen. F. Stäblein und R. Schläfer. 439 


Über Massenpunktverteillungen kon- 
stanter Trägheit auf der Kugel. Die vor- 
liegenden Bemerkungen beziehen sich auf das 
Problem, die Masseneinheit diskontinuierlich so 
auf der Oberfläche der Einheitskugel zu verteilen, 
daß die axialen Trägheitsmomente für alle durch 
den Kugelmittelpunkt hindurchgehenden Achsen 
gleich groß ausfallen. Bei einer vorgegebenen 
Massenverteilung gibt es mindestens eine Achse, 
für die das Trägheitsmoment den größten Wert 
erreicht. Jede Lösung der oben gestellten Auf- 
gabe liefert eine Verteilung der Masse 1 auf der 
Einheitskugelfläche der Eigenschaft, daß das oben 
bezeichnete größte Trägheitsmoment seinen klein- 
sten noch möglichen Wert annimmt. Der Minimal- 
wert dieses Trägheitsmaximums ist 


Ist der Kugelmittelpunkt Ursprung eines ortno- 
gonalen Koordinatensystems, und bezeichnen wie 
üblich für eine Verteilung der Massen m 


A=!:m(y’-+ 2°), 
B=Zm (2?+ 12°), 
C=Zm(2°’+y’) 


die axialen Trägheitsmomente der Koordinaten- 
achsen, 


D=Zmyz, 
E=Zmzz, 
Fzznsg 
die Deviationsmomente, so sind die drei Momente 
,<shsh, 


der Trägheitshauptachsen, bekanntlich durch die 
drei reellen Wurzeln der charakteristischen 
Gleichung 
| A—iA—F —E | 
B-i—D Io 


C—, 


—HF 

I-E 
gegeben. Für eine Verteilung der Masseneinheit 
auf der Einheitskugel ist 


—D 
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> +y+2=1, 
A+B+C=2, 
und nach der bekannten Beziehung 


A+B+C=ht+h+%, 


schließt man auf 


also 


> 


»Sh= 3 

für irgendeine beliebige Massenverteilung. Eine 
Lösung unserer Aufgabe liegt dann vor, wenn alle 
drei Wurzeln der charakteristischen Gleichung zu- 
sammenfallen, d. h. wenn 


P) 


3 

ist. Das Trägheitsellipsoid ist in diesem Falle eine 
Kugel, und alle axialen Trägheitsmomente sind 
gleich. - Gibt es nun diskontinuierliche Massen- 
verteilungen, welche das leisten? — Es dürfte eine 
bekannte Tatsache sein, daß eine gleichmäßige 
Verteilung der Masse auf die Eckpunkte eines der 
Kugel einbeschriebenen regulären Körpers eine 
Lösung darstellt. Dies kann auf Grund geeigneter 
Sätze über die regulären Körper leicht eingesehen 
werden !). So ergeben sich Lösungen des Problems 
für die Massenpunktzahlen 


n=4, 6, 8, 12, 20. 
Wir zeigen jetzt, daß es für alle n>2 sehr ein- 
fach geartete Lösungen gibt. Damit wird eine 
Eigenschaft des regulären Vielecks im Raum auf- 
gezeigt. 
Für jedes rn >2 wird durch 


1 
m,=— 


„=== 


n 
< 
7,= v: cos .. 





eine Lösung gegeben. Es handelt sich um die 
Eckpunkte eines regulären n-Ecks, das dem 
Schnittkreis der Einheitskugel mit einer Ebene im 
Abstand 1/Y3 vom Mittelpunkt einbeschrieben 
u age Kreis hat den Radius Y2/yY3, vgl. 
i 2 





1) Vgl.H.Hadwiger: Über ausgezeichnete Vektorsterne 
und reguläre Polytope. Comm. math. helv. Bd. 13 (1940. 4)), 
S. Wu bis 107, bes. S. 101 bis 103. 


In der Tat ist in diesem Fall 


n 


. 2vn 
> RER 1 dich 
(2sin - +1) 


„==; 


Die charakteristische Gleichung hat die Gestalt 


2 3 

(5-2) =% 
und liefert so die weiter oben erwähnten zusam- 
menfallenden drei Wurzeln. 

Man kann formal n>oo gehen lassen, und 
strebt so einer Massenverteilung längs des oben 
bezeichneten Kreises auf der Kugel entgegen. 
welche eine Lösung des Problems für kontinuier- 
liche Massenverteilungen darstellt. 


Bern. H. Hadwiger. 446 


1 > / 
Der Näherurgswert ı — 3 (y141 — y6). 
Im Folgenden wird eine vereinfachte Darstellung 
des vom Verfasser gefundenen Näherungsweıtes 
für zn gegeben). 
Man zeichne den Einheitskreis um den Nullpunkt 


und den Kreis durch den Nullpunkt um den Punkt 
,=+23,9,=0 


N a 5 \ 


der den Einheitskreis in A e . V 4 schneidet. 


Mit Ai= v: als Radius beschreibe man um 


den Punkt z_,=— 1, y_,=0 den Kreis 


at +y=3. N re © 





Bild 1. 


1) Soc. Seient. Feunica, Com. Phys. Math. V, !7 (1930); 
Allg. Vermess.-Nachr. Bd. 45 (1933), 8. 421; Acta Acad. 
Abvensis, Math. Phys. VIII, 8 (1935). 
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Dieser schneidet den Kreis (1) in einem Punkt mit 
der Abszisse zg = 2. Aus einem bekannten Satz 
über rechtwinklige Dreiecke, angewandt auf das 
Dreieck OB4, folgt dann, daß 


DE DE 1 1 
TB=]/ 4.5=3 Y141 und oB=)/ 45-4 


ist. Ferner berechnet man aus dem rechtwinkligen 


Dreieck OXB die Strecke BX = 5 V 6. Somit ist 


1 — 3 (Y141— y6 )=3,1416174 - --. 
Die Konstruktion ergibt somit die Rektifikation 
des halben Einheitskreises mit einem Fehler, der 


kleiner als er - 10° * ist. 
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BUCHBESPRECHUNGEN 


MARCEL SEDILLE, Ingenieur A. M. et E.C. P., 
Docteur des Sciences, Pr&cis de me&cani- 
que des fluides & l’usage des inge- 
nieurs. VIII+154 S. m. 102 Fig.. Paris 1942, 
Verlag Dunod. Preis brosch. 105 Fr. 


Dies Buch beabsichtigt, die französischen In- 
genieure mit der modernen Strömungslehre be- 
kannt zu machen. Dabei soll der Gebrauch der Ma- 
thematik möglichst eingeschränkt sein, ohne aber 
das Strömungsfeld in ein zu enges Schema, wie 
es etwa in der hydraulischen, eindimensionalen 
Behandlung geschieht, zu pressen. Der Verfasser 
kommt bei diesem schwierigen Unternehmen da- 
durch zum Ziel, daß er sich auf die Angabe der 
Ergebnisse beschränkt, sobald die mathematische 
Herleitung zu verwickelt erscheint. Beispielsweise 
werden die Navier-Stokesschen Differentialglei- 
chungen und die durch Reibung dissipierte Ener- 
gie einfach hingeschrieben. Dafür wird der phy- 
sikalische Sinn der Formeln gründlich erörtert. 
Durch weitgehende und folgerechte Anwendung 
der Ähnlichkeitsmechanik wird der Überblick über 
den verarbeiteten Stoff wesentlich erleichtert. 
Dieser ist in der Tat sehr umfangreich. Er um- 
faßt sowohl die klassischen Potentialströmungen 
mit ihren Anwendungen auf die ebene Tragflügel- 
theorie durch Kutta und Joukowsky als auch 
die Theorie der laminaren und turbulenten Strö- 
mungen durch Rohre und über Platten, selbst 
wenn diese rauh sind. Bei der Besprechung der 
turbulenten Strömungen ist recht viel Platz der 
durch Taylor begründeten Theorie der isotropen 
Turbulenz eingeräumt. Daran schließen sich Aus- 
einandersetzungen über den Formwiderstand, den 
induzierten Widerstand und den Wellenwider- 
stand. Das vorletzte Kapitel bringt eine gute Ein- 
führung in das Verhalten kompressibler Strömun- 
gen, während das letzte Kapitel den verwickelten 
Fragen des Wärmeaustausches gewidmet ist. 

Der Verfasser bringt keinerlei Hinweise auf 
Originalarbeiten, anhangsweise ist nur eine Zu- 
sammenstellung weiterführender Lehrbücher oder 
zusammenfassender Werke gegeben. Trotz des 
Fehlens des wissenschaftlichen Apparates merkt 
der Sachkenner auf Schritt und Tritt, mit welcher 
Umsicht die neueren Fortschritte berücksichtigt 
und zu knappen Sätzen verdichtet sind. Auf diese 
Weise ist ein flüssig geschriebenes lLesebuch über 
den eindrucksvollen Aufstieg der Strömungslehre 
in unserem Jahrhundert entstanden. Freilich wird 
damit der Leser noch nicht an die Schwelle selb- 
ständiger Entwicklungsarbeit geführt, aber er be- 
kommt jedenfalls Bescheid, worum es bei dem 
Ben Aufbau der modernen Hydrodynamik 
geht. 


Dresden. W. Tollmien. 447 


Dr. phil. HANS REICHARDT, Gesetzmäßig- 
keiten der freien Turbulenz, (VDI-For- 
schungsheft 414, Beilage zu „Forschung auf dem 
Gebiete des Ingenieurwesens“. Ausgabe B, Bd. 13, 
Mai/Juni 1942.) 22 S. m. 15 Bildern. Berlin 1942, 
VDI-Verlag G.m.b.H. Preis brosch. 5 M. 


In dieser. Zeitschrift Bd. 21 (1941), S. 257 bis 264, 
hatte H. Reichardt eine neue Theorie derjeni- 
gen turbulenten Strömung vorgeschlagen, die 
unter den Begriff der freien Turbulenz zusammen- 
gefaßt werden, da sie ohne Beeinflussung durch 
begrenzende Wände erfolgen. Indem Reichardt 
die turbulente Ausbreitung des mittleren Impulses 
in Analogie zur molekularen Wärmeleitung setzte, 
gelang ihm eine mathematisch einfache und da- 
bei verblüffend genaue Darstellung der ein- 
schlägigen Turbulenzerscheinungen. Zur weiteren 
Nachprüfung seiner Ansätze hat er in dem 
hier angezeigten Forschungsheft ein umfang- 
reiches eigenes Beobachtungsmaterial zusammen- 
getragen, da die vorhandenen Meßergebnisse 
seinen Genauigkeitsansprüchen nicht. genügten. 
Schon der Aufsatz in dieser Zeitschrift brachte 
eine neue Vermessung der ebenen Strahlgrenze, 
die sich zwischen einem homogenen Luftstrom 
und der angrenzenden ruhenden Luft ausbildet. 
Dazu sind nun noch Messungen des ebenen und 
achsialsymmetrischen Kreisstrahls sowie des 
ebenen Windschattens gekommen. Bemerkens- 
wert ist dabei die gute Bestätigung der theoreti- 
schen Ansätze auch im axialsymmetrischen Fall, 
wo die Nachprüfung bisher noch ausstand. Beim 
ebenen Freistrahl wurde erstmalig auch die Strö- 
mung im Außenraum, die durch das Einsaugen 
der Luft in den Strahl entsteht, theoretisch dar- 
gestellt. 

Durch die wertvollen neuen Beobachtungs- 
ergebnisse und deren gründliche theoretische Dis- 
kussion ist dies Forschungsheft geeignet, der Tur- 
bulenzforschung neue Anregung zu geben. 
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Dr.-Ing. E. MELAN, o. Prof. a. d. Techn. Hoch- 
schule Wien, Dipl.-Ing. R. SCHINDLER, Obering. 
der Waagner-Biro A.-G., Wien, Die genaue 
Berechnung von Trägerrosten, mit 
39 Tafeln für die praktische Anwendung ausgear- 
beitet von Dipl.-Ing. R. Rothmayer. V + 134 S. m. 
101 Abb. u. 9 Zahlentafeln. Wien 1942, Springer- 
Verlag. Preis brosch. 21,10 M. 


Der Trägerrost ist die wichtigste Idealisierung 
eines Tragwerks des konstruktiven Ingenieurbaus, 
dessen Festigkeit jedoch zur Vereinfachung der 
Lösung meist ohne Rücksicht auf den zweidimen- 
sionalen Zusammenhang, also ohne Rücksicht auf 
die Kontinuität der Querträger nachgewiesen wird. 
Diese ist bisher, abgesehen von Näherungslösungen, 
nur bei der Beschreibung des Spannungszustandes 
ausgewählter Bauaufgaben mit Symmetrie nach 
zwei Achsen streng behandelt worden, da dann 
die Lösung der Aufgabe durch Belastungsumord- 
nung weitgehend vereinfacht werden kann. Die 
statische Untersuchung dieser hochgradig statisch 
unbestimmten Tragwerke bereitet formal weder 
bei frei drehbarer noch bei drehsteifer Verbindung 
der beiden sich rechtwinklig oder schiefwinklig 
kreuzenden Trägerscharen Schwierigkeiten. ‚Diese 
ergeben sich erst bei der numerischen Auflösung 
der zahlreichen linearen Verträglichkeitsbedingun- 
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gen. Der konstruktive Ingenieurbau kann sich je- 
doch gegenwärtig in seinem Streben nach wirt- 
schaftlicher Ausbildung des Tragwerks dıeser Auf- 
gabe nicht mehr entziehen. Die zahlreichen schwer- 
belasteten Rostbrücken im Zuge der -Reichsauto- 
bahn sind dafür gute Beispiele. Dies Bedürfnis 
nach einer strengen Berechnung ist nunmehr durch 
die ausgezeichnete Arbeit der beiden Verfasser 
befriedigt worden. Die Berechnung gelingt mit 
einigen plausiblen Annahmen über die Eintragung 
der äußeren Kräfte und die elastischen Eigenschaf- 
ten des Trägerrostes, ohne dabe: besondere Be- 
stimmungen über die allgemeine Anordnung und 
die Stützung der Hauptträger und über die Anzahl 
der Querträger zu treffen. Die Lösung 'st also 
umfassend und im Sinne der Baustatik streng. Sie 
erstreckt sich auf vorgeschriebene Lastengruppen 
und auf die Ableitung der Einflußlinien zur Be- 
rechnung der Grenzwerte für ausgezeichnete 
Schnittkräfte. Die statische Untersuchung der 
Verfasser führt auf ein Simultansystem von 2 ho- 
mogenen Differenzengleichungen zweiter Ordnung, 
deren Eigenwerte und Eigenlösungen allgemein 
abgeleitet und für zahlreiche Sonderfälle numerisch 
berechnet werden. Die Trennung der Veränder- 
lichen des Ansatzes führt auf die Berechnung des 
Rostes mit Hilfe eines statisch unbestimmten 
Hauptsystems, dem durch die Querträger zusätz- 
lich elastisch gestützten Hauptträger, der als ein- 
facher oder durchlaufender Balkenträger' oder als 
Bogenträger vorgegeben sein kann. Diese Unter- 
suchungen werden klar, ausführlich und _mit 
großer. Sorgfalt vorgetragen und an zahlreichen 
Beispielen für alle wichtigen Belastungsarten des 
konstruktiven Ingenieurbaus erläutert, so daß die 
Anwendung der Theorie keine Schwierigkeiten be- 
reitet. Um die damit verbundene Arbeit zu er- 
leichtern und die Zahlenrechnung abzukürzen, 
sind wertvolle Rechenvorschriften angegeben und 
Zahlen- und Kurventafeln für die statische Unter- 
suchung häufig verwendeter Trägerroste ausge- 
arbeitet worden. Damit ist eine nicht nur statisch und 
mathematisch interessante und lehrreiche, sondern 
auch für die künftige Entwicklung des konstruk- 
tiven Ingenieurbaus wertvolle Arbeit entstanden, 
welche die befriedigende Lösung eines wichtigen 
statischen Problems bedeutet. Ihr Studium wird 
angelegentlichst empfohlen. 


Dresden. K. Beyer. 461 


rerner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 
bleibt vorbehalten): 


Dipl.-Ing. WERNER WIRZ, Die Querruder- 
wi'kung bei hoher Fluggeschwindig- 
ke't. (Sonderdruck aus der Zeitschrift „Flugwehr 
un« -Technik, Nr. 5 und 6, 1942.) 17 8. Zürich 1942, 
Buchdruckerei a. d. Sit! A.-G. Preis geh. Schweizer- 
fran' en 1,50. 


Dr.- ıg. habil. HANS PFRIEM, Nichtstatio- 
näre Wärmeübertragung in Gasen, 
insbesondere in Kolbenmaschinen. 
(VDI-Forschungsheft 413, Beilage zu „Forschung auf 
dem Gebiete des Ingenieurwesens“, Ausgabe B, 
Bd. 13, März/April 1942.) 29 S. m. 9 Bildern, Berlin 
1942, VDI-Verlag G.m.b.H. Preis brosch. 5 M. 


Prof. Dr. A. WITTING, Repetitorium und 
Aufgabensammlung zur Integralrech- 
nung. (Sammlung Göschen Bd. 147.) 121 S. m. 
52 Fig. u. 309 Beispielen. Berlin 1942, Verlag 
Walter de Gruyter & Co. Preis geb. 1,62 M. 


Prof. Dr.-Ing. ©. KRISCHER, Der Wärme- 
und Stoffaustausch im Troceknungs- 
gut, die analytische und graphische 
Behandlung der Trocknung poriger 
hygroskopischer Güter. (VDI-Forschungs- 
heft 415, Beilage zu „Forschung auf dem Gebiete 
des Ingenieurwesens“, Ausgabe B, Bd. 13, Julij 
August 1942.) 22 S. m. 18 Bildern u. 5 Zahlentafeln. 
Berlin 1942, VDI-Verlag G.m.b.H. Preis brosch. 
5M. 


Dr. WILHELM MÜLLER, Prof. a. d. Univ. Mün- 
chen, Einführung in die Mechanik des 
Fluges, 2. neubearb. Aufl. 146 S. m. 76 Abb. 
Leipzig 1942, Dr. Max Jänecke Verlagsbuchhand- 
lung. Preis brosch. 3,80 M. 


Dr. techn. ERWIN PAWELKA, Übungen aus 
der Mechanik. IV +70 S. m. 59 Abb. im Text. 
Wien 1942, Springer-Verlag. Preis brosch. 3,60 M. 


Prof. Dr. BERNHARD BAULE, Die Mathe- 
matik des Naturforschers und In- 
senieurs, Bd. 1, Differential- und 
Inteeralrechnung. VIII + 158 S. m. 161 Abb. 
Leipzig 1942, Verlag S. Hirzel. Preis kart. 7 M. 


NACHRICHTEN 


Zum 60. Geburtstag von F. Pfeiffer. 

Am 2. Mai 1943 beging Prof. Dr. F. Pfeiffer in 
Stuttgart seinen 60. Geburtstag. Prof. Pfeiffer ge- 
hört zu den Mathematikern, die in ihren Bestre- 
bungen der ZAMM ganz besonders nahestehen: 
Bei vielseitigem Interesse für die Fragen der 
reinen Mathematik gilt seine Arbeit doch in erster 
Linie den Beziehungen der Mathematik zu ihren 
Anwendungsgebieten. Diese Einstellung kommt in 
seinen Veröffentlichungen zum Ausdruck, die sich 
überwiegend mit Problemen der angewandten Ma- 
thematik und der Mechanik befassen und von 
denen hier nur sein. schöner Artikel über Elasto- 
kinetik im Handbuch der Physik genannt sei. Sie 
gibt aber. auch seiner Vorlesungstätigkeit das be- 
sondere Gepräge. Nach vorangehender Tätigkeit 
an der Techsiischen Hochschule Danzig und an 
den Universitäten Halle und Heidelberg vertritt 
Prof. Pfeiffer seit 1922 an der Technischen Hoch- 
schule Stuttgart die Mathematik. In der allge- 
meinen Vorlesung über „Höhere Mathematik“ so- 
wie in zahlreichen Spezialvorlesungen versteht er 


es, den angehenden Ingenieuren das notwendige 
mathematische Rüstzeug in einer Form zu vermit- 
teln, die zwar'stets den mathematischen Kern der 
Betrachtungen in den Vordergrund stellt, dabei 
aber nie auf die lebendige Bezugnahme zu den 
Anwendungen verzichtet. Daß von seinem unüber- 
trefflich klaren Vortrag, der Anschaulichkeit und 
mathematische Gewissenhaftigkeit auf das glück- 
lichste verbindet, auch der junge Mathematiker 
größten Gewinn hat, bezeugen zahlreiche dankbare 
Schüler. 

Seine Freunde und Schüler verehren aber in 
Prof. Pfeiffer nicht nur den Forscher und Lehrer, 
sondern auch den gütigen und hilfsbereiten Men- 
schen, der ihnen stets verständnisvoll entgegen- 
kommt und wo irgend möglich seine Förderung 
zuteil werden läßt. Mögen Frof. Pfeiffer noch viele 
Jahre rüstigen Schaffens zum Nutzen unserer 
Wissenschaft und zur Freude aller, die ihn 
kennen, beschieden sein! 


Rostock. F. Lösch. 471 
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